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Abstract 

We consider particle systems that are perturbations of the voter model and 
show that when space and time are rescaled the system converges to a solution 
of a reaction diffusion equation in dimensions d > 3. Combining this result 
with properties of the PDE, some methods arising from a low density super- 
Brownian limit theorem, and a block construction, we give general, and often 
asymptotically sharp, conditions for the existence of non-trivial stationary dis- 
tributions, and for extinction of one type. As applications, we describe the 
phase diagrams of three systems when the parameters are close to the voter 
model: (i) a stochastic spatial Lotka-Volterra model of Neuhauser and Pacala, 
(ii) a model of the evolution of cooperation of Ohtsuki, Hauert, Lieberman, 
and Nowak, and (iii) a continuous time version of the non-linear voter model 
of Molofsky, Durrett, Dushoff, Griffeath, and Levin. The first application con- 
firms a conjecture of Cox and Perkins [8] and the second confirms a conjecture 
of Ohtsuki et al [38] in the context of certain infinite graphs. An important 
feature of our general results is that they do not require the process to be 
attractive. 



AMS 2000 subject classifications. Primary 60K35. Secondary 35K57, 60J68, 60F17, 92D15, 92D40. 
Keywords and phrases. Interacting particle systems, voter model, reaction diffusion equation, evo- 
lutionary game theory, Lotka-Volterra model. 



* Supported in part by NSF Grant DMS-0505439 Part of the research was done while the author 
was visiting The University of British Columbia 

^Partially supported by NSF Grant DMS-0704996 from the probability program 

* Supported in part by an NSERC Discovery Grant 



1 



1 Introduction and Statement of Results 

We first describe a class of particle systems, called voter model perturbations in 
[7]. The state space will be {0, l}^"*, or after rescaling {0, lY^'' , where throughout 
this work we assume d > 3. The voter model part of the process will depend on 
a symmetric (i.e, p{x) = p{—x)), irreducible probability kernel p : Z'^ — t- [0, 1] with 
p{0) = 0, covariance matrix a^I, and exponentially bounded tails so that for some 
(0,1], 

p{x) < Ce-'^l^l. (1.1) 

Here and in what follows \x\ = supj \xi\. 

For 1 > e > 0, x € Z'^ and ^ € {0, 1}^'' define rescalings of p and ^ by Pe{£x) = 
p{x), and ^e(ex) = ^(x), so that G {0,1}'^^''. Also define rescaled local densities 
/f by 

/f(ex,e.) = 5^ Pe(y-ex)lU,(y) =i}, ^ = 0,1. (1.2) 

We will write fi{x,S^) if e = 1. For x,^ as above, introduce the voter fiip rates and 
rapid voter fiip rates given by 

c^(x, = [(1 - a^))fiix, + e(x)/o(x, 0], c^si^x, Ce) = e"\^{x, 0- (1.3) 

The processes of interest, G {0, l}^'*, are spin-fiip systems with rates 

c°{x,0 = c\x,0+e^4{x,0>0, (1.4) 
where c*(x,^) will be a translation invariant, signed perturbation of the form 

4{x,^) = {l-ax))hl{x,0+ax)ht>{x,0. 

Therefore the rescaled processes, (^^~2f(ex) G {0, 1}'^^'^, we will study have rates 

Ce{ex,Ce) = cUex,Ce)+4{x,0 > 0. (1.5) 

We assume there is a law q of (y^ . . . ,y^o) G Z'^^o and gf on {0, l}^o, i = 0, 1, 
and ei G (0,oo],eo G (0, 1] so that 

g! > 0, (1.6) 

and 

hlix, = -e^^hix, O+Evigmx + Y^), . . . C(x + F^"))), i = 0, 1 (1.7) 

for ah ^ G {0, l}^^ X G Z'^, e G (0, Eq] . 

Here Ey is expectation with respect to q and in practice the first term in the above 
allows us to take > 0. It is important to have gf non- negative as we will treat 
it as a rate in the construction of a dual process in Section 2. On the other hand, 
in the particular examples motivating the general theory will often be negative 
(see, e.g., (1.33) below of the Lotka-Volterra models). 
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We also suppose that (decrease k > if necessary) 

P{Y* >x)< Ce-^"" for X > (1.8) 

where Y* = max{|y-'^|, . . . |y^"|}, and there are hmiting maps gi : {0, 1}^" — t- M+ 
such that 

hm||(7f-5i||oo =0, i = 0,l. (1.9) 

£4,0 

The condition (1.7) with eq < ei (without loss of generality) easily implies the 
non-negativity in (1.5). 

We now show that the conditions (1.7)-(1.9) hold for general finite range conver- 
gent translation invariant perturbations without any non-negativity constraint on 
the gl- 

Proposition 1.1. Assume there are distinct points yi, . . . ,yNo £ 
gl.gi : {0, l}^o R such that 

hl{x,C)=gl{ax + yi),...,ax + yN,)), x G Z^e G {0, 1}"^°, (1.10) 

{x : p{x) > 0} C {yi, . . .,yNo}, 1™ H^f - 9i\\oo = z = 0, 1. 

£4,0 

Then (1.7)-(1.9) hold for appropriate non-negative gf , gi satisfying \\gf — gi\\oo = 
\\9i - 9i||oo, and = yi. 

The elementary proof is given in Section 2.1. In terms of our original rates (1.4) 
this shows that our class of models include spin-flip systems G {0,1}^'', t > 0, 
with rates 

c",{x,0 = c"{x,0 + e'^c*{x,O + e^o{e)>0, (1.11) 

where p (governing c") is now finite range, c*(x,^) = h{^{x),^{x-\-yi), . . . ,C{x-\-^]\fo)) 
is a finite range, translation invariant perturbation and o(e) means this term goes 
to zero with e uniformly in (x,.^). 

On the other hand, the random y*'s will also allow certain natural infinite range 
interactions. The formulation in terms of random locations will also simplify some 
of the arithmetic to come. 

Let ^f, t > be the unique {0, ly^ -valued Feller process with translation 
invariant flip rates given by Ce{x,£^) in (1.5) and initial state G {0,1}^^''. More 
formally (see Theorem B.3 in [33] and Section 2 of [14]) the generator of is 

the closure of 1^,(7(0 = c,(x, ^(^(r ) - 5(0), (1-12) 

on the space of g : eU^ — )■ M, depending on finitely many coordinates. 

Here is with the coordinate at x fiipped to 1 — ^(x). The condition (B4) of 
Theorem B.3 in [33] is trivial to derive from (1.7). 

We stress that conditions (1.5)-(1.9) are in force throughout this vi^ork, 
and call such a process a voter model perturbation. 
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Given a process taking values in {0, 1}^'*, or more generally in {0, l}^^'', we say 
that coexistence holds if there is a stationary distribution with 

^ ( E ^(^) = E 1 - ^(^) = ) = 1- (1-13) 

For voter model perturbations it is easy to see this is equivalent to both types being 
present z^-a.s.-see Lemma 6.1 at the beginning of Section 6. 
We say the i 's take over if for all L, 

P{^t{x) = i for all x G [—L, L]'^ for t large enough) = 1 (1-14) 

whenever the initial configuration has infinitely many sites in state i. 

Our main results, Theorems 1.15 and 1.16 in Section 1.6 below, give (often sharp) 
conditions under which coexistence holds or one type takes over, respectively, in a 
voter model perturbation for small enough e. Of course these results then hold 
immediately for our originally unsealed processes, again for small enough e. 

1.1 Hydrodynamic limit 

As d > 3, we see from Theorem V.1.8 of [32] the voter model with flip rates 
c'"{x,(^) = Cj'(x,^) has a one-parameter family of translation invariant extremal 
invariant distributions {P„ : u G [0,1]} on {0,1}^'* such that E„(^(2;)) = u. We 
write {g)u for Eu{g{0)- (1-7) and (1.9) imply 

lim ||/if — /ijIIoo = where (1-15) 

e4-0 

h,ix, = -efh{x, + E{gMx + Y'),...,^{x + y^°))). 

Define 

f{u) = ((1 - m)hi{0,0 - mho{0,0)u- (1.16) 

Then / is a polynomial of degree at most A'^o + 1 (see (1.26) and Section 1.8 below). 
The non- negativity condition (1.5), the fact that 

(c^(0,0)o = (c^r(0,0)i=0, (1.17) 

and the convergence (1.15) show that 

/(0)>0, /(1)<0. (1.18) 

Our first goal is to show that under suitable assumptions on the initial conditions, 
as e — 7- the particle systems converges to the PDE 

— = -Au + f{u), u{0,.)=v{-), (1.19) 

The remark after Proposition 2.1 in [2] implies that for any continuous v : ^ [0, 1] 
the equation has a unique solution u, which necessarily takes values in [0, 1]. 
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For a continuous v as above we will say that a family of probability measures 
{Xe} of laws on {0, 1}^^ has local density v if the following holds: 

There is an r G (0,1) such that if = \e''-^]e, Qe = [0, ae)^neZ"', \Qs\ = card(Qe), 
and 

D{x, e) = ^ ^(^ + y) for X G a,Z^ ^ G {0, 1}^^' , (1.20) 
then for ah R,5 > 0, 

lim sup A£(|D(x,0 - ^^(x)! > -5) = 0. (1.21) 

\x\<R 

The family of Bernoulli product measures given by 

n 

X^{^{wi) = l,i = l,...,n) = Yl v{wi) for all n G N and Wi G eZ'^ . (1.22) 

i=l 

certainly satisfies (1-21) for all r G (0, 1). 

Theorem 1.2. Assume v : ^ [0, 1] is continuous, and the collection of initial 
conditions {£,q} have laws {A^} with local density v. Let x^ G M"' and x^ G eZ'^, 
k = 1, . . . K satisfy 

x^ — )• x^ and e~^\x^ — | — )• oo as e — )• for any k / k' . (1.23) 

Ifu is the solution of (1.19), then for any rj G {0, i}{iv,i}x{i,...A'}^ yi, . . . , y/, G Z'' 
and r > 0, 

lim PieAxs + eyi) = Vi,k, i = l,...,L, k = l,...K) 

e-5>0 

= n^i{^(y^) = ^^.'^'^ = 1' • • • '^})«(T,x^)- (1-24) 

fc=l 

/n particular, if x^ G eZ"^ satisfies — t- x as e — t- 0, t/ien 

lim P(CT(xe) = 1) = n(r, x) /or allT >0,x e R'^. (1.25) 

De Masi, Ferrari and Lebowitz [11], Durrett and Neuhauser [20] and Durrett 
[15] have proved similar results for particle systems with rapid stirring. The local 
equilibrium for rapid stirring is a Bernoulli product measure, but in our setting it 
is the voter equilibrium. As a result there is now dependence between nearby sites 
on the microscopic scale. However, there is asymptotic independence between sites 
with infinite separation on the microscopic scale. 

It is easy to carry out a variance calculation to improve Theorem 1.2 to the 
following L^-convergence theorem (see the end of Section 3 for the proof). If 5 > 
and X G M'^, let Is{x) be the unique semi-open cube Y[i=i[^i^Af^i + ^)^): ^ 
which contains x. 
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Theorem 1.3. Assume the hypotheses of Theorem 1.2. 

Assume 6{£) G eN decreases to zero so that (5(e)/e — t- oo as e J, 0. If 

u%t,x)= Yl etiy){e/S{e)f, 

then as e ^ 0, u'^{t, x) — t- u{t, x) in LP' uniformly for x in compacts, for all t > 0. 

A low density version of this theorem, in which the hmit is random (super-Brownian 
motion with drift), was proved in [7] and is discussed in Section 1.8. 

To apply Theorem 1.2 to the voter perturbation we will have to evaluate f{u). 
This is in principle straightforward thanks to the duality between the voter model 
and coalescing random walk which we now recall. Let {B^ : x G Z'^} denote a rate 
1 coalescing random walk system on with step distribution p and Bq = x. For 
A,B C Z'^, let 1/^ = {Bf -.xG A}, t{A) = mf{t : = 1} and t{A, B) be the first 
time it'^if + (it is oo if either ^ or is empty). The duality between B and the 
voter model (see (V.1.7) and Theorem V.1.8 in [32]) implies for finite A^B (Z Z*^, 

|A| \B\ 

= j;j;^^(l-n)^P(||^|=j,||«| = fc,T(Ai?) = oo). (1.26) 

j,=0 A:=0 

The = term is non-zero only if i? = in which case the above probability is 
P(|,^^| = j), and similarly for the j = term. It follows from (1.26) and the form 
of the perturbation in (1.7) that f{u) is a polynomial of degree at most A'^o + 1 with 
coefficients given by certain coalescing probabilities of B (see (1.86) below). 

1.2 PDE results 

As in Durrett and Neuhauser [20], Theorem 1.2, in combination with results for the 
PDE and a block construction, lead to theorems about the particle system. Durrett 
[16] surveys results that have been proved by this method in the last 15 years. To 
carry out this program we will also need some low density methods taken from the 
superprocess limit theorems of Cox and Perkins [7, n]. 

To prepare for the discussion of the examples, we will state the PDE results on 
which their analysis will be based. The reaction function / : M — )• M is a continuously 
differentiable function (as already noted, in our context it will be a polynomial). 
Assume now, as will be the case in the examples, that /(O) = /(I) = 0. We let u(t, rr) 
denote the unique solution of (1-19) with continuous initial data u : M*^ — t- [0, 1]. 

We start with a modification of a result of Aronson and Weinberger [2] . 

Proposition 1.4. Suppose /(O) = f{a) = 0, /'(O) > 0, /'(a) < and f{u) > for 
u G (0, a) with < a < 1. There is a w > so that if the initial condition v is not 
identically 0, then 

liminf inf u(t,x) > a. 

t^oo \x\<2wt 
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We also will need an exponential rate of convergence in this case under a stronger 
condition on the initial conditon. We formulate it for / < on (0, 1). The brief 
proofs of Propositions 1.4 and 1.5 are given at the beginning of Section 4. 

Proposition 1.5. Assume / < on (0, 1) and /'(O) < 0. There is a w > 0, and 
if 6 > there are positive constants Ls, c = as, and C = Cs so that if L > Ls and 
v{x) < 1 — 5 for \x\ < L, then 

u{t, x) < Ce"^* for |x| < L + 2wt. 

There are different cases depending on the number of solutions of f{u) = in 
(0, 1). In all cases, we suppose that /'(O) ^ and /'(I) ^ 0. 

Case I: f has zero roots in (0, 1). In this case we can apply Propositions 1.4 (with 
a = 1) and 1.5, and their obvious analogues for — /. 

Case II: f has one root p G (0, 1). There are two possibilities here. 

(i) /'(O) > and /'(I) < and so the interior fixed point p G (0, 1) is attracting. 
In this case we will also assume f'{p) ^ 0. Then two applications of Proposition 1.4 
show that if f ^ and v ^ 1 

lim sup \u{t,x) — p\ = 0. (1-27) 

\x\<wt 

(ii) /'(O) < and /'(I) > 0, so that and 1 are locally attracting and p G (0, 1) 
is unstable. In this case the limiting behavior of the PDE is determined by the speed 
r of the traveling wave solutions, i.e., functions w with w^—oo) = p and w{oo) = 
so that u(t, x) = w{x — rt) solves the PDE. The next result was first proved in d = 1 
by Fife and McLeod [25]. See page 296 and the appendix of [20] for the extension to 
d> 1 stated below as Proposition 1.6. The assumption there on the non-degeneracy 
of the interior zeros are not necessary (see Fife and McLeod [25]). These references 
also show that ^ 

sgn(r)=sgn^y f{u)duj. (1-28) 

|x|2 will denote the Euclidean norm of x and the conditions of Case II (ii) will apply 
in the next two propositions. 

Proposition 1.6. Suppose Jq f{u)du < and fix rj > 0. If 6 > there are positive 
constants L^, cq = co{6), and Co = Co{6) so that if L > and v{x) < p — 6 when 
\x\2 < L, then 

u{t,x) < Ce""* for \x\2 < (|r| - r])t. 

For the block construction it is useful to have a version of the last result for the L°° 
norm, and which adds an L to the region in which the result is valid. 

Proposition 1.7. Suppose Jq f{u)du < 0. There is a w > 0, and if S > there are 
positive constants Lg, c = cs and C = Cs so that if L > and v{x) < p — 5 for 
\x\ < L, then 

u{t, x) < Ce"^* for |x| < L + 2wt. 
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The short derivation of Proposition 1.7 from Proposition 1.6 is given at the beginning 
of Section 4. 

In the next three subsections we first motivate our main results on coexistence 
and domination by a single type (Theorems 1.15 and 1.16) by illustrating their use 
in three distinct families of examples. 

1.3 Lotka-Volterra systems 

We now apply Theorem 1.2 to a stochastic spatial Lotka-Volterra model introduced 
by Neuhauser and Pacala [-37] . In addition to the kernel p for the voter model, the 
flip rates depend on two non- negative competition parameters, ao and ai, and are 
given by 

CLv{x,0 =/i(/o + ao/i)(l - + /o(/i + ai/o)e(x) 

=c"{x,0 + (ao - l)/i'(l -e(^)) + (ai - l)/oe(^)- (1-29) 

In words, a plant of type i at x dies with rate fi{x, ^)-|-Qj/i_j(x, ^) and is immediately 
replaced by the type of a randomly chosen neighboring plant, which will be 1 — i 
with probability fi-i{x,S,). The death rate reflects the effects of competition from 
neighboring sites. The constant Oj represents the effect of competition on a type 
i individual from neighbors of the opposite type. If aj < for both i = 0, 1, then 
individuals prefer to be surrounded by the opposite type and ecological arguments 
suggest coexistence will hold. Conversely if both > 1, we expect no coexistence. 

We refer to the associated Feller process {Ctit > 0) as the LV{ao, ai) process. 
Proposition 8.1 of [8] implies that 

if ao Aai > 1/2 then LV{ao,Oii) is monotone (or attractive). (1.30) 

Write LV{a) < LV{a') if LV{a') stochastically dominates LV{a), that is, if one 
can define these processes, ^ and respectively, on a common probability space so 
that ^ < pointwise a.s. Then, as should be obvious from the above interpretation 
of ai (see (1.3) of [8]), 

< a'o — ^0,0 < ai < a'l, and either a^ A ai > 1/2 (1-31) 
or Oq a a'l > 1/2, implies LV{a) < LV{a). 

If 6i € M, let ai = ai = l + e^Oi and consider the rescaled Lotka-Volterra process 

^tix) = 4-2t(e"^a;), x G eZ'^. (1.32) 

From (1.29) we see that this process has rates given by (1.5) with 

h|{x,O = h^{x,O = 0l^ifi{x,O^i = 0,l, xGZ^Cg{0,1}^'. (1.33) 

To verify (1.7) take < ei < (6*0 )"^/^ A ((9^)"^/^ A^o = 2, Y\Y'^ chosen indepen- 
dently according to p, and define for i = 0, 1, 

glivum) = 9i{m,m) = ^fmi^ - m) + i^f + Oi^,)i{r]i = r]2 = i) > o. (i.34) 
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Then 

as required. Therefore is a voter model perturbation. 

To calculate the limiting reaction function in this case consider the system of 
coalescing random walks {B^ : x G Z"^} used in the duality formula (1.26). Let 
{61,62} be i.i.d. with law p{-) and independent of the {B^ : x G Z'^}. If we abuse 
our earlier notation and let denote expectation on the product space where 
(ei, 62) and are independent, and is given the voter equilibrium with density u, 
then from (1.16), (1.33) and the fact that fi{0,£,f = PeiCi^i) = '^(^2) = i), we have 

/(n) = eoiil - mnei)ae2))u - Oiim^ - e(ei))(l - C{e2)))u- 

In view of (1.26) we will be interested in various coalescence probabilities. For 
example, 

p{x\y, z) = P{3t > such that Bf = B^, and Vt > 0, B^ / B^ and B^ / Bf) 
and 

p[x\y\z) = P{Bf,B^ and B^ are all distinct for all t). 

In general walks within the same group coalesce and those separated by at least one 
bar do not. If we define 

P2 =p(0| 61,62), P3 =p(0|ei|e2), (1.35) 

where the expected value is taken over ei, 62, then by the above formula for / and 
(1.26), 

/(n) = 6*0^(1 - u)p2 + Oou^{l - u)p3 - 6*1(1 - u)up2 - 6*1(1 - ufup3 

= U{1 - u)[eoP2 - Ol{p2+P3) + upsiOo + 61)]. (1.36) 

To see what this might say about the Lotka-Volterra model introduce 

" ^''"'^ = P3(^l+^0) ^'-''^ 

SO that f{u) = for u = 0, 1 or 'u*(6li/6lo). If 

mo = (1.38) 

then u*{m) increases from to 1 as m increases from mg to tuq^. We slightly abuse 
the notation and write u* for u*{0i/9q). 

To analyze the limiting PDE we decompose the ^0 ~ ^1 plane into 5 open sectors 
drawn in Figure 1 on which the above O's are all simple. 
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Figure 1: Phase diagram near (1,1) for the Lotka-Volterra model with the shape of 
/ in the regions. 
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• If G / > on (0, ti*), / < on {u* , 1), so u* £ (0, 1) is an attracting fixed 
point for the ODE. Then (1.27) shows the PDE solutions will converge to u* 
given a non-trivial initial condition in [0, 1]. As a result we expect coexistence 
in the particle system. 

• If ^ G i?2, / < on (0, 1), is an attracting fixed point for the ODE. Propo- 
sition 1.4 implies solutions of the PDE will converge to given a non-trivial 
initial condition and we expect O's to win. 

• If G i?3, / > on (0, 1), 1 is an attracting fixed point for the ODE and so 
by the reasoning from the previous case we expect I's to win. 

• On i?4Ui?5, u* G (0, 1) is an unstable fixed point, while and 1 are attracting 
fixed points for the ODE. This is case 2 of Durrett and Levin [19], so we expect 
the winner of the competition to be predicted by the direction of movement 
of the speed of the decreasing traveling wave solution u{x,t) = w{x — rt) with 
w{—oo) = 1 and u;(oo) = 0. If r > then I's will win and if r < then O's will 
win. Symmetry dictates that the speed is when 6q = 6i, so this gives the 
dividing line between the two cases and the monotonicity from (1.31) predicts 
O's win on R4 while I's win on R^. Alternatively, by (1.28) r has the same 
sign as f{u) du which is positive in and negative in i?4. 

Our next two results confirm these predictions for a close to (1,1). For < r/ < 1, 
define regions that are versions of i?i , i?2 U R4, and R^ U i?5 shrunken by changing 
the slopes of the boundary lines: 

= I (ao, ai) e [0, 1]^ : < a^ - 1 < "^"("""^H , 

L mo 1 — rj J 

Aq = |(ao, ai) G (0, 00)^ : < ao < 1, mo(l — r/)(ao — 1) < ai — 1, 

or 1 < ao, (1 + r]){ao - 1) < ai - l|, 

A^ = |(ao,ai) G (0,00)^ : < ao < 1, ai - 1< "° ~ ^ , , 
L mo(l-7?) 

or 1 < ao, ai — 1 < (1 — r?)(ao — 1)|. 

Theorem 1.8. For < t] < 1 there is an ro(?7) > 0, non- decreasing in r], so that 
for the LV{a): 

(i) Coexistence holds for (ao,ai) G and 1 — ao < ^^{7]). 

(a) 7/(ao,ai) is as in (i) and Va is a stationary distribution of with i/a{^ = or ^ = 
1) = 0, then 

r.,(C(x) = l)-n*(^^) 
V ao — 1 / 



sup 



< 



(i) is a consequence of Theorem 4 of [8], which also applies to more general 
perturbations. The main conditions of that result translate into /'(O) > r] and 
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Figure 2: Coexistence on C, type i takes over on Aj. 



/'(I) > ?7 in our present setting (see (1.88) in Section 1.8 below), (ii) sharpens (i) 
by showing that if r] is small then the density of I's in any nontrivial stationary 
distribution is close to the prediction of mean-field theory. Durrett and Neuhauser 
[20] prove results of this type for some systems with fast stirring and /(I) < 0. 
Neuhauser and Pacala [37] conjectured that coexistence holds for ah = ai < 1 
(see Conjecture 1 of that paper) and proved it for Oi sufficiently smah. Hence (i) 
provides further evidence for the general conjecture. 

The next result is our main contribution to the understanding of Lotka-Volterra 
models. It shows that (i) of the previous result is asymptotically sharp, and verifies 
a conjecture in [8] (after Theorem 4 in that work). We assume p has finite support 
but believe this condition is not needed. 

Theorem 1.9. Assume p{-) has finite support. For < t] < 1 there is anrQ{rj) > 0, 
non- decreasing in r], so that for the LV{a): 

(i) O's take over for (Qo,ai) G Aq and < |ao — 1| < ^^{rj), 

(ii) I's take over for (ao,ai) G A^ and < |ao — 1| < ^^{rj). 

Conjecture 2 of [37] states that I's take over for a G A]', ao > 1 and O's take 
over for a G Ag, qq > 1. Theorem 1.9 establishes this result asymptotically as a 
gets close to (1,1) at least for d > 3. 

Together, Theorems 1.8 and 1.9 give a fairly complete description of the phase 
diagram of the Lotka-Volterra model near the voter model. In Figure 2 C is the union 
over r] G (0, 1) of the regions in Theorem 1.8 (i) on which there is coexistence, and Aj, 
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i = 0, 1, is the union over r] of the regions in Theorem 1.9 (i) and (ii), respectively, on 
which i's take over, as weh as other parameter values for which the same result holds 
by monotonicity. For example, if (0:0,01) S Ai, with qq A ai > 1/2, and {aQ,a[) 
has Oq > cto and a[ < ai, then by (1.31), (q:q,q;'^) E Ai. Theorem 1.8 (i) and 
Theorem 1.9 show that the three mutually exclusive classifications of coexistence, 
O's take over, and I's take over, occur on the three regions, C, Aq and Ai, meeting 
at (1, 1) along mutually tangential lines with slopes mg, 1 and m^^. 

1.4 Evolution of cooperation 

Ohtsuki et al [oS] considered a system in which each site of a large (A^ vertex) graph 
G is occupied by a cooperator (1) or a defector (0). Simplifying their setting a bit, 
we will assume that each vertex in G has k neighbors. The interaction between 
these two types is governed by a payoff matrix with real entries 



This means that a cooperator receives a payoff a from each neighboring cooperator 
and a payoff f3 from each neighboring defector, while for defectors the payoffs are 7 
and 6 from each neighboring cooperator or defector, respectively. The terminology 
is motivated by the particular case in which each cooperator pays a benefit 6 > to 
each neighbor at a cost c > per neighbor, while each defector accepts the benefit 
but pays no cost. The resulting payoff matrix is then 



In this case the payoff for D always exceeds that for C irregardless of the state 
making the payoff. As a result in a homogeneously mixing population cooperators 
will die out. The fact that such cooperative behavior may nonetheless take over in a 
spatial competition is the reason for interest in these kind of models in evolutionary 
game theory. In a spatial setting the intuition is that it may be possible to the C"s 
to form cooperatives which collectively have a selective advantage. 

If ni{y) is the number of neighboring i's for site y £ G, i = 0,1, and (,{y) G {0, 1} 
is the state at site y, then the fitness Pi{y) of site y in state i is determined by its 
local payoffs through 



Here w S [0, 1] is a parameter determining the selection strength. Clearly for some 
wo{a, (3,^, 6, k) > 0, Pi > for w G [0,i/;o], which we assume in what follows. For 
the death-birth dynamics in [.■iS] a randomly chosen individual is eliminated at x 
and its neighbors compete for the vacated site with success proportional to their 



C D 

C a /3 
D 7 5 




(1.39) 



Pi{y) = l-w + w{ani{y) + f3no{y)) if ^(y) = 1 
Po{y) = l-w + w{-fni{y) + (5no(y)) if i{y) = 0. 



(1.40) 
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fitness. We consider the continuous time analogue which is the spin-flip system 
(,tix) € {0, 1}, X £ G, with rates (write y ~ x if and only if y and x are neighbors) 

c{x,0 = {'^ - ^{x))ri{x,0 + ^(.x)ro{x,0, 

n{x,i) = — ' . ^^TT -^G 0,1. (1.41 

In [-■JS] the authors use a non-rigorous pair approximation and diffusion approx- 
imation to argue that for the cooperator-defector model in (1.39), for large popula- 
tion size N and small selection w, cooperators are "favored" if and only if &/c > k. 
Here "favored" means that starting with a single cooperator the probability that 
cooperators take over is greater than 1 /N ^ the corresponding probability in a selec- 
tively neutral model. They also carried out a number of simulations which showed 
reasonable agreement foi N ^ k although they noted that b/c > k appeared to be 
necessary but not sufficient in general. It is instructive for the reader to consider the 
nearest neighbor case on Z starting with cooperators to the right of and defectors 
to the left. It is then easy to check that the C/D interface will drift to the left, and 
so cooperators take over, if and only if b/c > 2. This was noted in [38] as further 
evidence for their b/c> k rule. 

Our main result here (Corollary 1.12 below) is a rigorous verification of the 
b/c > k rule for general symmetric translation invariant graphs with vertex set Z"^ 
when w is small. More precisely, choose a symmetric (about 0) set M of neighbors 
of of size k, not containing 0, and consider the graph with vertex set Z"^ and x ^ y 
if and only \i x — y € A/". Assume also that the additive group generated by J\f is Z"^ 
and Ylix&N^i^il^ ~ ^'^^ij^ ^o that p{x) = k~^l{x G M) satisfies the conditions on 
our kernel given in, and prior to, (1.1). Set w = e^. For x G Z*^ and ^ G {0, l}^'', let 

/f (x, i) = k~' my) = ^)Mx, e) = E E = ^(^) = ^) ^ 

e, (x, = (/3A: - l)/i (x, C) + A:(a - /3)/f ^ (x, C) , 
0o(x, = ilk - l)/o(x, + k{6 - j)fS'\x, e), 
(l){x,0 = {Oo + Oi){x,0. 
Using (1.40) in (1.41), we get 

r^{x,0 = ^^^^{x,0■ (1-42) 
1 + e^fo 



Note that 



and 



V |6'o| V \(p\{x,C) < 2k{l + \a\ + \/3\ + \-f\ + \5\) = R, (1.43) 
1 + e' 



2, 



f + e\e- fci>) + eV(/<A - 0) [E(-^'"'' 
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It follows that Q{£x) = C£-2t(x), X G Z*^ has rates given by (1.41)) has spin-flip 
rates given by (1.5) with 

/if (x, = diix, i) - Mx, 04>{x, + e^Aihix, i)A{x, 0). (1.44) 

If 

< (2i?)-\ (1.45) 

then one easily from (1.43) that 

mUiA,Oi){x,m<mR+^)- a-^e) 

From this and (1-44) it is clear that the hypotheses of Proposition 1.1 hold with 

\\9t-gi\\oo<e'2R{R+l) (1.47) 

and 

hi{x,o = {ei-M){x,o, i = 0A. 

Hence is a voter model perturbation. This also implies 

ho + hi = 00 + 01- (I) = 0. (1.48) 
Some elementary arithmetic and (1.48) lead to 

ho{x,0 = (7 - mfofiix,0 + k{5 - ^)f^^\x,0 

- A:/o(x,e)[(a - /3)/f ^ + {5- ^)f^\x,0, 
hi{x,0 = -ho{x,0. (1.49) 

As before, let 61,62,63 denote i.i.d. random variables with law p. If Pg denotes 
averaging over the e^'s then we have 

/,(0,e) = Pe(e(ei) = i), fP{0,0 = Pe{aei) = h^{ei + e2) = i), (1.50) 

/*i(0,O/f (0,6 = Peiaei) = ii,C{e2) = ^2,^(62 + 63) = ^2), 

and similarly for higher order probabilities. We also continue to let (•)„ denote 
expectation on the product space where (61 , 62 , 63 ) and the voter equilibrium ^ are 
independent. If ^ = 1 — ^, then starting with (1.16) we have, 

^ =k~\mh{o,o - mho{o,c))u = k-\hi{o,^))^, 

where in the last equality (1.49) is used to see that what appears to be a quartic 
polynomial is actually a cubic. Using (1.50) and some arithmetic we obtain 

^ =(/3 - 7)(e(ei)e(e2))« + (a - /3)(|(6i)e(62)e(e2 + 63))^ (1.51) 

+ (7-5)(e(ei)e(e2)|(e2 + 63))„. 

To simplify further we will use a simple lemma for coalescing random walk 
probabilities (Lemma 1.17 in Section 1.7 below) together with the duality formula 
(1.26) to establish the following more explicit expression for / in Section 1.7. 
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Lemma 1.10. 1^ = [{^ - 5) + k-^-y - 6)]p{0\ei)u{l - u) 

+ [(a-/3)-(7-5)]['u(l-n)(p(ei|e2,e2+e3)+up(ei|e2|e2+e3))]. 

Rather than try to analyze this cubic as in Section 1.2, assume a — (3 = ^ — 6 
(which holds in our motivating example) so that / becomes a quadratic with roots 
at and 1. li (3 — 5 > k~^{5 — 7), then / is strictly positive on (0,1) and so 
Proposition 1.4 shows the PDE solutions will converge to 1. If /3 — (5 < k~^{5 — 7), 
then / is strictly negative on (0,1) and so by symmetry the PDE solutions will 
converge to 0. As a result for w = small, in the former case we expect I's to take 
over and in the latter case we expect O's to take over, and this is in fact the case. 
The following result is proved in Section 1.7. 

Theorem 1.11. Consider the spin- flip system on (d > 3) with rates given by 
(1.41) where a — fi = ^ — 5. If 'y — 6 > k{6 — (3), then 1 's take over for w > 
sufficiently small; ifj — 6 < k{6 — 13), then O's take over for vu > sufficiently small. 

The particular instance of (1.39) follows as a special case. 

Corollary 1.12. Consider the spin-flip system onU^ (d > 3) with rates given by 
(1.41) where the payoff matrix is given by (1.39). Ifb/c > k, then the cooperators 
take over for w > sufficiently small, and ifb/c < k, then the defectors take over 
for w > sufficiently small. 

Proof. In this case a — (3 = ^ — 5 = b, 5 — f3 = c, and so we have 
7 - 5 > A;((5 - /?) iff 6 > /cc iff b/c > k. 

□ 

1.5 Nonlinear voter models 

Molofsky et al. [3G] considered a discrete time particle system on in which each 
site is in state or 1 and 

P{^n+i{x,y) = l\Cn) =Pk 

if k of the sites (x, y), {x + 1, y), [x — 1, y), [x, y + 1), (x, y — 1) are in state 1. They 
assumed that po = and = 1, so that all O's and all I's were absorbing states and 
= 1 — P4 and p2 = 1 — P3 1 so that the model was symmetric under interchange of 
O's and I's. If the states of adjacent sites were independent then the density would 
evolve according to the mean field dynamics 

xt+i = h(xt) = pi ■ 5xt{l - xt)^ +P2 • 10xj(l - xtf 

+ (1 - P2) ■ 10x?(l - xtf + (1 - pi) ■ 5x^(1 - xt) + xl 

Based on simulations and an analysis of the mean-field equation, Molofsky et 
al [36] predicted the phase diagram given in Figure 3. To explain this, h{x) = x is a 
fifth degree equation with 0, 1/2, and 1 as roots. h'{0) = h'{l) = 5pi so and 1 are 
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Figure 3: Conjectured phase diagram for the discrete time two-dimensional nonlinear 
voter model of [36]. 



locally attracting if 5j3i < 1 and unstable if 5pi > 1. h'{l/2) = (15 — 15pi — 10p2)/8) 
so 1/2 is locally attracting if 15pi + 10p2 > 7 and unstable if 15pi + 10p2 < 7. 
From the stability properties of 0, 1/2, and 1, it is easy to determine when there 
are additional roots a and 1 — a in the unit interval and whether or not they are 
stable. The four shapes are given in Figure 4. To make the drawing easier we have 
represented the quintic as a piecewise linear function. 

The implications of the shape of f{u) (= h{u) — ti in the above) for the behavior 
for the system will be discussed below in the context of a similar system in continuous 
time. There we will see that the division between 4A and 4B is dictated by the speed 
of traveling waves for the PDE. Here we have drawn the "Levin line" 6pi + 2p2 = 2 
which comes from computing the expected number of I's at time 1 when we have 
two adjacent I's at time 0. Simulations suggest that the true boundary curve exits 
the square at (0.024, 1), see page 280 in [36]. 

For our continuous time model the perturbation rate from the voter model is 
determined by four points chosen at random from x+N where N is the set of integer 
lattice points in {[—L, — {0}). Let a{i) > be the flip rate at a given site when 
i (randomly chosen) neighbors have a type disagreeing with that of the site and 
suppose a(0) = 0. Let (Yi, . . . 14) be chosen at random and without replacement 
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Figure 4: Four possible shapes of the symmetric quintic /. Black dots indicate the 
locations of stable fixed points. 



from M . Then we consider the spin-flip system £^{x),x G with rates 

c(x, = c^(x, + e\l - i{x))EY{gmx + Y^), . . . ,i{x + Y^))) 
+ C{x)EY{go{C{x + ¥'),..., C{x + Y')))], 

where 

4 4 

1 1 

Then the rescaled system ^f{£x) = ^f^-2{x), x G Z'^ is a voter model perturbation 
since the required conditions are trivial (clearly (1.7) holds with ei = oo). We call ^ 
the nonlinear voter model. General models of this type were introduced and studied 
in [4]. 

We abuse our notation as before and incorporate expectation with respect to an 
independent copy of y = {Y^, . . . , Y^) in our voter equilibrium expectation (•)„. If 
Y^ = 0, then a short calculation shows that our reaction function in (1.16) is now 

4 

f{u) = ^a{j){q,{u) - q,{l - n)), (1.52) 



where 



Clearly 



/4\ ^ 

^•^^ 1=0 i=5-j 



/(O) = /(I) = /(1/2) = and f{u) = -/(I - u). (1.53) 



It does not seem easy to calculate / explicitly, but if L is large, most of the sum comes 
from y* that are well separated and so the ^ values at the above sites should be nearly 
independent. To make this precise let A = {Y^~^ , . . . ,y^}, B = {Y^, . . . ,Y^~^}, 
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(1 < i < 4), and note by (1.26) that 
qj{u) : 



^) E E ^'(1 - = ^' I?- 1 = ^' ^(^' B) = oo) 

i=i k=i 



J 
4 
J 

where 

+ E E + ^ < 5)^^(1 - u)' Filial = i, ICll = k, t{A, B) = oo] 

i=l fc=l 

5 

= ^di(j,L)'u*. 

i=l 

If riQ{L) = P{\i^^\ < 5), that is the probabihty that there is a coalescence among 
the random walks starting at Y^, . . . , y^, then it follows easily from the above that 
\di{j,L)\ < co7]q{L). Use this in (1.52) to conclude that f{u) = fi{u) + f2{u), where 
/2 includes the (smaller) contributions from the qj^s. That is 

5 

f2{u) = ^e(j,L)n^ 

where 

sup \e{j,L)\ < cir/o(L), (1.54) 

i<i<5 

and 

/i(n) = - u[a(4)(l - + a(3) • 4^(1 - uf + a(2) • 6u'^{l - uf + a(l) • 4u^{l - u)] 
+ (1 - u)[a(4)n^ + a(3) • 4u^{l - u) + a(2) • 6^2(1 - uf + a(l) • 4ii(l - uf] 
=6iu(l - n)^ + 62^2(1 - uf - b2U^{l - uf - biu^{l - u) 

where bi = 4a(l) — a(4) and 62 = 6a(2) — 4a(3). By symmetry we have 

/i(0) = /i(l) = /i(l/2) = and /i(n) = -/i(l - n). (1.55) 

Clearly r]Q{L) — )• as L — )• 00, in fact well-known return estimates (such as Lemma 2.6(a) 
below with to = 0, rg = 1 and p large) and a simple optimization argument show 
that 

r?o(L) < C5L-['^('^-2)/(2('^-i))]+5, 5 > 0. (1.56) 
To prepare for the next analysis we note that 

f[{u) =bi[{l - uf - 4n(l - uf] + 62[2n(l - uf - 3u\l - uf] 
- b2[3u^{l - uf - 2n^(l - u)] - 6i[4n^(l - u) - u'^], 
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Figure 5: Phase diagram for the continuous time nonhnear voter model with large 
range in d > 3. 



and so we have /{(O) = /{(I) = bi and /{(1/2) = -(661 + 262)/16. A little calculus, 
left for the reader, shows 

r^'h{u)du = = - f h{u)du. (1.57) 

JQ -Ly2 7i/2 

We are now ready to describe the phase diagram for the nonlinear voter. Consult 
Figure 5 for a picture. Note that in what follows when L is chosen large, it is 
understood that how large depends on a = (a(l), . . . , a(4)). 

(1) /((O) > 0, /{(1/2) < and so by (1.54) for large enough L the same is true for 
/. In this case 0, 1/2, and 1 are the only roots of / (all simple) and 1/2 is an 
attracting fixed point for the ODE. An application of Proposition 1.4 on [0, 1/2] 
and a comparison principle, showing that solutions depend monotonically on 
their initial data (see Proposition 2.1 of [2]), to reduce to the case where 
V E [0, 1/2], shows that any non-trivial solution u of the PDE (1.19) satisfies 
liminft_5.oo inf|j;|<2«,t 'u(i, x) > 1/2 for some w; > 0. The same reasoning with 
and 1 reversed shows the corresponding upper bound of 1/2. Therefore any 
non-trivial solution of (1-19) will converge to 1/2 and we expect coexistence. 

(2) /{(O) > 0, /((1/2) > and so by (1.54) for large enough L the same is 
true for /. In this case 0, 1/2, 1 are unstable fixed points for the ODE 
and there are attracting fixed points for the ODE at a and 1 — a for some 
a G (0, 1/2). All are simple zeros of /. Another double application of Proposi- 
tion 1.4 now shows that any non-trivial solution u{t,x) to the PDE will have 
liminft^ooinf|^|<2t„t^x(i,x) > a and limsupt^o^ sup|^|<2^„t'u(t,x) < 1 - o, so 
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we expect coexistence. Simulations in [3G], see Figure 7 and the discussion on 
page 278 of that work, suggest that in this case there may be two nontrivial 
stationary distributions: one with density near a and the other with density 
near 1 — a. The symmetry of / about 1/2 is essential for this last possibility 
as we note below (see Theorem 1.14). 

(3) /((O) < 0, /((1/2) > and so by (1-54) for large enough L the same is true 
for /. In this case 0, 1/2, and 1 are the only roots of / (all simple) and 
1/2 is an unstable fixed point while and 1 are attracting. In this bistable 
case the winner is dictated by the sign of the speed of the traveling wave, but 
by symmetry (recall (1.53)) the speed is 0. One would guess that clustering 
occurs in this case and there are only trivial stationary distributions, but our 
method yields no result. 

(4) /((O) < 0, /((1/2) < and so by (1-54) for large enough L the same is true 
for /. In this case 0, 1/2, 1 are attracting fixed points and there are unstable 
fixed points at a and 1 — a for some a € (0, 1/2) (all simple zeros of /). By 
the discussion in Case II in Section 1.2 (with [0, 1/2] and [1/2, 1] in place of 
the unit interval) there are traveling wave solutions Wi{x — Cit), i = 1,2 with 
wi{—oo) = 1, wi{oo) = W2{—oo) = 1/2 and u;2(oo) = 0. Symmetry implies 
C2 = — ci, but we can have ci < < C2 (Case 4A) in which case Proposition 1.7 
and its mirror image show that solutions to the PDE will converge to 1/2 
providing that the initial condition is bounded away from and 1 on a large 
enough set. We again use the comparison principle as in Case 1 to assume 
the initial data takes values in the appropriate interval, [0, 1/2] or [1/2, 1], and 
assume L is large enough so that the integrals of /i and / on [0, 1 /2] (and 
hence on [1/2,1]) have the same sign. Hence we expect coexistence in Case 
4A and all invariant distributions to have density near 1/2. If ci > > C2 
(Case 4B) and L is large enough, there is a standing wave solution wo{x) of 
the PDE in d = 1 with 'Wo{—oo) = 0, wo{oo) = 1 (see p. 284 in [26]), and our 
method yields no result. 

Theorem 1.13. Assume (61,62) are o^s in Case 1, 2 or 4 A. If L is sufficiently large 
(depending on a) then: 

(a) Coexistence holds for e small enough (depending on L and a). 

(b) In Case 1 or 4 A if r] > there is an eo{r],L,a) so that if < £ < Eq and v is 
any stationary distribution for the nonlinear voter model satisfying 

v[^ = Q or i = l) = 0, then 



sup 

X 



Ke(x) = 1) - ^ 



< rj. 



Remark. The proof is given in Section 1.7.3 below. Case 4A is of particular interest 
as there is coexistence even though /'(O) < 0. Here the low density limit theorem 
in [7] shows convergence to super-Brownian motion with drift /'(O) < (see the 
discussion in Section 1.8 below). From this one might incorrectly guess (after an 
exchange of limits) that there is a.s. extinction of O's for the nonlinear voter model, 
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while our proof of coexistence will show that there is positive probability of survival 
of O's even starting with a single 0. 

We next show that more can be said in Case 2 if we break the symmetry. This 
also demonstrates how one can handle higher degree reaction functions in the pde 
and still apply the general results in the next Section. Consider the nonlinear voter 
model ^ as before but now for A > replace gi with 

4 

5i,A(6,---,e4)) = (1 + A)a(j^e*), (1-58) 

1 

while go is unchanged. To avoid trivialities we assume Yli <^ij) > 0- A short 
calculation now shows that if / is as in (1.52), then our reaction function in (1.16) 
becomes 

4 

f^x){u) = /(n) + A^a(j>,(n) = f{u) + Xfo{u) > f{u) on (0,1). (1.59) 
i=i 

Decomposing qj{u) as before we get 

f(\){u) = fi,x{u) + f2,x{u), 

where 

fi,x{u) = ih + 4Aa(l))n(l - + (62 + 6Xa{2))u\l - uf (1.60) 
- (62 - 4Aa(3))u3(l - uf - (61 - \a{A))u\l - u) 
= fi{u) + Xh{u) > h{u) on (0,1), 
5 

f2,x{u) =y2e{j,L,X)u^ , and sup |e(j, L, A)| < C2(A + l)7/o(i). (1.61) 
We also have 

flxiO) = 61 + 4Aa(l), /;,;,(1) = 61 - Aa(4). (1.62) 
Theorem 1.14. Suppose 61 > and < A < 6i/a(4). 

(a) Coexistence holds for large L and small enough e (depending on L, A and a), 
(h) Assume 36i + 62 < and let 1 — a' denote the largest root of fi{u) = in (0, 1). 
If r] > 0, L > £1(77,0)^^, < A < ei{r],a), < e < eo{7], L, X,a) and u is any 
stationary distribution satisfying = or ^ = 1) = 0, then 



sup 



v{i{x) = 1) - (1 - a') 



< rj. 



See Section 1.7.3 for the proof. For a concrete example, consider a(l) = a(2) = 1, 
a(3) = a(4) = 3, which is a version of the majority vote plus random flipping. Then 
5i = 4a(l) - a(4) = 1, 62 = 6a(2) - 4a(3) = -6, and 36i + 62 = -3, and so the 
hypotheses of (b) hold for small A > 0. As A ^ 0, the density of any invariant 
measure approaches 1 — a the largest root of fi{u) = 0, while symmetric reasoning 
shows that if A f 0, the densities will approach a. Of course the closer A gets to 0, 
the smaller we must make £ to obtain the conclusion of Theorem 1.14, so we are not 
able to prove anything about the case A = 0. 
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1.6 General Coexistence and Extinction Results 

Our results for the three examples will be derived from general results with hy- 
potheses concerning properties of the limiting PDE. In this subsection, we state 
those results and give an overview of the contents of the rest of the paper. We work 
with a voter model perturbation ^f{x),x G 811^,1 > throughout. 

In Section 2 we first introduce a family of Poisson processes ( "graphical repre- 
sentation") which we use to define our process (x) on x G eZ'^. Using this and 
working backwards in time we define a "dual process" X which is a branching coa- 
lescing random walk with particles jumping at rate according to p^{x) = p{x/e) 
and with a particle at x giving birth to particles at x -|- eY^, . . .x + eY^° when 
a reaction occurs at x. The ideas in the definition of the dual are a combination 
of those of Durrett and Neuhauser [20] for systems with fast stirring and those of 
Durrett and Zahle [23] for biased voter models that are small perturbations of the 
voter model. 

Duality allows us to compute the value at z at time T by running the dual process 
backwards from time T to time starting with one particle at z at time T. Most 
of the work in Section 2 is to use coupling to show that for small e, X^ is close to 
a branching random walk X^. Once this is done, it is straightforward to show (in 
Section 3) that as e — t- the dual converges to a branching Brownian motion X'^, 
and then derive Theorem 1.2 which includes convergence of (x) = 1) to the 
solution u(t,x) of a PDE. 

Our general coexistence result will be based on the following assumption about 
solutions to the PDE. The coexistence results for the models discussed in the pre- 
vious section are obtained by verifying this assumption in the particular cases. 

Assumption 1. Suppose that there are constants < vq < u^, < u* < vi < 1, and 
w, Li > 0, so that 

(i) ifu{0,x) > vq when \x\ < Lq, then liminft_5.oo inf|a;|<t„tti(t,x) > u^. 
(a) ifu{0,x) < vi when \x\ < Li, then limsup^_i.oo sup|2.|<^( u(t, x) ^ 

We also will need a rate of convergence in (1.9), namely for some tq > 0, 

1 

Y.\\9t-9^\\oo<Cl,Q^e'\ (1.63) 

1=0 

Assumption 1 shows that the limiting PDE in Theorem 1.3 will have solutions 
which stay away from and 1 for large t. A "block construction" as in [14] will be 
employed in Section 6 to convert this information about the PDE into information 
about the particle systems. In effect, this allows us to interchange limits as e — )• 
and t — 7- oo and conclude the existence of a nontrivial stationary distribution, and 
also show that any stationary distribution will have particle density restricted by 
the asymptotic behavior of the PDE solutions at t = oo. 

Both Theorem 1.15 and 1.16 below will be formulated for the voter model pertur- 
bations on eZ'', but the conclusions then follow immediately for the original unsealed 
particle systems in (1.4). 
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Theorem 1.15. Suppose Assumption 1 and (1.63). If e > is small enough, then 
coexistence holds and the nontrivial stationary distribution v may he taken to he 
translation invariant. 

If 7] > and e > is small enough, depending on rj, then any stationary distri- 
bution V such that 

'^{Y,ax) = 0or^{l-ax)) = 0)=0 (1.64) 

X X 

satisfies z^(^(x) = 1) G (u* — rj,u* + rj) for all x. 

Note that in the second assertion we do not require that be translation invari- 
ant. 

Results that assert O's will take over will require a stronger pde input: 

Assumption 2. There are constants < tii < 1, C2,C2,w > 0, Lq > 3 so that for 
all L > Lq, if u{0, x) < ui for \x\ < L then for all t > 

u{t,x) < C2e"'^2i all \x\ < L + 2wt. 

Finally we need to assume that the constant configuration of all O's is a trap for 
our voter perturbation, that is, 

5^(0, ... ,0) = 0, or equivalently /i^(0,0) = 0, for < e < eo , (1-65) 

where is the zero configuration in {0, 1}^'' . This clearly implies /(O) = and is 
equivalent to it if gf does not depend on e, as is the case in some examples. Recall 
the definition of "i"s take over" from (1.14) and that q is the law of (y\ . . . , Y^"). 

Theorem 1.16. Suppose Assumption 2, (1.63), (1.65), p(-) and q{-) have finite 
support, and /'(O) < 0. Then for e small the O's take over. 

We believe the theorem holds without the finite range assumptions on p and q. 
By Proposition 1.1 in the above finite support setting, it suffices to assume (1.10) in 
place of (1.7), and also assume (1.63) holds for the ,5i appearing in (1.10) instead 
of the gl,gi. 

In order to show that O's take over, say, we will need several additional argu- 
ments. The pde results will only ensure we can get the particle density down to a 
low level (see Section 4) but clearly we cannot expect to do better than the error 
terms in this approximation. To then drive the population to extinction on a large 
region with high probability we will need to refine some coalescing random walk 
calculations from the low density setting in Cox and Perkins [8] -see Section 7 and 
especially Lemma 7.6. This is then used as input for another percolation argument 
of Durrett [I-!] to guarantee that there are no particles in a linearly growing region. 
Since there was an error in the original proof of the latter we give all the details 
here in Sections 5 and 7.3. 

Quantifying the outline above, the first step in the proof of Theorem 1.16, taken 
in Section 4 is to use techniques of Durrett and Neuhauser [20] to show that if has 
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density at most ui on [—L, L]'^ then at time Ti = ci log(l/e) the density is at most 
on [—L — wTi,L + wTi]'^ (see Lemma 4.2). Here /3 is a small parameter. The second 
step, taken in Section 7, is to show that if one waits an additional T2 = C2 log(l/e) 
units of time then there will be no particles in [— L — wTi + AT2, L + wTi — A2T2Y' 
at time Ti + T2. The first step here (Lemma 7.6) is to show that if we start a 
finite (rescaled) block of ones of density at most then with probability at least 
1 — e^/^ it will be extinct by time Clog(l/e). Here it is convenient that arguments 
for the low density regime of [7] (density e'^"^) continue to work all the way up to 
and also that the PDE arguments can be used to reduce the density down to 
. In short, although the precise limit theorems, Theorem 1.2 and Corollary 1.8 
in [7] apply in disjoint regimes (particle densities of 1 and e*^"^, respectively) the 
methods underlying these results apply in overlapping regimes which together allow 
us to control the underlying particle systems completely. Of course getting I's to be 
extinct in a large block does not give us what we want. The block construction in 
[20] is suitably modified to establish complete extinction of I's on a linearly growing 
set. A comparison result of [-34], suitably modified to accommodate our percolation 
process, is used to simplify this construction. 

1.7 Application to the examples 
1.7.1 Lotka-Volterra systems 

Proof of Theorem 1.9. (i) Let < < 1 and consider first 

ao = ckq = 1 — e^, ai = a\ = 1 — mo(l — ??)e^, 

so that in the notation of Section 1.3 we have set = —1) ^1 = —"^0(1 — v)- The 
rescaled Lotka-Volterra process is a voter model perturbation and from (1.36) we 
have 

f{u) = -n(l - u)[rip2 + tip3(l + "io(l - fi))] < on (0, 1). 

Proposition 1.5 verifies Assumption 2 in Theorem 1.16 and /'(O) < is obvious. 
(1.63) is trivial {gl = gi) and (1.65) is immediate from (1.34). The finite range 
assumption on g = p x p is immediate from that on p. Theorem 1.16 implies O's take 
over for e small. Therefore when < 1 — ao < ''o(??) a-nd ai = 1 + mo(l — ??)(ao — 1), 
then O's take over for LV{a). The monotonicity in (1.31) shows this is also the case 
for ai > 1 + ?7T.o(l — 7?)(ao — 1) and ao as above. 
Next consider 

aQ = = 1 + , ai = af = 1 + (1 + ?/)e^ that is, ^0 = 1; = 1 + r?. 

In this case we have 

f(u) = u{l - u)[p2 - (1 + ■ri){p2 +pz) + np3(2 + r/)], 

and so, assuming without loss of generality (by (1.31)) 1 + rj < rriQ^, from (1.37) / 
has a zero, and an unstable fixed point for the ODE, at 

* (1 +r?)(p2 +P3) -P2 ^ f'^ A 

U = -. r G -, 1 . 

P3(2 + ??) V2' ; 
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It follows that Jq f{u)du < and Proposition 1.7 establishes Assumption 2. As 
above, Theorem 1.16 and (1.31) show that O's take over if < ao ~ 1 is sufficiently 
small and ai > 1 + (1 + r/)(ao — !)• 

(ii) Interchange the roles of and 1 in (i). □ 

Proof of Theorem 1.8. We slightly extend the setting in Section 1.3 and for 77 E 
(0, 1 — mo) consider 

a^ = l-e2,a^ = l + £2^1, where - 0^ G [-^, ^— ^1 , lim 0f = 0i. (1.66) 



1 — ?7 mo 



£4,0 



Then the rescaled Lotka-Volterra model in (1.32) remains a voter model pertur- 
bation but now may now depend on e. From (1.36) we have 

f{u) = n(l - u)[-p2 - 9i{p2 +P3) + 'WP3(-1 + ^'i)], 
which has a zero, and attracting fixed point for the ODE, at 

n(-^i) = ^^^^^-^G(0,l). (1.67) 

Proposition 1.4 and its mirror image, with and 1 reversed, establish Assumption 1 
with u* = = u*{—6i) (see (1.27)). Theorem 1.15 therefore shows that for < 
e < eo{r]) 

coexistence holds, and if zy is a stationary distribution satisfying (1.68) 
i/(^ = or ^ = 1) = 0, then sup |i/(^(x) = 1) - u*{-ei)\ < r]. 

X 

Suppose first that (ii) of Theorem 1.8 fails. Then there is a sequence e„ | 0, 
{al",al") and k > so that (1.66) holds with £ — £n, and there is a stationary 
measure Vn for S^^" satisfying f„(^ = or ^ = 1) = and such that 

sup\i^n{({x) = l)-u*{-9l")\ > K. 

X 

Since u*{—9l") — )■ u*{—6i), if we choose r] < k this contradicts (1.68) for large n, 
and so proves (ii). The proof of (i) is similar using the first part of (1.68). That is, 
if (i) fails, there is a sequence En i so that coexistence fails for a^" as in (1.66), 
contradicting the first part of (1.68). □ 



1.7.2 Evolution of cooperation 

Here is the result on coalescing probabilities which will help us simplify the formula 
(1.51) for the reaction function /. The notation is as in Section 1.4. 

Lemma 1.17. (a) p{ei\e2) = p{0\ei) . 
(b)piei\e2 + e3) = (l + i)p(0|ei). 
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Proof. Let Bf denote a rate 2 random walk with kernel p starting at x, and note 
that for X 7^ 0, — IB^ has the same law as B^ until it hits 0. Note also that for 
X / 0, P{Bf / V t > 0) = EyPiy)P{Bt^^ / V t > 0). 
For (a), 

p(0|ei) = / for ah t > 0) 

XI 

= ^ ^p(2;i)p(2;2)P(Bf 1+^2 _A for all t > 0) (use xi / 0) 

= '^'^p{xi)p{x2)P{Bt^~'^'^ / for all t > 0) (by symmetry) 
= p{ei\e2). 

For (b), let Tj{x) be the time of the jth jump of B^. Then using symmetry, 
p{ei\e2 + eg) = ^ p{xi)p{x2)p{x3)P{B^^+''^+''' / for ah t > 0) 

= ^p{xi)P{Bt' / for all t > T2{xi)). 

XI 

Now using the above and first equality in the proof of (a), 

p(ei|e2 + 63) -p(0|ei) = J^p(xi)P(^^i = 0, B^' / for ah t > Ti{xi)) 

Xl 

= k-^^p{xi)P{B^ / for all t > Ti) 

Xl 

= fc-yoiei). 

The result follows. □ 
Proof of Lemma 1.10. We first rewrite (1-51) as 

^ =(/3 - 7)(e(ei)e(e2)). + (7 - '^)(e(ei)|(e2)|(e2 + eg) + i{ei)i{e2)i{e2 + eg)). 

+ ((a - /3) - (7 - mi{ei)i{e2)i{e2 + eg)). 
=1 + 11 + 111. (1.69) 

Some elementary algebra shows that 

II={j- 5)(e(ei) - e(ei)e(e2) - ^(61)^(62 + eg) + ^{62)^(62 + eg)).. (1.70) 

Note that (1.26) and Lemma 1.17(a) imply 

(C(ei)^(e2))« =^^^(61162) + u{l -p(ei|e2)) 
=n2p(0|ei) + u(l-p(0|ei)) 
=(^(0)e(ei))« = (^(e2)^(e2 + eg))„, 
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the last by translation invariance. Using this in (1.70) and again applying (1.26), 
we get 



/ + // = (/3 - 7)n(l - n)p(ei|e2) + (7 - <5)[n - (e(ei)C(e2 + eg))^] (1.71) 
= (/3 - 7)n(l - u)p{ei\e2) 

+ (7 - - v?p{ei\e2 + 63) - n(l -p(ei|e2 + 63))] 
= n(l - n)[(/3 - 7)p(ei|e2) + (7 - 5)p(ei|e2 + 63)] 
= n(l - n)p(0|ei)[(/3 - 7) + (1 + A:-^)(7 - <5)], 

where Lemma 1.17 is used in the last equality. A straightforward application of 
(1.26) allows us to find the coefficients of the cubic III in (1.69) and we obtain the 
required expression ior f{u)/k. □ 

Proof of Theorem 1.11. This is now an easy application of Theorem 1.16. Assume 
7 — 5 < k{6 — /3). Then Lemma 1.10 shows that f{u) = c\u{\ — u) for ci < 0. 
Proposition 1.5 shows that Assumption 2 of Theorem 1.16 is valid for any ui G (0, 1). 
The condition (1.63) holds with ro = 2 by (1.47) (recall \\gl — gi\\oo = H^f — 5j||oo by 
Proposition 1.1). The condition (1.65) is clear from the expression for /if in (1.44). 
Since /'(O) < is clear from the above, and w = , Theorem 1.16 completes the 
proof in this case. The case where the inequality is reversed follows by a symmetrical 
argument, or, if you prefer, just reverse the roles of and 1. □ 



1.7.3 Nonlinear voter models 

Proof of 1.13. Consider Case 4A first. As pointed out in this Case in Section 1.5, for 
L sufficiently large we may employ the mirror image of Proposition 1.7 on [0, 1/2] 
with p = a, the unique root of / in (0,1/2), and Proposition 1.7 on [1/2,1] with 
p = 1 — a, along with the comparison principle (Proposition 2.1 in [2]), to see that 
Assumption 1 holds for e < eo(^) with = ^ — 77, u* = ^ + vq = S, and vi = 1 — 6. 
(1.63) is trivial because gf = gi. Theorem 1.15 now implies (a) and (b) in this case. 
The proofs in Cases 1 and 2 are similar using Proposition 1.4 (note all the zeros are 
simple in these cases) to verify Assumption 1 (see the discussion in these cases in 
Section 1.5). □ 

Proof of 1.14. (a) The conditions on bi and A imply that /( -^{0) > and f[ ^(1) > 
0. Coexistence for large L and small e is now established as in Case 2 (or 1) of 
Theorem 1.13-see the discussion in Section 1.5. 

(b) By taking A and small, depending on {rj,a), we see from (1.60), (1.61), 
(1.62), and our conditions on the bi that /(a) (^i) = will have 3 simple roots in 
(0, 1), pi(A) < P2(A) < P3(A), within ry/4 of the respective roots 

a < 1/2 < 1 - a' 

of /i(ti) = 0. As (1.63) is again obvious, we now verify Assumption 1 of Theorem 1.15 
with u* = 1 — a' + ^ , ti* = l — a' — ^,vo G {p2,u*), and vi £ {u* ,1) (r/ is small so these 



28 



intervals are non-empty). The result would then follow by applying Theorem 1.15. 
The upper bound (ii) in Assumption 1 is an easy application of Proposition 1.5, 
with the interval {p^, 1) in place of (0, 1), and the comparison principle. 

For the lower bound (i) in Assumption 1, we use a result of Weinberger [42]. To 
state the result we need some definitions. His habitat T-L will be M'^ in our setting 
and his space B is the set of continuous functions from Ti to [0, vr+j. In our case 
7r+ = 1. His result is for a discrete iteration Un+i = Q{un), where in our case Q{u) 
is solution to the PDE at time 1 when the initial data is u. His assumption (3.1) 
has five parts: 

(i) ifueB then Q{u) G B. 

(ii) If Ty is translation by y then Q[Tyu) = TyQ[u). 

(iii) Given a number a, let Q{a) be the constant value of Q{ua) for = a. There 
are < ttq < vri < 7r+ so that if a € (vro,vri) then Q{a) > a. Qiiro) = ttq and 

Q{tTi) = VTi. 

(iv) u < V implies Q{u) < Q{v). 

(v) \i Un ^ B and n.„ — )• u uniformly on bounded sets then Q{un){x) — >■ Q{u){x). 

Clearly (i) and (ii) hold in our application. For (iii) we let ttq = P2(A) and tti = 
P3(A). (iv) a consequence of PDE comparison principles, see, e.g., Proposition 2.1 
in Aronson and Weinberger (1978). (v) follows from the representation of solutions 
of the PDE in terms of the dual branching Brownian motion (see Lemma 3.3). 
The next ingredient for the result is 

S = {xeM.'^ ■.x-i< c*(0 for all ^ G S'^'^], 

where S'^~^ is the unit sphere in M*^. c*(^) is the wave speed in direction ^ defined 
in Section 5 of [42]. Due to the invariance of the PDE under rotation, all our speeds 
are the same, c* = p, and 5 is a closed ball of radius p or the empty set. Here is 
Theorem 6.2 of [42]. 

Theorem 1.18. Suppose (i)-(v) and that the interior of S is nonempty. Let S" he 
any closed and hounded suhset of the interior of S. For any 7 > ttq, there is an 
so that if uo{x) > ^ on a hall of radius r^ and if Un+i = Q{un) then 

liminf min Un{x) > vri. (1-72) 

n-5>oo x&nS" 

To be able to use this result, we have to show that p > 0. Note that here 
we require lower bounds on the wave speed of solutions to the reaction diffusion 
equation in one spatial dimension. This is because traveling wave solutions in the 
direction of the form w{x ■ ^ — pt) correspond to traveling waves w in one spatial 
dimension. Recall that in the decomposition (1.59) f{u) is odd about u = 1/2, and 
for large L has /'(1/2) > by 3bi + 62 < 0. The latter shows / has 3 simple zeros in 
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(0,1) at a < 1/2 < 1 — a. The strict inequality in (1.59) on (0, 1) now easily implies 
(compare the negative and positive humps separately) 

/ /(A)(n)dn> / f{u)du = 0. (1.73) 

J pi J a 

So by the discussion in part II (ii) of Section 1.2 there is a one-dimensional decreasing 
traveling wave solution to (1.19) (with / = /(a)) over (pijPs) with positive wave 
speed r2(A). 

To consider traveling waves over (0,pi(A)), we note that Kolmogorov, Petrovsky, 
and Piscounov [31] have shown that if we consider 

— = — An+ f(u) 

in one dimension where / satisfies 

/(O) = /(I) = 0, f{u) > for < w < 1, f'{u) < /'(O) for < u < 1 (1.74) 

then there is a traveling wave solution with speed y^2(T^/'(0) and this is the minimal 
wave speed. For this fact one can consult Bramson [-'5] or Aronson and Weinberger 
[1]. However, the intuition behind the answer is simple: the answer is the same as 
for the linear equation 

dt 2 ^ ' 

which gives the mean of branching Brownian motion. For more on this connection, 
see McKean [.'>5]. 

Now let gi < §2 be functions on [0, 1] such that 

< 5'2 < /(A) on (0,pi), gi= g2 = f(^x) on [pi, 1], 

52(0) G (o, ^-^^ , g',{u) < g'M on [0,pi], (1.75) 
and for some < j^o < Pi j 

51(0) = 0, 51 < on (0,po),5i > on (po,Pi), 9i{t) dt > 0, g[{0) < 0. (1.76) 

Jo 

The existence of such functions is elementary. By the KPP result above, the minimal 
wave speed over (0,pi(A)) for the 52 equation is 

C2 = ^2^25^(0) < r2(A). (1.77) 

By Theorem 2.4 and Corollary 2.3 of [25] (or the discussion in part Il(ii) of Sec- 
tion 1.2) there is a unique traveling wave solution u{t, x) = w^x — cit) {w decreasing) 
to the gi equation with unique wave speed ci > (since the integral in (1.76) is 
positive) and range (0,pi). Note here and elsewhere that the traveling waves w in 
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[25] are increasing and so our wave speeds have the opposite sign. By a comparison 
theorem for wave speeds (Proposition 5.5 of [42]) we may conclude that 

C2 > ci. (1.78) 

The hypothesis of the above comparison result is easily verified using gi < g2 and 
the standard comparison principle (e.g. Proposition 2.1 of [2]). It follows from 
(1.77) and (1.78) that ci < r2(A) which are the wave speeds of the gi equation over 
(0,pi) and (pijPs), respectively. We can therefore apply Theorem 2.7 of [25] to 
conclude the existence of a traveling wave over (0,p3(A)) for the gi equation with 
speed ri(A) G (ci,r2(A)). The wave and its speed are both unique by Corollary 3.3 
of [25]. Since f(^x) — 9i on [0,^3], another application of Proposition 5.5 of [42] shows 
that p > ri(A) and in particular p > 0. 

Using (1.72), we have proved that for < = p, 

liminf inf u{n,x) > P3(A) > 1 — a' — 7, 

n-!>oo \x\<2wn 4 

providing that u{0,x) > vq for < r^^. The same reasoning gives the same 
conclusion with rer in place of n for any r > 0. Taking r small enough, a simple 
interpolation argument (use the weak form of the reaction diffusion equation and 
smoothing properties of the Brownian semigroup) now gives Assumption l(i) with 
u* = 1 — a' — ^ where the 2w in the above helps a bit in this last interpolation step. 
□ 

1.8 Comparison with low density superprocess limit theorem 

To make a comparison between our hydrodynamic limit theorem (Theorem 1.3) and 
the superprocess limit theorem of Cox and Perkins [7] we will write our perturbation 
terms in a different form, which will also be useful in Section 7. Define 

^siv) = n for r/ = (?7i, . . . , rjNo) G {0, 1}"^° ,SeVNo= {subsets of {1, ... , A^o}}, 
ieS 

and 

A G Vno = {subsets of Z'^ of cardinality at most Nq}. 
By adding an independent first coordinate to Y we may assume has law p. If 

fff(6, ■■■,Cno) = -ei^m = i) +91(^1, ■ ■ ■,Cno), (1-79) 
and (ji is as above without the superscript e, then 

lim||gf -5,||oo =0, (1.80) 



31 



and we may rewrite (1.7) as 

hl{x,i) = EY{gt{x + Y\...,x + Y''^)), i = 0,l, (1.81) 

and similarly without the e's. It is easy to check that {Hs'(-) : S G Vno) is a basis for 
the vector space of functions from {0, 1}^'' to M and so there are reals I3^{S), 5^{S), 
S G Vnq-, such that 

~9{{ri)= J2 ^e{S)Es{v), Uri)= E US)'^s{ri). (1.82) 

and similarly without the e's. If 5 G Vno, let = {Y^ : i e S], where Y £ Z'^^o 
has law q as usual. Let Ey denote expectation with respect to Y. It is easy to use 
(1.81) to check that 

hl{x,0= ^ Pe{S)EY{xiY',x,0)= J2 f3s{A)x{A,x,0 (1.83) 
h'o{x,0= E SeiS)EY{x{Y^x,0)= Yl SeiA)x{Ax,0, (1.84) 



where for A £ V 



No, 

I3,{A)= E ^e{S)P{Y'^ = A), 6e{A)= E S,{S)P{Y'' = A). (1.85) 

Analogous equations to (1.83), (1.84) and (1.85) hold without the e's. 
Now use (1.83) and (1.84) without the e's, and (1.26) to see that 

f{u) ^ ((1 - mhi{o,o - mho{o,o)u 

= E [(3{A)[Yu^{l-u)P{\ii\=j,r{A,{0})=^)\ (1-86) 
a&Pnq i=i 

|AU{0}| 

+ ^(0)(l-n)-5(A)[ E u^P{\ii''^'M=j)\\, (1.87) 

which is a polynomial of degree at most A''o + 1 as claimed in Section 1.1. If /3(0) = 0, 
then /(O) = and 

/'(0)= E PiA)PiT{A)<oo,T{A,{0})=oo)-5{A)PiT{Au{0})<oo). (1.88) 

From (1.82) one easily derives 

4(5) = E (-l)l^l-l^lfff (ly), Us) = E i-lf^-^'^Wv), (1.89) 
vcs vcs 
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and similarly without the e's. Therefore 

mS) - /3(5)| + \US) - SiS)\ < 2^»(||5f - 51II00 + \\g'o - Moo), (1.90) 
and so 

J2 |/3.(A)-/3(A)| + |<5,(^)-<5(A)| <22^"(||5f-9i||oo + ||5^-ffo||oo). (1.91) 
Our spin-flips are now recast as 

c,(ex,ee) = e-V(x,e)+ Yl xiA,x,omA)ii-a^))+6e{A)ax)], 

A£Vno 

which is precisely (1-17) of [7] with e = N~^/'^ . If we assume 

(0) = (and hence g\{<d) = 4(0) = 0) for smah e, (1.92) 

and the voter kernel p has finite support, then using the fact that the right-hand side 
of (1.91) approaches as e — )• (by (1.80)), it is easy to check that all the hypotheses 
of Corollary 1.8 of [7] hold. Alternatively, in place of the finite support assumption on 
p one can assume the weaker hypothesis (P4) of Corollary 1.5 of [S], and then apply 
that result. These results state that for e as above if Xf = Z^xseZ'^ Q{x)^x and 
— Xq weakly in the space Mpi^"^) of finite measures on W^, then X^ converges 
weakly in the Skorokhod space of Mi?(]R'^)-valued paths to a super-Brownian motion 
with drift 6 = /'(O) (as in (1.88)). In this result we are starting 0(e~^) particles 
on a grid of £~'^ {d > 3) sites per unit volume, so it is a low density limit theorem 
producing a random limit, whereas Theorem 1.3 is a high density limit theorem 
producing a pde limit. The latter result gives a natural explanation for the drift 
9 in the super-Brownian limit which was defined by the right-hand side of (1.88) 
in [7]. Namely, under (1.92), in the low density limit we would expect a drift of 
lim„^o f{iJ-)/u = /'(O), which of course happens to equal the summation in (1.88). 
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2 Construction, Duality and Coupling 



In this section, we first introduce a family of Poisson processes which we use to 
define on eZ,'^, a "dual process" X and a "computation process" C- The duality 
equation (2.17) below gives a representation of ^ti^) in terms of (X, ("). Next we 
show that for small e, {X, C) is close to the simpler {X, where X is a branching 
random walk system with associated computation process C- Finally we show by a 
strong invariance principle that for small e, {X, (^) is close to a branching Brownian 
motion and its associated computation process. 

However, our first task will be to prove Proposition 1.1 and reduce to the case 
where ei = oo in (1.7). 

2.1 Preliminaries 

Proof of Proposition 1.1. Let p = mva{p[yi) : p{yi) > 0}, choose Eq > so that 
M = supq^j<;(,jj II^oIIoo V ll^f lloo < CO and then choose ei > so that 

e^^p > M. (2.1) 

For < e < £0 define on {0, 1}^° by 

No 

gli^i, ■ ■ ■ ,Cno) = e^^^HCj = i)piyj) + gfiCi, ■ ■ ■ ,Cno), i = o,i, (2.2) 
1 

and define gi by the same equation without the e's. Clearly \\g? — giWoo = Wgf — 
giWoc — )• as e — )• 0. We may assume yi = 0. By replacing gf with ^fl(^i = 1 — i) 
and redefining hf analogously (this will not affect (1-5)), we may assume 

gK^u...,CNo) = Oif^i = i. (2.3) 

We now show that gf > 0. Assume first 

No 

J^^iP{yi) = 0- (2.4) 
1 

Choose ^ G {0, 1}^" so that C(2/i) = Ci- If C(0) = 0, then by (1.5), (1.10), and (2.2), 

< c,(0,ee) = 5i(e(yi), • • • ,^(m,)) = 5f (6, • • • ,^JVo)- 
If ^(0) = 1, then ^1 = ^(0) = 1 and by (2.3), (^i, . . . , ^atq) = 0. Assume next that 

No 

1 

Then the above sum is at least p and so 

5i(6,---,e7Vo) >er'p-ii5iiioo >er'p-^>o, 
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the last by (2.1). This proves gf > and a similar argument shows gQ > 0. Finally 

(1.7) with y* = Ui is immediate from (1.10) and the definition of gf. □ 

We claim we may assume without loss of generality that ei = oo in (1.7), that is, 
the first term in the right-hand side of (1.7) is absent. To see why, let = e~^— e^^ 
for e < El, and use (1.7) in (1.5) to rewrite the spin- flip rates of as 

c%ex,Ce) = e-^c^x,0 + (1 - ^(x))/if (x,^ + e(x)/i§(x,e), 

where 

hiix, = Evigm^ + y^)^ ...,i{x+ y^°))). (2.5) 

So by working with in place of /if throughout, we may use (2.5) in place of (1.7) 
and effectively set ei = oo. Note first that this does not affect the definition of the 
reaction term f{u) in the PDE (1.19) since the terms involving e~'^fi{x,^) cancel 
in (1.16). The only cost is that is replace with The ratio of these terms 
approaches 1 and so not surprisingly this only affects some of the proofs in a trivial 
manner. Rather than carry this with us throughout, we prefer to use e and so 

henceforth set ei = oo in (17). (2.6) 

2.2 Construction of 

Define c* = c*{g) by 

c* = sup ll^llloo + llfi-olloo + 1- (2.7) 
0<e<eo/2 

To construct the process, we use a graphical representation. For x G eZ"^, introduce 
independent Poisson processes {T^,n > 1} and {Tn^,n > 1} with rates and c*, 
respectively. Recall Pe(?/) = p{y/e) for y S eZ*^ and let qe{y) = q{y/£) for y E eZ'^^o. 
For X € eZ*^ and n > 1, define independent random variables Zx^n with distribution 
Pe, Yx,n = {Yx,ni ■ ■ ■ '^ivn) with distribution q^, and Ux,n uniform on (0, 1). These 
random variables are independent of the Poisson processes and all are independent 
of an initial condition £ {0, 1}'"^ . 

At times t = T^,n > 1 (called voter times), we set S,t{x) = Ct-{x + Zx^n)- To 
facilitate the definition of the dual, we draw an arrow from (x,T^) — )• {x + Zx^mT'n)- 
At times t = Tn^, n > 1 (called reaction times), if ^t-(x) = i we set .^t(x) = 1 — i if 

Ux,n < gf-iiCt-ix + YlJ, . . . ,6-(x + yjy»))/c*, and otherwise Ct{x) = Ct-{x). 

At these times, we draw arrows from (x,T^'^) {x + Yx^^iTn^) for 1 < i < A^'o- 
We write a * next to (x,T^'^) and call these *-arrows. It is not hard to use ideas of 
Harris [29] to show that under the exponential tail conditions on p and q^ (1-1) and 

(1.8) , this recipe defines a pathwise unique process. This reference assumes finite 
range interactions but the proof applies in our infinite range setting as there are 
still finitely many sites that need to be checked at each reaction time. To verify 
this construction and to develop a useful dual process we now show how to compute 
the state of x at time t by working backwards in time. It is easy to verify that 
^ is the unique in law {0, 1}^ -valued Feller process with rates given by (1.3) and 
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(1.5), or more precisely has generator as in (1.12). For example one can recast the 
graphical representation in terms of SDE's driven by Poisson point processes and 
use stochastic calculus as in Proposition 2.1(c) of [8] (it is easy to verify condition 
(2.3) of that reference in our current setting). 

We use B'^'^ to denote a continuous time random walk with jump rate and 
jump distribution starting at x G eZ'^ and drop dependence on a; if x = 0. We 
also assume 

{S^'^ : X G eU^} are independent random walks distributed as above. (2-8) 

It will be convenient to extend the Poisson times to the negative time line indexed 
by non-positive integers, and hence have {T^,n G Z}, {Tn^,n £ Z} with the as- 
sociated {Zx.n,n G Z} and {{Yx^n,Ux,n),'n S 2}j respectively. At times it is useful 
to work with the associated independent Poisson point processes of reaction events 
A^{dt, dy, du) {x G eZ'^) on M x eZ'^^o x [0, 1]) with points {(r^'"", Y^^n, Ux,n)} and in- 
tensity c*dt X X du, and also the independent Poisson point processes of walk steps 
A^{dt,dz) {x G eZ"^) on M x eZ"^ with points {(T^, Zx^n)} and intensity e~'^dt x p^. 

2.3 The Dual X 

Fix r > and a vector of M -|- 1 distinct sites z = [z^, . . . , z*^), each Zi G eZ"^. Our 
dual process X = X^'^ starts from these sites at time T and works backwards in 
time to determine the values ^T(zj). X will be a coalescing branching random walk 
with Xq = (zq, . . . , zm, oo, . . . ) taking values in 

V = {(X°,X\...) G D([0,r],M'^U{oo})^+ : 

3Ko G Z+ s.t. Xj; = oo Vt G [0,r] and k > Kq}. 

Here oo is added to M"^ as a discrete point, D{[0, T], R'^Ujoo}) is given the Skorokhod 
Ji topology, and V is given the product topology. 

For X = {X^,X^,...) G V, let K{t) = max{i : XI / oo}, define i r^t i' 
iS XI = XI 7^ oo, and choose the minimal index j in each equivalence class in 
{0, . . . K{t)} to form the set J{t). We also introduce 

I{t) = {XI : i G J{t)} = {XI : Xj / oo}. 

Durrett and Neuhauser [20] call I{t) the influence set because it gives the locations 
of the sites we need to know at time T — t, to compute the values at z^, . . . , z^ at 
time T. 

To help digest the definitions, the reader should consult Figure 6, which shows 
a realization of the dual starting from a single site when A'^o = 3. If there were no 
reaction times Tn^ then the coordinates Xl,j G J{t) follow the system of coalescing 
random walks dual to the voter part of the dynamics. Coalescing refers to the fact 
that if Xi = Xi for some s and then Xl = Xf for all t G [•s,r]. Jumps 
occur when a particle in the dual encounters the tail of an arrow in the graphical 
representation. That is, if j G J{s—) and x = Xl_ has T — s = then Xi = x + Z^. 



36 



dual Xt 




Figure 6: An example of the dual with Nq = 3. 
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It coalesces with XI = x + Z^ \i such an i exists, meaning that i V j is removed from 
J{s—) to form J{s). If is a rate random walk on eZ'^ with step distribution 
Pfr then the coalescing random walks in the dual X follow coalescing copies of -B^. 

To complete the definition, we have to explain what happens at the reaction 
times. Put Rq = 0, and for m > 1 let Rm be the first time t > Rm~i that a particle 
in the dual encounters the tail of a *-arrow. If 

j G J{Rm — ) and x = Xj^^_ has T — Rm = Tn'* for some n, (2.9) 

we let = j denote the parent site index. In the example in Figure 6 ni = 0, 
fi2 = 1, ^3 = 3, and fi^ = 2. 

We create Nq new walks by setting = ^, 1 < i < Nq, 

KiR^) = K{R„,.i) + No, and 

The values of the other coordinates X^ , / G J{Rm—), j' ^ fim remain unchanged. 
Each "new" particle immediately coalesces with any particle already at the site 
where it is born, and we make the resulting changes to J(i?m— ) to construct 
J{Rm) ^ J{Rm — )- To compute £,t{z^)-, we will also need the random variables 

Um = Ux,n where x, m, and n are as in (2.9). (2-11) 

This computation is described in the next subsection. 

K{s) changes only at reaction times and always increases by exactly Nq, so 

K{s) = M + mNo, for s G [R^,R^+i). (2.12) 

Let J-t be the right-continuous (time reversed) filtration generated by the graphical 
representation restricted to [T — t,T), but excluding the {f/x,n}- More precisely Tt 
is the right-continuous filtration generated by 

{AU[T - s,T) X A) : s <t,x £ eZ'^, A C eZ'^}, (2.13) 
{A^([r- s,r) X B x[0,l]) ■■ s <t,B C elJ^^^x G eZ^}. 

The {Rm} are then (J^()-stopping times and X is (J^t)-adapted. Since 
P{Rm+i — Rm G -IJ^R^) is stochastically bounded below by an exponential random 
variable with mean (c*(M -|- mNQ))~^ , Rm t c« a.s. (recall our graphical variables 
were extended to negative values of time) and the definition of X is complete. 
Note that 

is Tr^ - measurable and Jy„,[/„ = Ar {{T - Rm} x •). (2.14) 

As the above time reversed Poisson point processes are also Poisson point processes, 
one may easily see that 

{Ym} are iid with law and Ym is J^r^ — measurable, (2.15) 

and 

{Um} are iid uniform on [0, 1] and are independent of Too- (2.16) 
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2.4 The computation process ( 

Given an initial time to ^ [Oi^)i the coalescing branching random walk {Xs,s S 
[0,T — to]}, the sequence of parent indices {fim}, the sequence of uniforms {Um}, 
and a set of initial values in {0, 1}, CtoO) — ?to(^T-io)' ^ J{T — to), we will define 
{Cr(fc),r G [0,r],0 < /c < K{{T-r)-)} so that 

Cr(A:) = ^•(^T-r) for all r E [to, T] and < K{{T - r)-) . (2.17) 

The left hand limits here reflect the fact that we have reversed the direction of time 
from that of X. 

In general we consider a general initial state (toU) S {0, 1}, j G J(T — to)- First 
we complete this initial state by setting Cto(^) = CtoO) if ^ ~r-to J ^ "^(^ ~ ^o)- 
Suppose that for some m > 1, Rm is the largest reaction time smaller than T — to- 
The values Cr{k) do not change except at times T — Rn, so Cr = Oo for r < T — 
We decide whether or not to flip the value of ( at /x^ at time t — i?m as follows. 
Define K^ G {0, 1}^° by 

Vi = CiT-R^)-{M + {m-l)No+j), j = l,...,No. (2.18) 

Letting i = C(T-ij^)-(/^m) we set 

C(T-ii^)(/^m) = < . ^, . (2.19) 

I I otherwise. 

To update the dual now, for k < M + {m — l)NQ = K{{T —{T — Rm))—) and k ^ fim, 

if k r^R^ fim set CT-Rrn{k) = (T-R^f^m) ■ (2.20) 

Otherwise we keep C(r-ij„)(^) = C(T-r„)-(^)- 

The values Cr{k) remain constant for r G [T — Rm,T — Rm-i)- Coming to 
r = T — Rm.-i, if m — 1 > 1 we proceed as above. When we reach r = T — Rq = T 
we end by setting (t = Ct~- If ^toU) = ^i-^T-to^ 3 ^ "^(^ ~ ^o); the verification 
of (2.17) is an easy exercise from the definitions of X and 

2.5 Branching random walk approximation X 

Due to the transience of random walk in dimensions d > 3, and the fact that the 
random walk steps are occurring at a fast rate in X when e is small, any coalescing 
in X will occur soon after a branching event and close to the branching site. As in 
Durrett and Zahle [23], such births followed quickly by coalescing are not compatible 
with weak convergence of the dual. Thus we need a way to excise these events from 
X. As in [23] we define a (non-coalescing) branching random walk X and associated 
computation process Later we will couple {X, C) and {X, C) so that they are close 
when e is small. 

For m S N, denotes the set of partitions of {0, . . . , m} and for each vr S 11^, 
Jq^tt) is the subset of {0, ... , m} obtained by selecting the minimal element of each 
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cell of vr. Write i ~,r J if i and j are in the same cell of vr. Let {B^\i = 0, . . . A'^o} 
be the rate one coalescing random walk system on Z"' with step distribution p and 
initial points at = 0, y^, . . . , where {Y^, . . . , Y^°) has law q. Let uq denote 
the law on Uno of the random partition associated with the equivalence relation 
« ~ j iff B^^ (t) = B^^ (t) for some t > 0. For e > let i^^ denote the law on 
Uno of the random partition associated with the equivalence relation i ~^ j iff 
= Note that £"3/2 = £1/2^-2 this is a short amount of 

time for the sped up process. For later use when we define the branching Brownian 
motion Z we note that since e~'^^'^ — )• oo, 

i^e converges weakly to z^o as e . (2-21) 

As before we will have a fixed T > and distinct sites zq, . . . , zm in eZ"^. Our 
branching random walk X will have paths in T> and an associated set of indices 
Ji't) = {j '■ 7^ °^}- I^^t ttq £ Um be defined by the equivalence relation i ~ j 
iff B^ ^^(e^-^/^) = B^ ^j (e^'^/2). In words, ttq will be used to "mimic" the initial 
coalescence in X of the particles starting at Zj before any reaction events occur. 

For n > 1 let 7r„ G Hatq be iid with law i/g and independent of ttq. From {7r„} we 
inductively define a sequence of nonempty subsets {Jn} of Z_|_ by Jq = Jol^ro) and 
for 77, > 

Jn+l = Jn U {M + 77iVo + J : j G JoK+i) \ {0}}. (2.22) 

Set Rq = and conditional on {vr^} let {-R^+i — Rn : n > 0} he independent 
exponential random variables with means (c*| and let {fin} be an independent 

sequence of independent random variables where fJ-n, n > 1, is uniformly distributed 
over Jn-i- fi-n is the index of the particle that gives birth at time Rn- 
To define X inductively we start with 

Xq = Zj if j £ Jo = J(0) and oo otherwise. (2.23) 

On [Rn, Rn+i), the X^ : j G J„, follow independent copies of B^ starting at Xi . 

Jin 

At Rn+1 we define 

'Xi if JG J„ = J(i?„), 

7ln+l — 

ifjGJ„+i-J„, 

Hn+l — 

OO otherwise. 

Note that offspring are no longer displaced from their parents and that coalescence 
reduces the number of particles born at time Rn+i, but otherwise no coalescence 
occurs as J{t) = Jn on [Rn, Rn+i)- Thus, conditional on the sequence {vr„}, X is 
a branching random walk starting with particles at zj, j G Jo(vro), with particle 
branching rate c* and giving birth to |7r„| — 1 particles on top of the parent Xg," 

Rn 

(who also survives) at the ?ith branch time Rn- 



tin+l 
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2.6 Computation process ( 

As we did for X, for tQ E [0,T) we now define an computation process {Cr(^) '■ < 
k < k{{T - r)-),r G [to,T]} for X. Here ^(s) = M + mNo if s G [^^,^^+1). 
Given are the branching random walks {Xs, s G [0, T — tg]}) the associated sequence 
{■Kn, Rn, P-n}, & Sequence of iid random variables {Un}, uniformly distributed on 
[0,1] and independent of {X , {-Kn, Rn, fj-n ■ n- G N}), and a set of initial values 
Cto(j),i G J{T — to). In the next section when we couple {XX) and {XX) we 
will set Un equal to C/„ defined in (2.11). Define an equivalence relation on 
{0,...,M + nNo}hY 

M+{m- l)No + j ^f^^ M + {m- l)No + i {I < i, j < Nq, I < m < n) iS j i, 
M+{m- l)No + j Am (1 < J < iVo, 1 < m < n) iff j 0, (2.24) 

J-rJ: {0<i,j<M) iSj^^.i. 

Finally if ^„ < t < Rn+i define i ss* j iff i j for < i, j < M + niVo- To 
prepare for the proof of Lemma 2.10, note that the definition of ^ that follows is 
just the definition of ^ with hats added and ~t used in place of 

First we complete the initial state Qo by setting Cto(^) = CtoO) if ^ ~T~tQ j S 
j(T — io), k < K{T — to) = K{{T — to) — ) a.s. Suppose that for some ?Ti ^ 1, Rm is 
the largest branch time smaller than T — to- The values Cr{k) do not change except 
at times T — i?„, so Cr = Go for t < T — Rm- We decide whether or not to flip the 
value of C at fim at time t — i?^ as follows. Define S {0, 1}^° by 

Vl = l^^_^^^_{M + {m-l)No+j), j = l,...,iVo. (2.25) 
Letting i = C,(t~r )-(/*m) we set 

^ ' \i otherwise. 

To update C, now, for k < M + {m — l)A'^o and k 7^ fim, 

^ ^Rm 1^'^ '^T-R^^) = '5^(T-ij™)(A™) ' (2.27) 
and for the remaining values of A; < Af + {m — l)A''o keep 

The values C,r{k) remain constant for r G [T — Rm,T — Rm-i)- Coming to 
r = T — Rm-i, if m — 1 > 1 we proceed as above. When we reach r = T — Ro = T 
we end by setting (t = Ct- ■ 
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2.7 Coupling of (X, C) and {X, () 

We now give a construction oi X,(^ which will have the property that with high 
probability for small e, (i) X and X are close and (ii) given identical inputs, C and 
( will compute the same result. As before, T > and z = {zq, . . . , zm), Zi £ e^'^ are 
fixed. Recall the reaction times Rm, the uniform random variables Um from (2.11), 
and the natural time-reversed filtration J-f used in the construction of the dual X 
given in (2.10). 

The following general definition will be used to construct the partitions {7r„ : 
n S needed to define X, (distributed as in Section 2.5) in terms of the graph- 
ical representation. Let V be an J^j-stopping time (think V = Rm), and let 
70 ... , 7a/' G eZ'^ G eZ^ be J"y-measurable. Let {S^'T' : i = 0, . . . , M'} C eZ"^ 
be the rescaled coalescing random walk system, starting at time V at locations 
7o,...,7m'5 determined by the {T^} in the graphical representation. That is, 
{B^''^* : i = 0, . . . , M'} are as described in Figure 6 but now starting at time T — V 
at sites 70 , . . . , 7^^ . For each t > let vry^^ (t) G Hm' be the random partition of 
{0, . . . , M'} associated with the equivalence relation i ~ i' iff B^''^'- (t) = B^''^*' (t) . We 
call TTv^-yit) the random partition at time V+t with initial condition 7 = (70, . . . , 7m') 
at time V. 

Let ttq = vro,^(-y/e) G Um be the random partition of {0, . . . , M} at time ^/£ with 
initial condition z = {zq, . . . , zm) at time 0, and note that its law is the same as the 
law of the ttq described just before (2.21). For m > 1 let 

and {7r[^,m G N} be an iid sequence with law i/^ and chosen independent of J-'oo. 
For m G N, define 

/^Rm,7m(\/^) a Rn> Rn^i + Ve iov alll <n<m 
I vr^ otherwise. 

By the translation invariance and independent increments properties of the Poisson 
point processes used in the graphical representation and also (2.15), is indepen- 
dent of V(7(7r^, n < m) = Fm-i, and has law Ve defined just before (2.21). 
It is also easy to check that vr^ is ^m-measurable (m > 0) and so {vTm, m > 0} are 
independent and distributed as in Section 2.5. 
For m G N let 

= inf{s > Rm-i : 3z / j both in J{Rm-i-), or 

i G J{Rm~i-) \ {Mm-i},i G J{Rm~i) \ J{Rm^i-), SO that Xi = Xi}, 
Tm = mf{s > Rm^i + y/^ : inf \Xi - X^ < e'/^}, 

Y:^ = max{\Y^\:i = l,...,No}. 



We introduce the time, Tb, that one of four possible "bad events" occurs: 
Tb = inm{Rm : m> l,Rm. < Rm~i + \/e or > - log(l/e)} 

K 

A min{rm : m > 1, < -Rm} A min{T^ : m > 2, < Rm-i + Ve}- 
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Here min0 = oo To see why the last two minima should be large, note that after 
a birth of A^o particles from particle Hm at time Rm, we expect some coalescence 
to occur between the parent and its children. After time ^/e, particles should all 
be separated by at least e^/^ and remain that way until the next reaction time 
when again there may be coalescing within the family producing offspring but no 
other coalescing events. The qualifier m > 2 is needed in the last minimum because 
we have no control over the spacings between particles at time 0. The collision 
of particles 2 and 7 in Figure 6 is an example of a bad event that enters into the 
definition of t^. We assume throughout that 

< e < ei(K) so that - log(l/e) < (2.29) 

Given {vTm} we now construct X and A{s) = {{fin-, Rn)^{Rn < ■5))n6N (with the 
law described in Section 2.5) initially up to time T = T{, A Tfy, where 

fh = mm{Rm : m>l,Rm.- Rm~i < V^}- 

Once one of the five bad events (implicit in the definition of T) occurs, we will give 
up and continue the definition of the branching random walk using independent 
information. The coupling of X and X will be through our definition of {vr^} and 
also through the use of the random walks steps of X^ to define corresponding random 
walk steps in X^ whenever possible, as will be described below. 

We begin our inductive construction by setting Rq = 0, j(0) = Jo(vro), and 
define Xq as in (2.23). Note that 

J(0) = J(Ve) = Mt^o) if ^1 > Ve. (2.30) 

Assume now that {X, A) has been defined on [0, Rm AT]. Assume also that Rm < T 
implies the following for all 1 < i < m: 

Ri = Ri, jj-i = /^i, (2.31) 

J{Ri) = J{R^-l) yj{M + {i- l)No Jo(7r,) \ {0}}. (2.32) 

J{Ri_i) = J(s) C J(s) for all s G [Ri^i,Ri). (2.33) 

J{s) = J{s) = J{Ri^i + Ve) for all s G [R^-l + ^e, i?»), (2.34) 

The m = case of the induction is slightly different, due for example to 
the special nature of vro, so let us assume m > 1 first. To define (X,yl) on 
(Rm A T,Rm.+i A T] we may assume Rmi^) < T{oj) and so (2.31)-(2.34) hold by 
induction. On {Rm,{Rm + \/e) A A T] let (-^,^4) evolve as in Section 2.5 

conditionally independent of Jqo given {vr^}. Here it is understood that the unused 
partitions {ttj : i > m} are used to define the successive branching events as in 
(2.22). 
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Next, to define {X, A) on {{Rm + V^) A Rm+i A T, Rm+i A T] we may assume 
Rm{oj) + -^/e < Rm+i A T{uj). By the definition of TJ, this imphes Rm+i > Rm + \/e 
and so for all s G [i?mi Rm + V^]) 

J(s) = J{Rm) = J{Rm-i) U {M + (m - l)iVo + J : J G Jo(vr„,) \ {0}} 

= J{Rm-i + V^)U {M +{m- l)No e Jo(^m) \ {0}}. 

(2.35) 

In the first equality we used (2.32) and in the second we used (2.33) and (2.34) with 
s = Rm-i + V^- The fact that < Rm (since Tf, > Rm + -y/e) shows there are no 
coalescings of X on [Rm-i + V^, Rm) and so 

J{Rm.^i + V^) = JiR^-). (2.36) 

Again use Tb > Rm + together with (2.29) to see that y„* < f log(l/e) < 
and so the spacings of the previously existing particles at time Rm < Tm ensures 
that none of the A^o new particles at time Rm will land on a previously occupied 
site. Therefore if 

MYm) = {1 < i < Ao : ^ {y^ : < i < j}}, 

then 

JiRm) = J{Rm-) U {M + (m - l)Ao + J : j G Ji{Ym)]. 

The fact that Rm+i A t^+i > Rm + means that X has no branching events in 
{Rm,Rm + and X has no particles coalescing on [Rm,Rm + except those 
involving + Y^ , i = 0, . . . , Aq. Therefore, the definition of iTm ensures that 

J{Rm + V^) = J{Rm-) U {M + (m - l)Ao + J : i G Jo(vr™) \ {0}}. 

= J(s) for all s G [i?^, iJ^ + (2.37) 

where in the last line we have used (2.35) and (2.36). For s G [Rm + Rm+i A T) 
we have s < Tm+i and so 

|Af -A,^| > e^/*^ for all j / /c both in J(s), for all s G [Rm+V^, Rm+i^f). (2.38) 

In particular X can have no coalescings on the above interval and so J{s) = J{Rm + 
^) for se[Rm + Ve, ^m+i AT). On {Rm + V^, Rm+i A f] let (a| , j G J(s)) follow 
the random walk steps and branching events of {X^ : j G J{s)} (of course there is 
at most one of the latter at time Rm+i providing Rm+i ^ T). In particular we are 
setting 

J{s) = J{s) = J{Rm + Ve) for s G [i?™ + Ve, Rm+i AT) or s = f < Rm+i- (2.39) 

(2.38) shows that the random walk steps and branching events for distinct particles 
of X on (i?m+-y/e, Rm+iAT] are independent. In addition, these steps and branching 
events are independent of the random walk increments used to define {tTu}- This 
shows that X evolves like the branching random walk described in Section 2.5 on 
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{Rm,Rm+i^T), and on {Rm,Rm+i/\T] if either T < Rm+i, or Rm + Ve > Rm+i/\T. 
(In the latter case the first part of the above construction did the job and in the 
former case there is no reaction event to define at T A Rm+i =T.) So to complete 
the construction at i = Rm+i /\T we may assume 

Rm + V^< Rm+l < f. (2.40) 

The above definition shows that Rm+i = Rm+i, we use (2.32) with i = m + 1 to 
define J{Rm+i) and we set fim+i = fJ-m+i- Clearly fim+i is uniform on J{Rm) = 
J{Rm + \/e) (given {'Kn}) and is independent of {fin '■ n < m}. In addition the 
branching events used to define {fin} are independent of the random walk steps used 
to define {7r„}. This completes our inductive definition of {X, A) on [0, Rm+i A T]. 

Next we complete the inductive step of the derivation of (2.31)-(2.34) for m + 1 
under (2.40) which is in fact weaker than the Rm+i < T condition. (2.39) implies 
(2.34) for i = m + 1, and (2.31) and (2.32) hold by definition. On {R^ + y/e < 
Rm+i < T} J can only decrease on [Rm, Rm + due to coalescings of the random 
walks, while J is constant on this interval by (2.37). The inclusion (2.33) therefore 
follows from the equality in (2.34). 

To complete the inductive construction of (X, A) on each [0, Rm A T] and proof 
of (2.31)=(2.34) it remains to give the m = step of the construction and verify 
the m = 1 case of the induction. Both follow by making only minor changes in 
the above induction step. For example, (2.30) is used in place of the (now non- 
existent) induction hypothesis (2.32) both in defining X on the initial interval and 
in obtaining (2.37) for m = 0. 

Since Rm t oo a.s. we have defined {X^A){s) on [0, T] and to complete the 
definition we let it evolve conditionally independently (given {vTn}) for s > T. 

The above construction and (2.16) show that 

{X, {7r„}, {fin}, X, {jln}, {Rn}) is independent of {J7„}, (2.41) 

where {Un} are the uniforms from (2.11). Therefore the computation process C for 
the above X may be defined as in Section 2.6 but with Un = Un- 

Lemma 2.1. (a) For all m G Z+, Rm < Ti, and Rm < Tfe imply Rm = Rm < T. 
(h) For all m eN, if 

Gm = {a; ■.[^ZlR^ - Ri-i) A (a™+V; - > V^, 

Ri <TiMi< m,maxyj* < - log(l/e), A™ - Ri^i > ^/e}, 

i<m K 

then Gm C {Rm. = Rm < T}. 

Proof, (a) The implication is trivial for m = so assume it for m and assume 
also Rm+i < Ti), Rm+i < Tf,. By induction we have Rm = Rm < T. Since 
Rm+i A Rm+i > Rm + V^, wc also know T > Rm + "v/e. The construction of X on 
{Rm + Rm+i A T] shows that the next reaction time of X on this interval must 
be Rm+i (if it exists) and so Tb > Rm+i- Since Tfe > Rm.+i by hypothesis we get 
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{Rm + \/e,Rm+i AT] = [Rm + \/^-,Rm.+i]- Hence our construction of X on this 
interval shows Rm+i = Rm+i and so the result follows for m + 1. 
(b) The first four conditions in the definition of imply 

Th > Rm+l A A Tm+l > Rm- 

The last condition implies TJ, > Rm- Now apply (a). 

□ 

As an immediate consequence of the above and our inductive proof of (2.31)- 
(2.34) we get the following: 

Lemma 2.2. 

Rm <f ^ for alll<i<m (2.31) - (2.34) hold. 

On Gm and on the intervals [Rm-l+^/£, Rm) our definition of X and Lemma 2.1(b) 
shows that the movement of particles in X and X are coupled (they take identical 
steps) but on [R„i~i,Rm-i + V^) they move independently. To bound the discrep- 
ancies that accumulate during these intervals we use: 

Lemma 2.3. //w G Gm., then 

snv{\Xi -Xi\:je J{s), s G [0,Rm)} (2.42) 

m— 1 

< _ 1) £ iog(l/e) + sup \Xi-Xl\ + \Xi - X|jJ . 

Proof. Suppose first that m > 1 and we are on Gm- By the coupling of the spatial 
motions noted above, for j G J{Rm-i) 

sup \Xi-Xi\ = 

S&[Rm-l,Rm) ; 

< \xi -xi I + 

' — 1 -'I'm— 1' 

+ 

On Gm-, a newly born particle to X-'^^__^_ may jump a distance at most | log(l/e) 

from its parent, while for X]j ^_ it will be born on its parent site, so the above is 
at most 

sup |X|j_^_-X^_^„| + ilog(l/.) 
feeJ(i?^-i-) 

+ sup \Xl-X^^^J+ sup |X^--X^_J. 

Things are simpler when m = 1 because there are no initial jumps to worry about 
and so the second term in the above is absent. The required bound now follows by 
induction in m and the fact that Gm is decreasing in m. □ 
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2.8 Bounding the probability of bad events 

Here and in what follows it is useful to dominate X with a branching random walk 
X, also with paths in T) and with the same initial state. Particles in X follow 
independent copies of and with rate c* give birth to A'^o particles located at 
Bl + y^, i = 1, . . . ^Nq, where Bf is the location of the parent particle. At the 
mth birth time Rm we use 

XM+(m^i)No+i ^ i = 1, . . . , iVo to label the new particles, 
so that if J{t) = {j : Xl ^ oo}, then J{Rm) = {0, . . . ,M + rriNo}. Coalescence 
is avoided in X by having the coalescing particle with the larger index have its 
future steps and branching events dictated by an independent copy of the graphical 
representation. This will ensure that J{t) C J(t) and {X^it) : j £ J{t)} C {X^{t) : 
j G J(t)} for ah t > 0. 

Let Nt = min{m : Rm. > T} and define Nt in the same way, using the branching 
times {Rm}- Let 

Cb = c*No > 1 . (2.43) 
We will also need to separate the particles in X and so define 

fm = inf{s > Rm~i + \/e : inf \Xi - Xi\ < e^/^}, m G N. (2.44) 

Lemma 2.4. There is a constant C2.4 so that for all T > and n G N 

(a) P{Nt >n)< P{Nt > n) < e^t^(M + l)(niVo)"^- 

(b) P{mmi<m<NT Pm - Rm~l < y/s Or mini<m<ArT Pm - Rm-1 < V^) 

< ca,^e^f^(M + l)eV6. 

Proof, (a) The first inequality follows from the domination of X by X. For the 
second one note that E{J{T)) = (M + l)e'^''^ and conclude 

P{Nt >n)< P{\J{T)\ > M + 1 + niVo) 

< (M + 1 + nNQ)-^{M + 1)6'="'^. 

(b) Let Z be a mean one exponential random variable. The domination of X by 
X shows that for any n > 1, 

P ( min Rm - Rm-i < \^] < P ( min Rm - Rm~i < 

\l<m<NT J \l<m<NT 

n „ 

m.=l 

n 

< e^''^{M + l)(niVo)^^ + ^ (M + 1 + (m - l)A^o)\/^, 

m=l 

by (a). Now set n = [e~"^/^] and note that the sum is at most {M +l)n^/e+v? NQy/e . 
A similar calculation gives the same upper bound for the Rm^s. □ 

Lemma 2.5. There is a constant C2.5 so that for all T > 

P{Y^ > - log(l/e) for some m < Nt) < cs.se'^'^iM + l)e^/2 
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Proof. By (a) of Lemma 2.4, P{Nt > n) < e''>''^{M + l)(n7Vo)~^. Using (1.8) gives 
P{Y^ > — log(l/e) for some m < n) < nCe. 

Taking n = [e~^/^] now gives the desired result. □ 

The following facts about random walks will be used frequently. 

Lemma 2.6. Let Zs denote a continuous time rate 2 random walk on jumping 
with kernel p, and starting at x ^ under P^ , and B'^ he our continuous time 
rescaled copy of Z , starting at z ^ eU^ under P^. 
(a) For any to ^ 0, tq > 1, x G Z'* and p > 2, 



P^'ilZsl < ro for some s > to) < C2,6 / \[{\x\-ro)+]-P{sP/\s)\A\{s\/l)-'^/^r^ 



to 



ds. 



(h) sup^P^dZ^I < for some s > e"3/2) < c^^^e3/8_ 

(c) For any z G eZ"', tq > 1 

Pz{\B's\ < roe for some s > 0) < c^,g(|z|e^i)-(2/3){<i-2)^2(d+i)/3_ 

Proof, (a) Use T{tQ, y) < oo to denote the time of the first visit of Z to y after time 
to, and let 

poo 

G= / P°{Z, = 0)ds 







be the expected time at (which is finite since d > 3). Then 

roo , roo 

/ P^(Z, = y)ds = EHl{T{to,y){u) < oo} / P^(^._T(to,j/)H = y)ds 

= GP^{Tito,y)<oo). 
Summing over \y\ < r^ for > 1 and rearranging, we get 

/•CO 

P^(|Z,| < ro for some s > to) < G'^ ^ / P^(Z, = y) ds 

|y|<ro"^*° 

POO 

= G-^ / P^(|Z,| < ro)ds. (2.45) 
Jto 

A martingale square function inequality shows that for p >2, 

P'd^.l < ro) < P°(|Z,| > {\x\ - ro)+) < c((|x| - ro)+)-n(s^/' V s). (2.46) 
A local central limit theorem (see, e.g. (A. 7) in [6]) shows that 

P^'ilZ.l < ro) < c{s V l)^'^/^^. (2.47) 
Use the above two inequalities to bound the integrand in (2.45) and derive (a). 
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(b) Set ro = e"'^^^ and to = e~^^^ in (a) and use only the second term in the infimum 
inside the integral. The right-hand side is ce"^'^/*^"^^/^^"''^^'^/^^ . To complete the 
proof we note that exponent is smallest when (i = 3. 

(c) We may assume without loss of generality that ro < \z\£~'^ /2 = M/2 (or the 
bound is trivial) and so ti = M'^l'^r^'^ > 1. Apply (a) with p = 2d and break the 
integral at ti to see that the probability in (c) is 

-P'^'' {\Zs\< ro for some s > 0) < c / M-'^'^{s'^ V s)ds + / g-'^l'^rf^ds 

'-Jo Jti 

<c(M-2'^tf^ + t}-('^/'Vo") 

< cM-(2/3)('^-2)^2{d+l)/3^ 



□ 

Lemma 2.7. P{Tm < Rm or fm < Rm for some I < m < Nt) 

< C2.7e'=''^(M+l)2e3/32. 

Proof. To bound P{Tm < Rm for some 1 < m < Nt), we start with 

P(,Tm < Rn^J'Rrr.-i) < P{R,n > Rrn-1 + V^, 3i / j both in J{Rm-l + V^), S.t. 

inf \Xl-Xi\<e'/'\J^n^_,). 

\/e + Rm-l<S<Rm 

Now i 7^ j both in J{Rm-i + \/e) and Rm > Rm-i + a/e imply i, j € J{Rm-i) and 
XI / Xi for all s € [Rm^i,Rm~i + -v/e]- Therefore, the above is at most 

P{Xl - xi / 0, Vs G + VF], (2.48) 

\Xi - Xi\ < 35 > + 

If Z as in Lemma 2.6, we may use (b) of that result to bound the above by 

\J{Rm-l)\^ sup P^»(|Z,| < 8-^/^ 3S > £^3/2) 

< (M + 1 + (m - l)iVo)2 • ce^/* 
Using Lemma 2.4(a), we conclude 

PiTm < Rm for some 1 < m < Nt) 

n 

< e^'^^iM + l)(nA^o)"^ + ^ (M + 1 + (m - l)iVo)^ce^^^ 

m=l 

To bound the sum we note that for a, 6 > 1, 

V(a + (m- 1)6)2 < /■ (^^2.5)2^^ ^ Ij^^^^^^s _^3] < ^^2^^^^3^ 
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Taking n = [e"3/32/^^^ 

gives the desired bound. A similar calculation (in fact there 
is some simplification) gives the same upper bound for 

Pi'T'm < Rm for some 1 < m < Nt-) 

□ 

Lemma 2.8. 

P{ min T^^^-Rm<V^)<C2.8e'''^{M + lfe^/^''. 

l<m<Nj' 

Proof. Define 5m D Gm as Gm (in Lemma 2.1) but without the lower bounds on 
A^^V| - Ri-i or A^^i^i - Ri-i. Note that Sm e J^r^ and if w e Sm, then 

- Xj^j > e'^/^ for all distinct i,j G J{Rm-)- (2.49) 
In addition, since < ■|log(l/e) we have that for all i E J{Rm—) — {/Um}, 3 £ 

J{Rm) — J{Rm-) 

\X'^^ - Xj,^ I > e'/' - - log(l/e) > e"'/2 (2.50) 

since e < ei(K) (recall (2.29)). 

If Tq is the return time to zero of the random walk Z in Lemma 2.6, we have 
(see P 26.2 in [ 11] for d = 2> and project down for d > 3) 

P'"'{Tq <oo) <c\zq\~^. (2.51) 

Use (2.49), (2.50), and (2.51) with scaling, and the bound 

\J{Rm-)\ <M + l + {m-l)No 

to see that on Sm G -^i?„i) 

P{Tm+i-Rm<yfe\^Rm) 

< c[{M + 1 + (m - l)iVo)2eV8 + (m + 1 + (m - l)No)Noe^/^] 

< c{M + ifni^N^e^^^. (2.52) 
Using the bound in Lemma 2.4(a), we conclude 

yl<m.<Nx J 
n 

< P{Nt > n) + V P{Nt = k,Sk, min t^^, - Rm < V^) 

^ — ' l<m<k 
k=l - - 

n k 

< C2Ae'''^{M + l)(niVo)-i + Wi -Rm<V^)- 

k=l m=l 

Using (2.52) now, the above sum is at most 



cN^{M + ifne^/^ ^ < c(M + ifn^e^ 



m=l 

Take n = \e~^/^^^ and use Lemmas 2.4, 2.5, and 2.7 to bound P{Sj^^) to get the 
desired result. □ 
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2.9 When nothing bad happens, (X, () and {X, () are close 

The next result gives a concrete bound on the difference between X and X and 
deals with the final interval [Rm, Rm A T]. Let 

G,n = Gmr\{f> R^}, men, 

and for < /? < 1/2, define 

= Gnt n {sup sup \Xi - Xi\ < e^/^} 

n{r^u^Lo'[^m,i?m + 2e^]}. (2.53) 

Allowing smaller /3 values will be useful in Sections 3 and 4, but for now the reader 
may take (3 = 1/2. 

Lemma 2.9. There is a C2.9 and £2.9{i^) > so that for any T > 2e^ , < e < 



1 a£ 



40 ' " 3 , 



On we have J{s) = J{s) for all s £ [T - e'^,T], and - > e'^/^ for all 
i ^ j in J{T) . 

Proof. Dependence on /3 will be suppressed. For s as above, Lemma 2.2 implies 
J{s) = J{s) on Gt since s G [Rn{t)-i + V^iRn(t)) Gt- The last assertion of 
the Lemma holds on Gt because on Gt, t^^^t) ^ ^n{T) 

T G [i?jv(r)-i + V^,Rn(t)) = [Rn(t)-i + V^,Rn{T))- 
Lemmas 2.4, 2.5, 2.7, and 2.8 imply 

PiG'N^) < ce^'^iM + l)2eV40_ (2.54) 
To deal with the first additional good event in Gt, we note that by Lemma 2.3 
PIGnt, Slip sup \Xi - Xi\ > e^/6) < P{Nt > n) 



s<T , 



n-1 

£ 



+ P((n- l)-log(l/e) + ^ sup sup \Xi - Xj^J + \Xi - xj^J > e^/^ 

By (a) in Lemma 2.4 the first term is at most e'^'''^{M + l)(nA'^o) If 

(n-l)-log(l/e) <ei/V2, (2.55) 

then it enough to bound 

^"-1 ,1/6 \ 



£ 



P ^ sup sup \Xl-X\^^\^\Xi-X\^^> 



i=0 ■ 



< J2^M + 1 + iNo)2P(^ sup \BI\ > —j < c{M + l)n^No ■ n^£-^/^£^ 

i=0 s<v^ 
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by the maximal inequality for martingales. If n = \e ^/^^^^ (so that (2.55) holds 
for e < e2.9{n)) the above gives 

pIgnt, sup sup \Xi -Xi\> e^/^ ) < ce^''^(M + l)e^/^\ (2.56) 

The domination of X by X ensures that 

Therefore (recall T > 2e^) for any ^ G N 

P{T eiJ^iro'[Rm,Rm + 2e^]) 

< P{Nt >£)+ P{T G Ut=\[i?™, Rn. + 2e^]). 

Lemma 2.4(a) shows that the first term is at most e'^^'^ {M + \)[iNQ)~'^ . Conditional 
on J-^^_^i Rm — R-m-i is an exponential random variable with rate (M + 1 + (m — 
1)Nq)c* , so the second term is at most 

e(j2p{T- 2e^ - Rm^i <Rm- Rm-i <T - 

\m=l / 

I 

< 2el^ ^ ((M + 1 + (m - l)iVo)c*) < ce"*''^ {M + 1)^8^^. 

m=l 

Taking I = [e-'^/^] then using (2.54) and (2.56) gives the desired bound on P{G^). 
□ 

The next ingredient required for the convergence theorem is: 

Lemma 2.10. Assume T > 2e^, to G [0,e^], and u e G^. If CtoU) = CtoU) for all 
j G J{T-to), then Ct(«) = Ct(«), i = 0, . . . ,M. In particular if CtoU) = Ctoi^T-to) 
for j G J(T - to), then Ct(«) = S.T{zi) for i = 0, . . . , M. 

Remark 2.1. By Lemma 2.9, J{T — to) = J(T — to) on G^, and so all the necessary 
inputs required for both computations are prescribed in the above result. 

Proof. The last statement is immediate from the first and (2.17) with r = T. 

By the definition of Gn^, D and Lemma 2.2 there is a unique n < Nt so that 

Rn + V^<T-e'^<t-to<T< Rn+i (2.57) 

and 

Rm = Rm and = for m < n + 1, K[s) = K{s) for s G [0, T]. (2.58) 

As was noted in Section 2.6 the inductive definitions of (" and <^ are identical except 
the latter has hats on the relevant variables and uses ~j in place of ~t. The above 
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shows that in our current setting the relevant variables are the same with or without 
hats (recall we are using Un = Un in our coupled construction of C) and so it remains 
to show the equivalence relations are the same and we do this now for the initial 
extensions. That is, we extended Cto to {0, . . . , K{T — to)} by Cto(^) — CtoO) 
k ~T-to J £ ~ *o) = — io) (see the above Remark) and extended Cto in 
the same way but if A; ~ j G J(T — Iq) which means ^T-to ~ ^T-to' 
now show these equivalencies are the same and hence so are the extensions. Note 
that in applying (2.24) to extend Qo "^e are using VTm = '^Rm,7ni,{V^) for m > and 
'''"0 = '^o,z{\/^)- This means two indices j, A: in a family which has branched at time 
Rm = Rm, 0<m<n(ifm = this means two initial indices) are equivalent (in the 
~ sense) at time T — if their corresponding X paths coalesce by time Rm + 
Lemma 2.2 implies that on there are no coalescing events in [0, T — to] (ii^ fact 
on [0, T]) except for those in [Rm,Rm + involving a common family born at 
Rm, for m < n. Therefore, the above condition is equivalent to ^^^-to ~ -^T-to 
the required result is proved. 

The Lemma now follows easily by induction up the tree of X. In place of the 
above we must show equivalence of the equivalencies used in (2.20) and (2.27) at 
times Rm- Note here that for the indices of interest in (2.20) and (2.27) this is 
equivalent to the corresponding equivalencies at times Rm~i + and this follows 
as above for m > 1. □ 



2.10 The branching Brownian motion and computation process 

We now define a branching Brownian motion starting at x G R'^ with paths in 
"D. Let {vrj^,n > 1} be an iid sequence of partitions with law uq (defined in the 
second paragraph of Section 2.5). Particles in X^ branch at rate c* and at the nth 
branching time, \tt^\ — 1 particles are born at the location of the parent who also 
remains alive. After birth, particles in X^ move as independent Brownian motions 
in R"^ with variance parameter o"^. To couple X^ with the branching random walk 
X^ from Section 2.5 we need two preliminary lemmas which allow us to couple 
the corresponding particle motions and offspring numbers, respectively, of the two 
branching processes. 

Lemma 2.11. We may define our scaled random walk and a d- dimensional 
Brownian motion B with variance a'^, starting at 0, on the same space so that for 
some constant C2.n 

P ^sup \Bf-Bt\>VI^ < cs.nTe. 

Proof. Apply Theorem 2.3(i) of Chapter 1 of [10] with H{x) = x^, to see we may 
define the unsealed random walk B^ (rate 1, step distribution p) and a Brownian 
motion as above, B' on the same space so that for all S > and r > 1 

P I sup \Bl -B',\>r] < cSr-^. (2.59) 
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Although the above reference applies to discrete time random walks, we apply it to 
the step distribution Xlill'i'' where {Xi} are iid p{-) and A^(l) is an independent 
Poisson(l) random variable. We arrive at the above after a short interpolation 
calculation for B^. 

To get the desired result from (2.59) we set Bf = eB^_2^, Bt = eB'^_2^ and use 
r = e~^l'^ to conclude that 

P(sup \Bl - Bt\ > Ve) < -P(sup \eBl-2t - eB[^2t\ > Ve) 

t<T t<T 

< P{ sup \Bl - B[\ > £-1/2) 

t<e-^T 

< ce-'^Te^ = cTe 

which proves the desired result. □ 

Lemma 2.12. For each e > we may construct the sequence {ir^ : n > 1} on the 
same space as {vr^ : n > 1} so that 

P«/vr0)<c,.,,e3/4. 

Proof. The obvious way to couple vr^ and ir^ is to use the same system of rate one 
coalescing random walks {B^^ : i = 0, . . . , Nq}. If Z is as in Lemma 2.6, then by 
(2.45) and (2.47) 

P(7r^ / 7r°) < supP^(Z, = for some s > e^^/^) 

X 

□ 

Let Xs G for e > and assume — )• rr G M'^. Our goal now is a joint 
construction of (X'^jX^) started from {xe,x), and associated computation processes 
(C"; C^) with the property that if C*", have the same inputs then they will have the 
same outputs with probability close to one. 

The branching random walk X^ starting with a single particle at x^, along with 
the associated index sets J^(-), branch times {Rm}^ and parent variables /l'^}, are 
constructed as in Section 2.5 using the sequence {vr^} in Lemma 2.12. There is 
no initial coalescing step now as we are starting with a single particle. We use the 
coupled sequence {vr^ : n > 1} to define the offspring numbers, branching times 
{iZj^ : n > 1}, index sets J^{-) and parent variables {/x^ : n > 1} with the same 
conditional laws (given {ttJ^}) as in the definition of X^. We may couple these two 
constructions so that for all n € on the set 

= {T^m = ^h < m < n} , 

we have 

Rf^ = /i^ = /^o^, and J°(s) = J^(s) for all s < for aU m<n. (2.60) 
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Define = inf{m : i?^ > t}. Using these sequences we follow the prescription in 
Section 2.5 for constructing X but substituting Brownian motion paths for random 
walk paths. Couple these random walks and Brownian motions as in Lemma 2.11 
at least as long as the branching structure of the two are the same. Note that if 
there are n branching events up to time T there are at most 1 + uNq independent 
random walk segments and Brownian motions of length at most T to couple (recall 
our labeling scheme from Section 2.5). In addition to the errors in Lemma 2.11 
there will be a small error from the difference in initial positions at time 0, and so 
we get 

P{G% , sup sup - Xl'^l > \xe -x\ + ,/e) (2.61) 

< P{N^ > n) + (1 + nNo)c2.nTe + C2,i2^e^^^. 

The first time vr^ ^ vr^ we declare the coupling a failure and complete the 
definition of for t > iiJl using random variables independent of X"^ . 

Fix T > and to G [0,T). Given X^,J^{t),0 < t < T, the sequences {tt^}, 
{Rn}, {fJ-n} ^u independent sequence of iid uniform [0,1] random variables {U^^}, 
and initial inputs {CtoU) '■ 3 £ J^iT — to)}) we define a computation process C?, to < 
t < T. The definition is analogous to that of Qt given in Subsection 2.6 for X'^ 
started at a single point, but we use gi-i in place of g\_^ in (2.26). That is, as in 
(2.25), (2.26), we have 

V^' = C(°T-Roj_((m - l)iVo + j), j = 1, . . . , iVo, 
and if 2 = C(7-„^o )„(Mm)) we have 

<lr,K)A'-' "f^^'-*''"'/^' (2.62) 
^ \i otherwise. 

We further couple CP and ("^ by using the same sequence of independent uniforms: 
{^m\ — {Um} in their inductive definitions. Just as in (2.41) we can show that 
this sequence is independent of all the other variables used to define X^ and Q', as 
required. We let JF^ denote the right-continuous filtration generated by X^ , X^ and 
A^{t) = ((i?^^, vrj^, Um)^{Rm < i))m.<^N as well as its counterpart for X^ . 

Notation. = G^o n {sup^<7. supjgjo(^) - Xf^\ < \x^ - x\ + 

G"/ = G^/n[Um^ for alUG {0,l}^«,m<ArO,i = 0,l}. 

Lemma 2.13. (a) On G^"^ , we have 

Rf^ = /z^ = fi^^jTrf^ = TT^, for all m < Nj^, and J'^{s) = J^{s) for all s <T. 
(h) P{{&^y)< 02,13^"'^^'" ■ 

(c) On G^rjf we also have for any to e [0, T), if CtoU) = CtoU) for all j G J^{T - to), 
then (0(0) = Ct{0). 



(d) P{{G'/)^) < cgjse^"^ e^/' + J^lo M " 9 
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Proof, (a) is immediate from (2.60) and the definition of G^'^. 

(b) follows from (2.61) and the now familiar bound P{Nj. > n) < by setting 
n=re-3/8^. 

(c) On G^^, we see from (a) and the inductive definitions of and C^, that all the 
variables used to define Cj'(O) and Cf'(O) coincide. Therefore these outputs can only 
differ due to the use of gi-^i in (2.62) and the use of gf_^ in (2.26). By induction we 
may assume Vm = additional condition defining now ensures that 
these two steps produce the same outputs. 

(d) The additional condition defining g!^^ fails with probability at most (recall 
c* > 1) 



i 

PiNT>n) + n2''^^[Y,\\9^-g^\\oo 

1=0 



< 



nNr 



+ n2 



No 



i=0 



\9i 9i lloo 



Now let n = I^^I^q \\gi — fff ||ool and use (b) to complete the proof. 



□ 
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3 Proofs of Theorems 1.2 and 1.3 
3.1 Proof of Theorem 1.2 

We start with a key estimate giving the product structure in Theorem 1.2. This 
relies on the fact that duals starting at distant points with high probability will not 
collide. For some results we will need a quantitative estimate. Let E eZ*^, iji G 
and = Zjfc = -\- eyi, for < i < L and 1 < k < K. Set 

Ae= min \zik - Zi>k'\e~^ ■ 

l<i,i'<L,l<k^k'<K 

The notation is taken to parallel that in Theorem 1.2 and hypothesis (1.23) of that 
result implies 

limAe = oo. (3.1) 

e 

Let X = X^'^ be the dual process starting at z for the time period [0,T], with 
associated computation process Ct which has initial inputs Co(j) = Co(^t)' ^ •^i'^)- 
Let = {zik,i = 0, . . . , L) and consider the duals X^''^'^, 1 < k < K defined as in 
Section 2 with their associated uniforms {Uj^} and parent variables {/^^}- These 
duals are naturally embedded in X^'^, and although the numbering of the particles 
may differ, we do have 

{X^'^ : j e J{t)} = uf=i{xf : j e r^m, t e [0,T]. (3.2) 

Define 

V,,T,e = mf{i G [0, T] -.X^"'^^^ = X^"'^^^'' (3.3) 

for some 1 < k ^ k' < K,j £ ■P''{t),j' G J^'='(t)}, 

where inf = oo. 

Lemma 3.1. P(l^^,T,e < oo) < C5__^(iv:,L)e'^''^(A£)-('^-2)/('^+3). 

Proof. We may dominate X^'"'^ by the branching random walks X^^yT fj-gm Sec- 
tion 2.8. By Lemma 2.4(a), if {Y^^} are iid, equal in law to Y*, and independent of 
in what follows, then for R> 1, 

P{V < oo) 

< P(max N^''- > n) 
k<K ^ 

+ Pil^t"'^ - ^t"''^) = 3j G J^''{t),j' G J'^'{t),t< T,N^- V N^/ < n). 

l<ky^k'<K 

The first term is bounded by Ke'^'''^n^^ and the second term is at most 

n 

^ il + nNofP,,,^-,^„m,\ < 5]e|y^| 3t<T) 

l<k^k'<K-fl<i,i'<L ■m=l 

< (1 + '^^o)' [^Ce--''' + P.,,-.y, (|i?ftl < neR 3t < T) 

l<k=/=k'<K;0<i,i'<L 
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where we used (1.8) in the last Hne. By Lemma 2.6(c), the probabihty in the last 
term is at most ca.gA^^^^^^^'^"^^ (ni2)2(^+i)/3, and so if 5 = A~^^^^'^^'^~^\ 

P{V < oo) < cK^{L + ife^^^N^ \n^^ + n^e"''^ + d{nRf^'^+^y^' 

Now, optimizing over n and R, set = 23TE^ ~ '^<i^og{l/5) and n = [e'''^/^]. 
Here we may assume without loss of generality that > M{k) so that i? > 1. A 
bit of arithmetic now shows the the above bound becomes 

P{V < oo) < c(A',L)e'=''^<5=^/(2d+6)^ 

and the result follows. □ 

We suppose now that the assumptions of Theorem 1.2 are in force. That is, 
r > is fixed, .^q has law satisfying the local density condition (1.21) for a fixed 
r G (0, 1), and (1.23) holds. It is intuitively clear that the density hypothesis is 
weakened by reducing r. To prove this, note that the boundedness of the density 
and uniformity in x of the convergence in (1.21) shows that the contribution to the 
density on larger blocks from smaller blocks whose density is not near v is small in 
L^. We may therefore approximate the mass in a large block by the mass in smaller 
sub-blocks of density near v, and use the fact that the contributions close to the 
boundary of the large block is negligible to derive the density condition (1.21) for 
the larger blocks. As a result we may assume that r < 1/4. 

By inclusion-exclusion, it suffices to prove for —l<Lk<L, 



lim P(^^(rE^ + eyi^) = 1, j = 0, . . . , L^, A; = 1, . . . iT) 

K 



= n^l{^(%) = l' i = 0,...,Lfc})„(^,,.). (3.4) 
fc=l 

Allowing /c-dependence in and general subsets of the yiS is needed for the 
inclusion-exclusion, but to reduce eyestrain we will set ij = j and = L in what 
follows. The general case requires only notational changes. By the duality equation 
(2.17), 

PieAzik) = l,i = 0,...,L,k = l,...,K) (3.5) 
= P{CT{i,k) = l,i = 0,...,L,k = l,...,K), 

so (3.4) is then equivalent to 
P{CT{i,k) = l,i = 0,...,L,k = l,...,K) 

K 

-^llm{yi) = '^,i = 0,...,L})^^T,.>^^ ase^O. (3.6) 

k=l 

The proof of (3.6) uses the approach of [20], pp. 304-306. 

To work with the left-hand side of (3.6) we need the following preliminary result 
to simplify the initial inputs Co- Define 

13 = 1.9r and t, = . (3.7) 
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Lemma 3.2. Assume is independent of the reseated random walks {B^'"^ : w G 
eZ°'} as in (2.8). Then for any n G N and k > 0, 



0. 



lim sup k(ne§(^r'))-n^(^«) 

'^^"\wi\,...,\w„\<k,w,(^eZ''' ^i^i ^ 
w^j^Wj for i^^j 

Proof. For zi, . . . ,Zn G a^l/ define 

r(zi, . . . , z„) = {Bl^^' ezi + Qs for 1 < i < n} , 

and -f{wi, Zi) = P{Bl'^' G Zi + Qg), so that P{T{zi, . . . , z„)) = HILi li^i^Zi). Let 
G be the union of the events T[zi, . . . , Zn) over distinct zi, . . . , z„ G Oe^'^ such that 
\zi — Wi\ < ky/tl for 1 < i < n. We claim that P[G'^) is small for /c large enough. 
To see this, fix (5 > and choose k large enough so that 

P(|i?j"f ' - > {k - l)^e) = P{\Blf\ >{k- l)Vre) < 6/n . (3.8) 

By a standard estimate (and also since r < 1/4), for Wi as above and i ^ j, 

P(|i?,^f » - Bt:'"^\ < 2a,) < c|Q,|P(i?g = 0) < cmie-\r''/' < ce^'^'^'/' , 

which implies 

P{\Bl'^' - Bl'p\ < 2ae for some l<i<i<n)< Cn^e^'^^^'^'^l^ . (3.9) 
By (3.9) and (3.8), 

P{G') < Cn^e^'^- 1)3/4 ^ ^ _ (3 -^Q) 

Now consider the decomposition 

n n 

i = l Zl,...,Zn i = l 

where the sum is taken over only those (zi, . . . , z„) used in the definition of G. A 
typical term in this sum takes the form 

n 

E E(Jleo{z^ + e^)l{Bt^'"^ = z. + ei 

ei,...,enGQ£ i=l 



n n 

E E(jl^^,{z. + e.))llP{Bt^'"^=z, + eO . 



ei,...,e„eQ£ i=l i=l 

Since ^/t^ » a,, the probabilities P{B^^^'' = Zi + Cj) = P{B^f = Zi — Wi + 
Cj) are almost constant over G Qe- In fact, a calculation, using the version 
of the local central limit theorem in the Remark after P7.8 in [41] to expand 
P{Blf = Zi - Wi + Cj) = P{B^_2^^ = {zi - Wi + ei)/e), shows that 

P{Bl'^ = Zi — Wi + e) 
lim sup 7i = 1 (3.12) 

e,e'eg. P{Bl'^ = Zi-W, + e') 
\zi—Wi\<.k^/t^ 
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The continuous time setting is easily accommodated, for example by noting that 
along multiples of a fixed time it becomes a discrete time random walk. 

Consequently, for all sufficiently small e > 0, we have k^/tl < 1 and uniformly 

\n \ wi\ < k, \zi — Wi\ < ky/t^ and e £ Qg, 



1-6 < ' ^ <l + 6 . 

l{Wi,Zi) 

Using this bound and the fact that the Zi are distinct we have 



(3.13) 



n n 



Zi ~1~ €i 



i=l 



ei,...,e„6(3e 1=1 
n 



j{Wi,Zi) 



i=l 



l + 6rE(llD{z,,C^))p{r{zi,...,Zn)). 



i=l 



The continuity of v implies that for small enough e, for all \w\ < k and |z — < 
/c-v/tT, \v{w) — v{z)\ < 5. Also for sufficiently small e and z G ael'^, \z\ < k + 1, we 
have P{D{z,C^) > v{z) + 6) < 6/n. Thus 

n n n 

E(jlD{zi,eo)) <S + Yl{v{z^) + 6)<6 + ll{v{wi) + 25) . 



i=l 



1=1 



i=l 



Returning to the decomposition (3.11), the above bounds imply that for sufficiently 
small e, 

n n 

E(j[eoiBtn-^G)<{l + Sr[6 + ll{v{wi) + 26)] ^ W^i, •••,-«)) 

" ' 4 = 1 

n 

< {l + 6r\6 + ll{v{wi) + 26) 



i=l 



Zi,...,Zn 



i=l 



Let e — )• and then 6^0 above and in (3.10) to obtain 

n n 

hmsup sup (E(neo(i?iD) - n^(^^)) ^ 0- 

e->0 \wi\,...,\w„\<k^ Y=0 i=l 

A similar argument gives a reverse inequality needed to complete the proof. □ 
We break the proof of (3.6) into three main steps. Introduce 



S = T-te = T-eP. 
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Step 1. Reduction to Bernoulli inputs and K = 1. 

Let X = X^''^ be the modification of the dual in which particles ignore reaction 
and coalescing events on [5", T], and let Ct be the associated computation process 
with inputs Co(j) = ^oi^r)- That is, Xf = Xf for t G [0,5], and during the time 
period [S, T], Xj , j £ JiS) follows the same path as X/ until the first time a reac- 
tion or coalescence occurs, at which time all the X^ switch to following completely 
independent random walks. 

On the event defined in (2.53) there are no reaction or coalescing events 
during [5*, T]. Thus, Xt = Xt for all t G [0,T] on Gj,, so it follows from Lemma 2.9 
that 

P{Q / for some t G [0,r]) < C2.9[(^ + l)Kfe'"'^ei>^3 . (3.14) 
Let fpeix) = Pf^^Q^x), where 

P^^f{x) = E{f{x + Bf)), X £ eZ'^, is the semigroup of B' , (3.15) 

and let Wi, W2, ... be an iid sequence of uniforms on the interval [0, 1], indepen- 
dent of and the random variables used in Section 2. We will use this sequence 
throughout the rest of this section and also in Section 4. Define a second computa- 
tion process Ct I'^e ^ i ^ T, for X, with inputs 

Cl{j) = HWj<v{X^^)}, j€J{S). (3.16) 

It is clear that conditional on 0"(^q)V7oo5 the variables Ct* (j)i J £ '^(5'), respectively 
Ctsij)-, 3 £ are independent Bernoulli with means u(X^), respectively tp^^Xg). 

Let X = X^''^ be the branching random walk dominating X which was introduced 
in Section 2.8. If we fix 5 > 0, then using Lemma 2.4(a) it is not hard to see that 
there exist n, k such that for all e sufficiently small, 

P{\J{S)\ < n and |X|(5)| < k for all j G J{S)) > I - 5. (3.17) 

It now follows from (3.17), Lemma 3.2 and the definitions of X, Co and Q that for 
any 6 : Z+ ^ {0, 1}, 

\P{CtAj) = bjj e J{S)) - PiClU) = bj,j G J{S))\ 

< ^(|P(4(i) = bj,j G JiS)\Ts V a{eo)) - PiClU) = bj,j G AS))\Ts V ai^oW 
— ^ as e — ^ 0. 

As a consequence, since both (tyQ^'^e ^ t < T are defined relative to X with 
identical {Um}, {/^m} and {Rm}, by conditioning on the input values, the above 
implies 

P{CT{i,k) = l,i = 0,...,L,k = l,...,K) 

- PiCUi, k) = l,i = 0, . . . , L,k = 1, . . . , K) ^ as e ^ 0. (3.18) 

Let Cj*'^'' be the computation process associated with X'^*"'^ , 1 < k < K with 
inputs as in (3.16). That is, for j G J^''{S) there exists a j' G J{S) with = X^'"-' 
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(by (3.2)) and we set Ctf'ij) = H^f < viXg"'^)}. Up to time V = V^^T,e the duals 
X^kyT ^ k < K, use independent random walk steps and branching mechanisms, and 
on {V = 00} the computation processes also use independent uniforms and 
parent variables as well as independent inputs at time t^. It follows that (see below) 

\PiC^{i,k) = l,i = 0,...,L,k = l,...,K) 

K 

-llP{C/'{i) = 1,0 <i<L)\ 
k=i 

< P{V < 00) ^ as e ^ 0. (3.19) 

The last limit follows from Lemma 3.1 and (3.1). Perhaps the easiest way to see 
the first inequality is to extend X^''- to t G [V, T] by using independent graphical 
representations and define the corresponding computation processes C using in- 
dependent collections of {Vl^j}'s for the inputs at time t^. The resulting computation 
processes are then independent, each *'^'= is equal in law to and the 

two are identical for all k on {V = 00}. On this set we also have 

{C^{i,k):i,k} = {CT"={i):i,k}, 

and so (3.19) follows. It is therefore enough to set K = 1 and drop the superscript 
k. Altering our notation to z = (zi), Zi = Xs + syi where — )■ x, it suffices now to 
prove 

PiCHi) = l,0< i<L)^ (l{e(y.) = 1, < i < as e ^ 0. (3.20) 

Step 2. Reduction to L = 0. Let X = X^'^,0 < t < T he the branching random 
walk started at z, with associated computation process Ct,te <i I^T. We suppose 
that X and X are coupled as in Subsection 2.7, and that Qt has initial inputs 

QtAj) = m3<v{xi)},j(^j{s). 

On the event G^, J{S) = J{S) and all the differences j e J{S) are small. 

It therefore follows from (3.17) (if we take y = X will stochastically dominate 
X), the continuity of f , the definitions of and Qt^-, and Lemma 2.9 that 

P(G?, CI (j) = 4 (i) for all j G J{S)) ^ 1 as e ^ 0. (3.21) 

By Lemma 2.10, on the event in (3.21), the outputs Ct and Ct agree, and conse- 
quently 

P((* (i) = 1, < i < L) - PiCrii) = 1, < i < L) ^ as e ^ 0. (3.22) 

Using the branching structure we can now reduce to the case L = 0. To see this, 
let X^'-''^ be the branching random walk started from Zi = + eyi, with associated 
computation process Ct'^'^e l^t <T with initial inputs Ct'iJ) which, conditional on 
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X^\0 < t < S are independent with means v{X^^'-'). The branching property and 
definition of Xq in (2.23) imply (recah z/g from just above (2.21)) 

P(CT(i) = l,i = 0,...,L)= 5^ i/.(7r) H P(C^^(0) = 1). 

Since — t- x as e — >• for z = 0, . . . , L, if we can estabhsh 

P{C,^- (0) = 1) ^ n(r, x) as e ^ 0, (3.23) 
for some u : x R'^ — )• [0, 1], then the convergence =^ ^0 iniphes 

P(CT(i) = l,i = 0,...,L)^ ^ z.o(vr)(n(r,x))W 

vren^ 

= (l{^(yO = 1, < i < L}),,(r,,) as e ^ , (3.24) 

where (1.26) is used in the last line. Combining this with (3.22) gives the desired 
result (3.20) but with u in place of u, that is, we get 

P(Ct(^) = 1, < i < L) ^ {l{i{yi) = 1, < i < as e ^ 0. (3.25) 

We first turn now to the proof of (3.23). 

Step 3. Convergence and identification of the limit. Let X^ be the branching Brow- 
nian motion started at x G M'^ run over the time period [0, T], with associated 
computation process Q ,t£ < t < T with inputs 

Cl{j) = l{W,<v{x'/)},jeJ'{S). 

Using the obvious analogue of (3.17) for X^, the continuity of v and the definitions 
of Ct^ and Q^, Lemma 2.13 (and the uniform convergence of gf to gi) implies 

PiG"/ , 4 (j) = C° (i) for ah j G jO (5) ) ^ 1 as e ^ 0. (3.26) 

By Lemma 2.13(c), on the event in (3.26), Ct(0) = Ct(0)' thus 

P{C^'{0) = 1) - P(C°(0) = 1) ^ as e ^ 0, (3.27) 

where we note that both quantities in the above depend on e. If we take the initial 
inputs for the computation process = C^'* at time to be 

C°'*(j) = 1{W, < vix'^'^)},j G J0(r), (3.28) 
it is now routine to see that P(C£ = Ct^*) — ^ 1 as e — )• 0, and so by (3.27) 

limP(C^^(0) = 1) = P(C^'*(0) = 1) ^ u{T,x). (3.29) 

£ — 

This proves (3.23), hence (3.24) and so to complete the proof of (3.20), and hence 
Theorem 1.2, we only need show u = u: 
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Lemma 3.3. Let Xjr,0 < t < T be the branching Brownian motion started at 
X £ W^, with associated computation process Ct'* ^t <T with initial inputs as in 
(3.28). Then 

P(C°'*(0) = l) = n(r,x), 
where u is the solution of the PDE (1.19). 

Proof. This is very similar to the proof in Section 2(e) of [20]. Recall Pf is the 
semigroup of B'^ . Let x & and G el/ satisfy |x — x^l < e and let be om' 
rescaled particle system where {^^{sy) ■ y G Z'^} are independent Bernoulli random 
variables with means {v{ey) : y € Z*^}. If 

de{ey,Ce) = -ayWo{y,0 + a-aym{y,0, eG{0,l}^', (3.30) 

then the martingale problem for shows that (cf. (2.25) of [20]) 

Eierixs)) = E{P^eoixe)) + r x E{d,{B'T^^,C)) ds, (3.31) 

Jo 

where Bq = x^ under P^^ ■ Our hypotheses on imply 

E{P^eo{xe)) = PtA^s) ^ Ptv{x) ase^O, 

where Pt is the d-dimensional Brownian semigroup with variance o"^ . Recall we have 
proved ((3.24) and the preceding results) that 

lim P{^^{xe + eyi) = r]i, i = 0, . . . , L) = {l{^{yi) = r]i, i = 0,... ,L})^it,x)- 

Now use the above with Fubini's theorem, the uniform convergence of gl in (1.9) 
and the coupling of and B in Lemma 2.11 to conclude that 

\iTnE,^xE{de{B'T-s.es)) 

= lim^,, xE^Ey (-es{B'T-s)9l{es{B'T-s + eY^), UB't-s + ^^^°)) 

+ (1 - esiB^T-^sWMiB^T^s + eY'), esiB^T_, + ey^o)) 

= i?,((-e(0)/lo(0,0 + (1 - m)hl{0,O)u{s,Br^,)) 
= E,{f{u{s,BT-s)), 

the last by (1.16). Now use the above to take limits in (3.31) to show that u solves 
the weak form of (1.19). As in Lemma 2.21 of [-' '] it follows that u solves (1.19) 
and so equals u. □ 

The following asymptotic independence result follows easily from Step 1 in the 
above argument. 
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Proposition 3.4. If K ^ N, there is a C3,^{K) so that if zi, . . . , zk S e^"^ satisfy 
infjjik \ zj — Zk\ > £^^^ and is deterministic, then 

K K 
\E{J{ iWk)) - n Em^k))\ < cs.me^^'^e^". 

k=l k=l 

Proof. Define V = Vz^T,e as in (3.3) but now with z^ G eZ*^, tliat is L = 0. Use tiie 
dual equation (7.11) and argue just as in the derivation of (3.19) to see that 

K K 

1^(11 ^t(^.)) -ll^mz^))\ < P{V < CX)). 
i=l i=l 

The fact that is deterministic makes the independence argument simpler in this 
setting. Now use Lemma 3.1 and the separation hypothesis on the z^s to bound 
the right-hand side of the above by 

C3.i(K,0)e'=''^e(3/4)(d-2)/{d+3) < C3.i(i^, 0)6^" V/^. 

□ 

3.2 Proof of Theorem 1.3 

Proof. Let t > and choose ri{e) | so that ?/(e)/e — )• oo and ri{e)/6{e) — )• 0. 
Recall Isix) is the semi-open cube containing x defined prior to Theorem 1.3. Write 

E{{u^{t,x) -u{t,x))'^) = 

= {^T ^ E{edxi)etix2)-uit,xmix,)+et{x2))+uit,xf). 

The contribution to the above sum from — < ^/(e) is trivially asymptotically 
small, uniformly in as e — >• 0. Theorem 1.2 shows that the expectation in the 
above sum goes to zero uniformly in xi,X2 G hix), \xi — X2I > ??(e), x G [—K,K]'^ 
as e — 7- 0. The result follows. □ 
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4 Achieving low density 

The first step in the proof of Theorem 1.16 is to use the convergence to the partial 
differential equation in Theorem 1.2, and more particularly the estimates in the 
proof, to get the particle density in (1.20) low on a linearly growing region. 

As we will now use the partial differential equation results in Section 1.2, we 
begin by giving the short promised proofs of Propositions 1.4 , 1.5 and 1.7. 

Proof of Proposition 1.7. Set r] = \r\/3 and let Ljj, Cq and cq be the constants in 
Proposition 1.6, and define 

Ls = Ll ci = co, ts = ■ 3Vd/\rl Ci = (Co V l)e'=«*^ 

Suppose t > ts, L > Ls and |x| < L + {\r\/3)t/^/d. Then we may write x = xq + y, 
where 

|y|<M* and \xo\<L-^<L-^-^ = L-Ll (4.1) 

For t > and z € define u{t,z) = u{t,xo + z). If \z\ < L^, then \xo + z\ < 
\xo\ + < L, which implies that u{0,z) < p — 6. Applying Proposition 1.6 to n, 
and recalling the bound on \y\ in (4.1), which implies \y\2 < we have that for 
t > ts, and |x| < L + {\r\/3)t/\/d 

u{t,x) = u{t,y) < Coe-'="* < Cie'^'K 

Since the right-hand side above is at least 1 if i < t^, the above bound follows for 
all t > 0, and we have proved the result with w = \r\/6^/d. □ 

Proof of Proposition 1.5. Extend /|[o,i] to a smooth function / on [0, 1 + 5o] so that 
/ > on (1,1 + So), /(I + 6o) = 0, /'(I + (5o) < and f{u)du < 0. The 

situation is now as in Proposition 1.7 with 0, 1 and 1 + Sq playing the roles of 0, p 
and 1. As the solutions take values in [0, 1] the extension will not affect the solutions 
and the Theorem follows from Proposition 1.7. □ 

Proof of Proposition 1.4- The version of Proposition 1.5 with the roles of and 1 
reversed, applied on the interval (0, a) shows there are positive constants L, c, and 
C so that if u{0,x) > a/2 for |x| < L, then 

u{t, x) > a - Ce'"^ for |x| < L + 2wt. 

It is here that we need /'(a) < 0, corresponding to /'(O) < in Proposition 1.5. By 
Theorem 3.1 of Aronson and Weinberger [2] there is a Tq so that 

u{T,x) > a/2 for |x| < L and T > Tq. 

Therefore we have 

u{t + To,x)>a- Ce'"^ for |x| < L + 2w{t + Tq) - 2wTo, 
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and so for t > 2Tq, 

u{t, x)>a- Ce'^^oe""* for \x\ < L + wt. 

The result follows as we may replace w with 2w. □ 

Recall the parameter r € (0,1), and definitions of a^, t^, Qe-, and D{x^£) in 
(1.20). We first show the density D{x,^j') is close to its mean. 



Lemma 4.1. Let T > and assume is deterministic. 

(a) If < r < then for all x G a^'Z'^, 

E{{D{x, er) - E{D{x, e^)))^) < C^je'^^^e'/^ 

(b) IfO<r < l/(16d) andC = y + [-L, L]"' for y G M"', then for all i] > 0, 

P( sup \D{x,^'t) - E{D{x,(l^))\ >v)<Cj ^e^'^^L'^e^^'^r]-'^. 

Proof, (a) Note that 

\{{zi,Z2)e{x + Q,f: \zi - Z2\ < e'/^}\ (4.2) 
< (2e-3/4 + < c,|Q,p(£^^)'^. 

If denotes the sum over 

z G {(zi, Z2) G (x + Q.)' : \zi - Z2\ > e'/^}, (4.3) 

then by (4.2) and Proposition 3.4 with K = 2, 

E{{D{x,eT) - E{D{x,i^))f) 

2 2 
k=l k=l 

<c,e(i-0'^ + 4c3.4(2)e^''V/8 

< C^ ie^^^e^l^, (4.4) 

where our condition on r is used in the last line, 
(b) Note that 

\C n OeZ'^l < CdL'^a-^ < CdL'^e-"^ < CdL'^e-^'^^ 
The result now follows from (a) and Chebychev's inequality. □ 

We recall the following hypothesis from Section 1.6: 

Assumption 2. There are constants < ui < 1, £2,6*2, w > 0, Lq > 3 so that for 
all L > Lq, if u{0, x) < ui for \x\ < L then for all i > 

u{t,x) < C2e~'^2t fQj. all \x\<L + 2wt. 
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We also recall the following condition from the same Section: For some ro > 0, 

1 

X^lbf -5i||oo < ci.ese^". (4.5) 

i=0 

We say that ^ € {0,1}'^^'* has density at most k (respectively, in [^1,^2]) on 
^ C M'^ iff D{x,C) < K (respectively D{x,C) G [ki,K2]) for all x G (a^Z'^) n A. We 
set (recall (3.7)) 

r = ^, hence /3 = ^,t, = ei-9/(i6'i),r = ^4.2 log(l/e), and 5 = T - t„ (4.6) 
if 1 



where ^4.2 = c, (job ^ f 

Lemma 4.2. Suppose Assumption 2 and (4.5) /lo/c/. Lei U2 G (0, ui) and 

'^'^.^ ~ ^ -*■) (i^w ^ ■5^)- There is an e^ 2 > 0; depending on (ui, U2, C2, C2) 
and satisfying 

e// < ^^2, (4.7) 

so that if0<e<e^2 ^''^^ 2 + Lq < L < £-•001/'^^ then whenever has density at 
most U2 in [— L, , 

P(Cf. has density at most e'^-^ in [-L -wT,L + wTf\^^) > 1 - e-°^ 

Note that (4.7) allows us to iterate this result and obtain the conclusion on succes- 
sively larger spatial regions at multiples of T. 

The proof of the above Lemma will require some preliminary lemmas. 

Lemma 4.3. Ifpfiy) = e~'^P{Bf = y), y e eTL'^ , then for < e < 1, 

\p\ix) -p\(x^y)\< c^.5|y|t"^''+^^/^ for all x,y £ eZ'^ and t > 0. 

Proof. This is a standard local central limit theorem; for d = 2 this is Lemma 2.1 
of [5] and the same proof applies in higher dimensions. □ 

Recall (from (3.15)) that is the semigroup associated with B^. 

Lemma 4.4. There is a ^ such that if 1 > a > (3/2, then for < e < 1, 

for all x,x' G eZ"' such that \x - x'\ < 2e" and all ^ G {0, ly^'' . 
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Proof. Let — oo <6<a, A = x — x' and assume |A| < 2e". Apply Lemma 4.3 to 
see that 

\PtW-P!M^')\ 

< \P{Bl = z) - P{Bl = z + A)\ 

< Y C4.3e'|A|t,-^'+'^/' + ^(l^t I > 3e') + P{\Bl\ > 3e' - A) 

If we use Chebychev to bound the last summand by ctge"^"^ = ce^""^^ and optimize 
over 6 (setting = f — °^~2+d'^ < ^ < a), we obtain the required upper bound. □ 

Lemma 4.5. For any i] > there is an e^.s^rj) > so that if < e < e^,5, 
u G [0, 1], L > 1, and G {0, 1}^^'' has density at most u in [—L,L)'^, then 

PfMx) <u + r] for allxe [-L + 1, L - 1]'' n eZ'^ . (4.8) 

Proof. By translation invariance it suffices to prove that for small enough e > and 
all x G [— Og, ae]*^ H eZ'^, if has density at most u in [— 1, 1)"^ then P^^^{x) < u + r]. 
(This addresses the uniformity in L.) Argue as in the upper bound in (3.13) to see 
that for e < £o{i]), 

'^'Jrif*' ^^J^^^ < 1 + ^ for all z G a,Z^ |z - x| < 1 and e G Q,. 

We therefore have 

piax) < p{\B:f\ > 1/2) 

+ Yl n\z-x\<3/4)Yaz + e)^^PiBtfGz + Q,) 

< Aa'^dte + Y ^(1-^ ~ ^1 - 3/4)u(l + ^)P{Blf ez + Q,) 

< Aa^dts + u + ^ <u + r], 

for e < £i{r]). □ 

We are ready for the Proof of Lemma 4.2. 

Proof. By the conditioning we may fix a deterministic .^g as in the statement of 
Lemma 4.2. In light of Lemma 4.1 our first and main goal is to bound E{D[x,^'^)) 
for a fixed x G ffleZ"^ U [-L - wT,L + wT]'^. Let z G x + Qe- Let X^'"^ be the 
modification of the dual X^''^ , starting with a single particle at z, in which particles 
ignore branching and coalescing events on [S, T] by following their own random walk 
and switching to independent random walk mechanisms when a collision between 
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particles occurs. Hence X'g = X'g for all i G J{S)^ and on [S,T] the particles 
in X^'-^ follow independent copies of . Let C,'^ be the associated computation 
process, defined just as ("^ is for X^''^, with initial values Co(i) = Co(^t)'-? ^ '^('S') 
(for X^'"^ the index set is constant on [S,T]). On G^, T ^ U^!:~^[i?m,, -Rm + 2e^], 
with /3 < 1/2, and so 

[S,T] n(U^Io^ + V^) = 0. (4.9) 

Therefore on G^, X^''^ has no branching or coalescing events on [5", T], and so 
X^'"^ = X^'"^ on [0,T]. This also means (by (2.17)) that, given the common inputs 
cm = cm = ^oi^i), j G J(T) = J{S) we have 

cm = cm = em on g^. (4.10) 

Let V'e(^) = PtAo^^)- Conditional on Fs, {X;^ — X^ : j G >/('S')} are iid with 
law Pq{B^^ G •), and so, conditional on {Ct-s^^) ~ ?o("^t) • 3 ^ ■^('S')} are 
independent Bernoulli rv's with means {ipsiXg) : j G J{S)}. Recall {VFj} is an 
iid sequence of uniform [0, 1] rv's independent of J-'oo (that is of our graphical con- 
struction). Let {Ct'*U) '■ j S "^(^ — t),T — S < t < T} he the computation process 
associated with X^'"^ but with initial inputs Ct-s(^) ~ ^O^j — i^ei-^s) '■ 3 ^ ^i^)}- 
Then {Q : T — S < t < T} and {Ct'* ■ T — S < t < T} have the same law 
because the joint law of their Bernoulli inputs and the processes X^''^,t < S and 
{{fj-m, Um)^{Rm 1^ t),t < S) uscd to define them are the same. Therefore by (4.10) 

iP(ei^(z) = i)-p(cr(o) = i)i 
= iPiem = 1) - picm = i)i < pug^tT)- (^-n) 

Consider now the branching random walk X starting with a single particle at z 
and coupled with X^'-^ as in Section 2.7, together with its computation process {("^ : 
t G [T-S,T]} with initial inputs Cr-s(j) = UWj < '0e(X^)),j G J(5). Conditional 
on Too, these inputs are independent Bernoulli rv's with means {ip{X-^g) : j G j(S')}. 
The computation processes C^'* and are identical on [T — S, T] if given the same 
inputs at time T — S. Therefore Lemma 2.10 shows that on G^ Ct(0) and C^*{0) 
will coincide if given the same inputs at time T — S. Therefore 

|p(CT(o) = i)-P(cr(o) = i)l 

< P((G^)^) + E{P{CmU) + Cr-sU) 3i G J(5)| J-oo)l(G^)) 
<P((G^)^) + i?( |V'.(X^)-V.(^^)|1( sup |X^-X^|<.V6)y (4.12) 
jeJ{s) ■'■^^(■^) 

Use Lemma 2.9 to bound the first term above and Lemma 4.4 with a = 1/6 > /3/2 
to bound the second, and combine this with (4.11) to conclude that (use d> 3) 

mm = 1) - PiCm = l)l < 2c2.9e=''^e^/3 + E{\JiS)\)c^,^e^'/'-^^/^''-^''^) 

< 2c2.9e^/(^°'^) + e^''^C4.4e(^/3-/3)/(2+d) 

< (2^2.9 + ^4.4)^'/^'°'^- (4-13) 
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To prepare with the couphng with the branching Browian motion we must extend 
i/^eix) = Pf ^q{x) from eZ'^ to in an appropriate manner. Since has density at 
most U2 in [—L,L]'^, Lemma 4.5 shows that for e < £4 5(1*1 — U2) we may extend 'i/'e 
in a piecewise hnear manner so that 

< uil{\x\ <L-2) + l(|x| > L-2) for all x G M'^. (4.14) 

In addition, using Lemma 4.4 with a = 1/6, we may assume the above extension 
also satisfies 

IV'e(x) - ipe{x)\ < 04^46 ' for x,x' £ R'^ such that \x - x'\ < e^^^. (4.15) 

Now consider the branching Brownian motion starting with a single particle 
at z, coupled with as in Section 3. Consider also its associated computation 
process on [T — S, T] starting with conditionally independent Bernoulli inputs 
{Ct-s'(j) = l(^j — ^^(^s"')) • J S J^i^)}- We may argue as in the derivation of 
(4.13), but now using Lemma 2.13, (4.5), and (4.15) in place of Lemmas 2.9 and 
4.4, to conclude after some arithmetic using d > 3, 

\P{Cm = 1) - PiC'AO) = 1)1 < e^'^C2,u{e'/' + ci.gge'^''/^) + 04.46^(40^). (4.I6) 

By Lemma 3.3 (we have shifted time by T — S), P{Cj^{0) = 1) = U£{S,z), where Ue 
is the solution of the PDE (1.19) with initial condition ^^(O, •) = V'e- Now combine 
this with (4.13) and (4.16) to see that for small enough e as above 

Eieri^)) < ce(V(40d))A(ro/4) ^ ^^(^^ 

Now use the bound on the initial condition (4.14) and Assumption 2 in the above 
to conclude that for \z\ < L — 2 + 2wS, and small e 

where t' = (^f^ (ilM ^) > 74.2 used the definition of S and some 

arithmetic. By taking e smaller if necessary we may assume 2wS — 3 > wT and so 
the above holds for \z\ < L + 1 + wT. This shows that 

E{D{x, Ct)) < eT' for x G a^Z'^ n [-L -wT,L + wTf. (4.19) 
Finally apply the above and Lemma 4.1 to conclude that for small enough e 
P( sup D(x,C^) >e^4.2) 

<P{ sup |£'(x,Ct) -^(Z)(x,^^))| > £^4.2/2) 

xe[-L-wT,L+wT]'ina.Z'i 

< C4jei/i6(L + it;r)'^e^''^4£-'^4.2 

^ (^^l/16^-.001g-l/100d^-l/120d ^ ^.05 

□ 
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Our next goal is to show that the dual process only expands linearly in time. 
The first ingredient is a large deviations result. Recall the dominating branching 
random walk {X^'^{t) : j G J'ii)} introduced at the beginning of Section 2.8 which 
satisfies 

{XI'^:jeJ%s)}c{XI'^:jer{s)}. 

If 1 1 X| I loo = sup{|X^'-^(s)| : j G J^{s)} and similarly for ||X|||oo, then the above 
domination implies 

ll^slloo > ll^f lloo for all s,e. (4.20) 

Recall c* is as in (2.7). 

Lemma 4.6. Assume starts from one particle at 0. For each R > there is an 
^4 di^*^-^) ^' nonincreasing in each variable, so that for < £ < £^ (j and t > 0, 

P{\\Xt\\oo > 2pt for some s < t) < {4d + 1) e^p{-t{j{p) - Cb)) for all < p < R, 

where 7(p) = min{p/2, p^/So"^}. Moreover, if p > max{4cfe, 2(T^}, then the above 
bound is at most (4d + 1) exp(— 1/9/4). 

Proof. The last assertion is trivial. Let be a random walk that starts at 0, jumps 
according to at rate and according to at rate c*. Since ^'l-^f | = exp(cbt) 
by summing over the branches of the tree, it suffices to show 

^'(ll'S'f lloo > pt for some s < t) < (4d + 1) exp {-j{p)t) . (4.21) 

As usual, Bf is the random walk that jumps according to Ps at rate By the 
reflection principle 

P ^supSf'^ > pt^ < 2P [pf > pt^ < 2e-'^P*E (^exp (ePfYj 

for any > 0. If </>e(^) = X^^, e^^^Pe{x) then a standard Poisson calculation gives 

E (exp (ePt^')) = exp{te-\MO) - !))• 
By scaling (pei^O) = (j)i{£9). Our assumptions imply (j)'i{0) = and </>'/(0) = cj^ so 

£-2(01 (e^) - 1) ^ ase^O. 
li < p < R and 6 = p/cr"^ in the above, it follows that for e < eQ{R), 
e-ipV<^'^)tE (exp ({p/a)Bf)) < eM-p't/^a^), 

and so, 

P(sup \Bf\ > pt) < 4 exp {-p'^t/3a^) . (4.22) 

s<t 

Let Jf be the one dimensional random walk that jumps according to the law of 
eY* = emaxj<7Vo \Y^\ at rate c* , and notice that this will bound the L°° norm of 
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the sum of the absolute values of the jumps according to Qs in up to time t. If 
we let (pj{0) = E{exp{6Y*)), then, arguing as above, we obtain 

P{J! > pt) < exp(-p^t + c*t{^j{e6) - 1)). 

The exponential tail of Y* (from (1.8)) shows that {(pj{£9) — l)/e6 — ^ EY* as e — ^ 0, 
and so, if we set = 1, then for small e, c*((/)j(e) — 1) < p/2. (The choice of e here 
works for all p because we may assume without loss of generality that p > po > 
as the Lemma is trivial for small p.) Therefore 

P{J! > pt) < exp{-pt/2). (4.23) 

To derive (4.21), write 

d 

P(sup lloo > 2pt) < P(sup \Bf\ > pt) + P( > pt), 

S<t s<t 

and use (4.22) and (4.23). □ 

Our next result uses the large deviation bound in Lemma 4.6 to control the 
movement of all the duals that start in a region. Recall that X^'^ is the dual 
for starting with one particle at x from time U. For x E M"^ and r > let 
Q{x, r) = [x — r, X + r]'^ and Q^{x, r) = Q{x, r) n eZ'^. Write Q{r) for Q{0, r) and 
Q'{r) for Q"(0,r). 



Lemma 4.7. Forc>0,b> Aca V la"^ , L > 1 andU >T' = clog(l/e), let 

Pe{b, c, L, U) = P{X^'^ is not contained in Q{L + 2bT') 

for some u £ [U — T', U],t < T' and some x G Q^{L)). 

Let = 12(4(i + 1)3^^. There exists e^ g{c*,b) > such that ifO<e<£^g 

Pe{b,c,L,U) < c'aL^{clog{l/e) + l)e'?-'^ 
where g = (^ - 2) Ae'^. 

Proof. By translation invariance it suffices to take U = T'. For x G eZ'^ let {Tj(x) : 
i > 0} be the successive jump times of the reversed Poisson process, starting at time 
T' , determined by the T^,Tn^ . Also let Nx be the number of such jumps up to time 
T' , so that Nx is Poisson with mean (c* +e~^)T'. The process ^f(x) is constant for 
t G (T' - Ti+i{x),T' - Ti{x)] and for such t the dual X'^'^v) is 

-^x,T-T,{x)^^ + (T' - Ti{x) -t)), 

that is, one is a simple translation of the other. This means for t as above 

U,<tX-'\v) C U,<r,_T,(.)X-'^'-^^(")(^;), (4.24) 
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(in fact equality clearly holds). As a result, in c, L, T') we only need consider t 
to be one of the times T' — Ti{x) for < i < and we may bound 1 —ps{b, c, L, T') 
by 

P(3x e Q^L) s.t. 7V^. > 3T'(e-2 + c*)) 

+ P(3x £Q%L), < < ST' (£-2 + c*) s.t. 

sup ||X^'^'-^'(^)(z;)||oo >26r') 

t;<r'-T,(a;) 

< {2Le-^ + l)'^exp{-3r'(e-2 + c*)}S(e^-) 

+ (2Le-^ + l)'^(3r'(e-2 + c*) + l)(4(i + 1) exp(-r'6/4). 

Here we are using Lemma 4.6 and the strong Markov property at Ti{x) for the 
filtration generated by the reversed Poisson processes J-t- Some arithmetic shows 
the above is at most 

3'^(L V e) [exp(-3T'(e"2 + c*)) exp((£^2 ^ _ 

+ (4d + l)(3T'(e^2 + c*) + i)eW4 

< S'^iL V e) [exp(-r'(e"2 + c*)) + 6(4(i + l)(clog(l/e)e~2 + 

< 3'^(L V Te^^^"") + 6(4(i + l)(clog(l/e) + l)e^^/^-^ 



<c'rf(LVe)'^(clog(l/e) + l)e-'^e 



(bc/4-2)A£' 



□ 
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5 Percolation results 



To prove Theorems 1.15 and (especially) 1.16 we will use block arguments that 
involve comparison with oriented percolation. Let D = d + 1, where for now we 
allow d > \, and let A be any D x D matrix so that (i) if x has xi + ■ ■ ■ + x^ = 
1 then {Ax)d = 1, and (ii) if x and y are orthogonal then so are Ax and Ay. 
Geometrically, we first rotate space to take {1/D, . . .1/D) to (0, . . . , 0, 1/y/D) and 
then scale x — )• x\fD. Let Co = {Ax : x G Z^}. The reason for this choice of lattice 
is that if we let Q = {Ax : x G [-1/2, 1/2]^}, then the collection {z + Q, z £ Co} 
is a tiling of space by rotated cubes. When d = 1, £2 = {('^j n) : m + n is even} is 
the usual lattice for block constructions (see Chapter 4 of [M]). 

Let T-Lj^ = {z G Cd '■ zd = k} = {Ax : x G = k} be the points on 

"level" k. We will often write the elements of T-Lk in the form {z, k) where z G W^. 
Let l-L'f^ = {z G M'^ : {z,k) G Hk]- When d = 2, the points in "H'o are the vertices 
of a triangulation of the plane using equilateral triangles, and the points in Ti'i are 
obtained by translation. One choice of A leads to Figure 7, where T-L\ and ^2 ^^'^ 
obtained by translating 'H'q upward by \/2 and 2\/2, respectively, and ^3 = T-L'^. 




Figure 7: Ti'o (black dots) and Ti'i (white dots) in £3 

In d > 3 dimensions (the case we will need for our applications in this work) 
the lattice is hard to visualize so we will rely on arithmetic. Let {ei, . . . , en} be the 
standard basis in R-^, and put Vi = Aet, i = 1, . . . , D. By the geometric description 
of A given above, Vi G 7ii has length VD, and writing Vi = {v'^,1), v'^ G has 
length yjD — \. For i j, \Wi — v'j\\2 = \\vi — Vj\\2 = V2-D, the last by orthogonality 
of Vi and Vj. The definitions easily imply that T-L'^j^i = v[ + H'^ = {v[ + x : x £ "H'^} 
for each i and k. Note that Dv[ G because Dvi — (0, . . . , 0, D) G T-Lq. This implies 

For X G T-L'i^ let Vx C be the Voronoi region for x associated with the points in 
i.e., the closed set of points in that are closer to x in Euclidean norm than to 
all the other points of H'^^ (including ties). If V = Vq (in d = 2, Vq is the hexagon in 
Figure 8 inside the connected six white dots), then the translation invariance of 
and fact that T-L'^ = kv[ + H'q show that Vx = x + V for all x G . It is immediate 
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Figure 8: T-L'q (black dots) and Voronoi region about (inside white dots) in 



from the definition of Voronoi region that for each A:, 

Furthermore, Vx is contained in the closed ball of radius D centered at x. (To see 
this we may set x = A; = and transfer the problem to Z"^ via . It then amounts 
to noting that if x G satisfies ^ Xj = and ||x||2 > V^, then there are i ^ j 
s.t. Xj > l,Xj < or Xj < — l,Xj > 0, and so ||x it (cj — ej)\\2 < ||2;||2-) From this 
inclusion we see that for any L > 0, 

if CL =L/{2D) then clV^ C c^x + Vf, (5.2) 
and so U^-g-^/^ c^x + [-L, Vf = M'^. 

The above also holds with 2cl in place of cl but the above ensures a certain overlap 
in the union which makes it more robust. Finally, one can check that for some 
positive C5.3(-D), 

if a G Vxi b ^ Vx and |a — 5| < C5.3 then b G Uj-^jV^.+,;'_^' . (5-3) 

For this, note that x + v'^ — v'j, 1 < i ^ j < D are the D{D — 1) "neighboring points 
to x" in Ti'f^, corresponding to the 6 black vertices of the hexagonal around in 
Figure 8 for x = and D = 3. The above states that the D(D — 1) corresponding 
Voronoi regions provide a solid annulus about Vx, as is obvious from Figure 8 for 
D = 3. 

Our oriented percolation process will be constructed from a family of random 
variables {r]{z),z € Cd} taking values or 1, where means closed and 1 means 
open. In the block construction, one usually assumes that the collection of r]{z) is 
"M dependent with density at least 1 — 0" which means that for any k, 

Piviz^) = l\r]iz,),j^i)>{l-e), (5.4) 
whenever Zi G Co, 1 < i < k satisfy \zi — Zj\ > M for all i 7^ j. 
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Our process will satisfy the modified condition 

P{rj{zh) = l\ri{zj),j < k) > {1 — 0) whenever Zj = {Zj,nj) G Cd, I < j < k (5.5) 
satisfy Uj < rik or (rij = and \z'j — z'^\ > M) for all j < k. 

It is typically not difficult to prove results for M-dependent percolation processes 
with 6 small (see Chapter 4 of [14]), but in Section 7 we will simplify things by 
applying Theorem 1.3 of [.34] to reduce to the case of independent percolation. By 
that result, under (5.4), there is a constant A depending on D and M such that if 

'-^'K'- (a-T)^Ua )('-W^-'»'^^)- 

we may couple {r]{z), z G Co} with a family {Ciz), z £ Co} iid Bernoulli random 
variables with P{(^{z) = 1) = 1 — 6' such that C(z) < r]{z) for all z £ Co- An 
examination of the proofs of Proposition 1.2 and Theorem 1.3 of [34], shows that the 
above result remains valid under our condition (5.5). [In their proof of Theorem 1.3 
we can order the vertices of a finite set in Cd so that the levels n of the vertices are 
non-decreasing, and then in the inductive proof of Proposition 1.2 we will only be 
conditioning {r?(zo) = 1} on vertices whose level is at most that of zq.] 

In view of the comparison, and the fact that 0' — )• as — >• 0, we can for the 
rest of the section suppose: 

r]{z) are i.i.d. with P{r]{z) = 1) = 1 — 9. (5.6) 

We now define the edge set S-^ for £/) to be the set of all oriented edges from z to 
z + Vi, z £ Cd, I < i < D. a sequence of points zq, . . . , Zn-i, Zn in Cd is called an 
open path from zq to 2„, and we write zq — )• z^, if there is an edge in £^ from zi 
to Zj+i and Zi is open for z = 0, . . . , n — 1. Note that Zn does not have to be open 
if n > 1 but Zq does. In Sections 6 and 7 we will employ a block construction and 
determine suitable parameters so that (x, n) G T-L^ being open will correspond to 
a certain "good event" occuring for our Poisson processes in the space-time block 
(cLX + [-KiT, KiTY^) X [nJiT, (n + 1) JiT] for appropriate L, Ki and Ji. 

Given an initial set of "wet" sites Wq C T-Lq, we say z G T-Ln is wet if — )• 2; for 
some initial wet site zq. Let Wn be the set of wet sites in Hn when all the sites in 
T-Lq are wet, and let be the set of wet sites in Hn when only G T-Lq is wet. Let 
0^ = {W° 7^ for all n > 0}. 

Lemma 5.1. (i) inf^6^„ P{x G Wn) > P{^lo) ^1 as9 ^G. 

(ii) Let Till = {{z,n) G Cd '■ z G [— r, r]''}. Then there are 65,1 > and r5,i > 
such that if 9 < 9.5,1 and r < r^^j then as N ^ 00. 

P{nl^ and n / VFn n ^7 for some n>N)^0. (5.7) 

Proof. The first result follows from well-known d = 1 results, e.g., see Theorem 4.1 
of [12]. The second result is weaker than a "shape theorem" for W^, which would 
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say the following, using the notation A' = {x' : {x',n) G A} for A C Tin- For 9 < 9, 
there is a convex set P C M"^, containing the origin in its interior, so that on il^. 



for all large n. More precisely with probability 1, if 5 > there is a random ns such 
that C n{l+6)V and D (l-(5)nPn(Tyn)' for ah n > ns. The technology 

exists to prove such a result for oriented percolation on Cd, but unfortunately no 
one has written down the details. The argument is routine but messy, so we content 
ourselves to remark that (ii) can be established by the methods used in Durrett and 
Griffeath [17] to prove the shape theorem for the d-dimensional contact process with 
large birth rates: one uses percolation in two dimensional subspaces A{mei + nej), 
^ ^ i < j ^ n and self-duality. □ 

Call sites in = Hn \ Wn dry. In Section 7, when we are trying to show that 
that S^f dies out, the block construction will imply for appropriate L and Ji, 



if (z,n) G then (clz + [-L,Lf) x [(n - l)JiT,7iJiT] is e-empty, (5.8) 



where a region is e-empty if ^f(x) = for all {x,t) in the region. This will not be 
good enough for our purposes because the space-time regions associated with points 
in = Hn \ might be occupied by particles. To identify the locations where 
there might be I's in we will work backwards in time. However in our coarser 
grid Cd, I's may spread sideways through several dry regions and so we need to 
introduce an additional set of edges for Cd- Let 8^ consist of the set of oriented 
edges from z to z — Vi for 1 < ^ < D, and from z to z + Vi — vj for 1 < i ^ j < D, 
z G Cd- 

We assume for the rest of this section that 



since we will in fact applying these results only for d > 3. Our next goal is to 
prove an exponential bound on the size of clusters of dry sites. Up to this point 
the definitions are almost the same as the ones in Durrett [13]. However, we must 
now change the details of the contour argument there, so that it is done on the 
correct graph. Let y ^ Cd with yn = n > (write y G ^/j)- In addition to P 
as in (5.6), for M > we also work with a probability P = Pn,M under which 
r]{z) = 1 for z = {z',m) G £^ satisfying m < n and \z'\ > M, and the remaining 
r/(z)'s are as in (5.6). Therefore under P the sets of wet sites {Wn} will be larger, 
although we will use the same notation since their definition is the same under either 
probability law. If y is wet put Dy = 0, and otherwise let Dy be the connected 
component in {Cd,8i) of dry sites containing y. That is, z G Dy iff there are 
zi = y, Z2, ■ ■ ■ , zk = z all in Cd so that the edge from Zi to Zj+i is in and each Zj 
is dry. Since all sites in Tio are wet, Dy C {z £ Cd : n > zd > 0}, and under Pn,A4, 
Dy C {z £ Cd '■ n > Zd > 0,\{zi, . . . , zd-i)\ < M}. We assume that uj satisfies 



d>2 



Dy{uj) is finite. 



(5.9) 
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The fact that (5.9) holds a.s. under -Pn,M is the reason this law was introduced. To 
make Dy into a solid object we consider the compact solid 

If Ry is the complement of Ry in X M+, we claim that both Ry and Ri are 
path-connected. For Ry, suppose for concreteness that D = 2> and note that for the 
diagonally adjacent points yifi) = ^(0,0,0) and y(l) = ^(1,-1,0), -Dj^(o) 
contains the edge ^({1/2} x {-1/2} x [-1/2,1/2]). For R^, if x G fi^ then there 
exists [x] G £j \ Dy such that x G [x] + Q and the line segment from x to [x] is 
contained in Ry. We first assume [x\ G T-Lk for some /c G {1, 2 . . . , If [x] is wet 
then there must be a path in Ry connecting [x] to Hq. Suppose [x] is dry, and let 
zq, zi, . . . , zk he a path in £i connecting zq = y to zk = [x]. At least one site on this 
path must be wet (else [x] G Dy), so let Zj be the first wet site encountered starting 
at Zk- Then for each i > j, Zi is dry and Zi ^ Dy (or else [x] would in Dy). Thus 
U^j{zi + Q) is path-connected, contained in Ry, and zj is connected to T-Lq by a 
path in Ry. Note that Ho C (M"^ x {0}) HRy = T-Lq which is path-connected because 
the rotated cubes making up Ry can only intersect M'^ x {0} in a discrete set of 
points (since Dy C {zd > 0}). It is here that we use d > 2. Now suppose [x] G Tik 
for some k > n. T-Lq is also connected to H-n+i by a path in Ry (assuming ^ < 1). 
This allows us to connect [x] to T-Lq and so conclude that Ry is path-connected. 

Let Ty be the boundary of Ry. To study Ty we need some notation. We define the 
plus faces of [-1/2, 1/2]^ to be [-1/2, 1/2]™ x {1/2} x [-1/2, 1/2]^^™"! , and define 
the minus faces to be [-1/2, 1/2]™ x {-1/2} x [-1/2, l/2]^-™-\ m = l,...,D. The 
images of the plus and minus faces of [—1/2, 1/2]"^ under A constitute the plus and 
minus faces of Q = A([— 1/2,1/2]*^), which are used to defined the plus and minus 
faces of Fj^ in the obvious way. Note that the plus faces of Ty will have outward 
normal Vi for some i while the minus faces will have outward normal —Vi for some 
i. 

Lemma 5.2. // (5.9) holds, then Ty is connected and bounded. 

Proof. For e > let i?^ = {x G M'^ : |x — w|oo < e for some w G Ry}. Since Ry is 
connected, so is Ry. If k{U) denotes the number of path-connected components of 
a set U, it is a consequence of the Mayer- Vietoris exact sequence with n = that 
for open sets U,V CR^ with ^7 U ^ = M^, 

k{u nv) = k{u) + k{v) - 1. 

See page 149 of [30] and also Proposition 2.7 of that reference. Applying this to the 
open connected (hence path-connected) sets Ry and Ry whose union is we find 
that Ry n Ry is path-connected. 

Finally, Ry H Ry is homotopic to F^, and therefore Ty is also path-connected. 
Boundedness is immediate from (5.9). □ 

For the next result we follow the proof of Lemma 6 from [13]. A contour will be 
a finite union of faces in Co which is connected. 



79 



Lemma 5.3. There are constants 0^,3 and 1x5,3 which only depend on the dimension 
D so that the number of possible contours with N faces, containing a fixed face, is 
at most 05,3(^5,3)^. 

Proof. Make the set of faces of Cd into a graph by connecting two if they share 
a point in common. Note that by the above definition a contour corresponds to a 
finite connected subset of this graph. Each point in the graph has a constant degree 
1/ = v{D). An induction argument shows that any connected set of N vertices has 
at most N{v — 2) + 2 boundary points. (Adding a new point removes 1 boundary 
point and adds at most v — 1 new ones.) Consider percolation on this graph in 
which sites are open with probabihty a and closed with probability 1 — a. Let be 
a fixed point of the graph corresponding to our fixed face, and Co be the component 
containing 0. If -Bjv is the number of components of size N containing 0, then 

1 > P(|Co| =N)> S^a^(l - a)^('^-2)+2. 

Rearranging, we get Bn < C^^ with C = (1 — a)"^ and /i = a~^{l — a)~^'^~'^\ 
Taking the derivative of — log a — (u — 2) log(l — a) and setting it equal to 0, we see 
that a = 1/(1^ — 1) optimizes the bound, and gives constants that only depend on 
the degree u. □ 

Lemma 5.4. If 65,^ = (Ifis^)""^^ , then 6 < ^5.^ implies that for ally = {y',n) G £j 
and all M > \y'\, PnM\Ty\ > N) < 20.5,32"^ for all G N. 

Proof. By Lemma 5.2 if Dy 7^ we see that under Pn,M, is a contour which by 
definition contains the plus faces of y + Q. Given a plus face in Ty if we travel the 
line perpendicular to M*^ x {0} and through the center of the face, then we enter and 
leave the set an equal number of times, so the number of plus faces of Ty is equal 
to the number of minus faces. Thus, if the contour has size N there are A/2 
minus faces. It is easy to see that a point of Wj adjacent to a minus face associated 
with a point in Vj+i must be closed for otherwise it would wet the point in t^+i 
(recall the outward normal of a minus face is —Vi for some i). The point of Wj that 
we have identified might be associated with as many as D minus faces, but in any 
case for a contour of size A there must be at least N/2D associated closed sites. 
Taking 9 < (2^/5 3)"^^^ using Lemma 5.3 to bound the number of possible contours 
containing a fixed plus face of y + Q, and summing the resulting geometric series 
now gives the result. □ 

It follows from the above and an elementary isoperimetric inequality that there 
are finite positive constants C,c such that for all y = {y',n) € £j and M > \y'\, 

ii0< 65.4 then Pn,M{\Dy \ > A) < C7 exp(-cA(^"^)/^) for all A G N. (5.10) 

Now fix r > and let Bn be the dry sites in Ti^^^ connected to the complement of 
^m=n/2^^^^ ^ path of dry sites on the graph with edges £1, where as for open 
sites the last site in such a path need not be dry. 
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Lemma 5.5. If < 6*5.^ then 

P{Bn 7^ infinitely often) = 0. 

Proof. Let M > n{r + \/2D). We couple the iid Bernoulli random variables {r){z) : 
z S Cd} (under P) with the corresponding random field f] (under P = Pn,M) so 
that 

r]{z) = fj{z) Vz = {z',m) where \z'\ < M or m > n. 

We claim that z E ^m=n/2'^n^ ^ ^® ^ ^' -'^^ clearly 

suffices to fix z = {z',m) G which is wet for fj and show it is wet for ij. A 
path of sites Zi = (z^,i), i = 0,...,m with edges in £^ from Tio to z satisfies 
maxj<m l-z^l < rn + \/2Dn < M. This is because the edges in £^ have length at most 
\/2D. Therefore if the sites in the path are open in fj, then they will also be open 
in rj. This proves the claim. 

Next note that if y G HlT^^, then y £ Bn for r] iS y ^ Bn for fj. This is because 
the path of dry sites connecting y to the complement of ^m=n/2^^'^^'^ taken 
to be inside Tin and so we may apply the claim in the last paragraph. It now follows 
from the above bound on the length of the edges in £^ that 

\Dy\ > -^j. 

The number of sites in T-Cn ^^ is at most Cn'^, and the bound in (5.10) shows that 
P{Bn 7^ 0) < ^ P {y G Bn) Is summable over n. □ 

Remark 5.1. We will prove in Section 7 that if wet sites have the property in 
(5.8), and the kernels p{-) and q{-) are finite range, then for an appropriate r > 0, 
Bn = 9 will imply that on all sites in [—CL,d'''n, CL^dfn]'^ will be vacant at times 
t G [(n — l)JiT,nJiT]. This linearly growing dead zone will guarantee extinction of 
the 1 's. 
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6 Existence of stationary distributions 



With the convergence of the particle system to the PDE estabhshed and the per- 
colation result introduced, we can infer the existence of stationary distributions by 
using a "block construction" . Recall that our voter model perturbations take values 
in {0, l}^^*^ and so our stationary distributions will be probabilities on this space of 
rescaled configurations. We begin with a simple result showing that for stationary 
distributions, having some I's a.s. or infinitely many I's a.s. are equivalent. Let 

ICI = E.C(^)- 

Lemma 6.1. If i' is a stationary distribution for a voter perturbation, then 

1^1 = oo u — a.s. iff \^\ > v — a.s. 
Proof. It suffices to prove 

z^d^l < oo) > implies z/(|^| = 0) > 0. (6.1) 

Assume first that the configuration is a trap. Then if |,^o| = K < oo, (1.5) 
shows the sum of the flip rates is finite and so it is easy to prescribe a sequence of 
K flips which occur with positive probability and concludes with the state. By 
stationarity we get the implication in (6.1). 

Assume next that is not a trap, which means (^f (0, . . . , 0) > 0. We claim that 
z^d^l < oo) = 0, which implies the required result. Intuitively this is true because 
configurations with finitely many I's have an infinite rate of production of I's. One 
way to prove this formally is through generators. Let Vt^ be the generator of our 
voter perturbation, Vt^ be the generator of the voter model in (1.3) and for i = 0, 1 

^Mi) = E = 1 - i)Ei9mx + Y^), . . . ,C(x + y^o)))(^(^x) _ ^(^)). 

Here if) will be a bounded function on {0, l}^'' depending on finitely many coordi- 
nates, and we recall that is ^ with the coordinate at x flipped to 1 — i{x). Recall 
that ie{ex) = i{x) for i G {0, l}'^\x G Z'^. For ip as above define tpe on {0, ly^" 
by i^eiie) = ^{i). Then by (1.5) and (1.12), 

n'^l>e{^e) = {e-^ - e^^)n.i^{i) + OoV^(0 + ^MO- (6-2) 
For < r < i?, let A{r, R) = {x e : r < \x\ < R] and 

i^rAO = l(eU(r,i?) =0), {0,1}^'. 

Considering two cases x G A{r, R) and x A(r, R) we have 

if e(x) = then Vr,fl(e) - A,r{0 < 0- (6.3) 

Since ipr,R{C^) — V'r,_R,(C) = 1 only if x is the only site in A{r, R) where ^(x) = 1, we 
have 

^vAMO < 1, ^O^rMO < Ikolloo- (6.4) 
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Choose A so that P{Y* < A) > 1/2, where Y* is as in (1.8). Fhpping a site from 
to 1 cannot increase tpr,R, and ipr,Ri(,) = 1 imphes ^(x) = for ah x G A{r,R), so 
we have 

niA,R{0<- J2 {l-^{xM{0,...,0)P{ax + Y') = OfoTl<i<No)^Pr,R{0 

x&A{r,R) 

< - ^^^"'^••'°V r,fi(OI^(r + A, - A)|. (6.5) 

The stationarity of i' implies, see Theorem B.7 of Liggett [.33], that if -0 = '>pr,R then 
jQ'^'iped^ = 0- Using (6.2), (6.4) and (6.5), and noting that 

J ip^dv = i/(^ = on A{er, eR)), 

we have 

< {e-' - ef) + ||<7g||oo - + A, - \)\u{i ^ on ^(er, eR)). 

Rearranging this inequality we get 

v{i = on A{er,eR)) < ~ ^ ^ ^, 

^ ^ - gf(0,...,0)|yl(r + A,ii-A)| 

(recall (/i(0, . . . , 0) > 0). Letting i? — t- cxo we conclude that = on A{er, oo)) = 
0. In words, for v-a.a. configurations there is a 1 outside the ball of radius er. As 
this holds for all r < oo, there are infinitely many ones with probability 1 under u. 
□ 

Assumption 1 and (1.63) are in force throughout the rest of this section and we 
drop dependence on the parameters w, Vi, Lq, Li, vq, etc. arising in those hypotheses 
in our notation. We continue to work with the particle densities D{x,S^) using the 
choice of r in (4.6). We start with a version of Lemma 4.2 which is adapted for 
proving coexistence. We let 

La = 3 + Lo VLi. 

Lemma 6.2. There is a Cg ^ > and for every r] > 0, there are T,^ > 1 and 
eg 2{'n) > so that for t S \Trj,Cg glog(l/e)] and < e < e g 

^0 has density in [vq +r/, ui — rj\ on [—L2,L2f, 

then 

P(^f has density in [u^ — r],u* + r]] on [—wt,wt]'^\S^Q) > 1 — e'*^^. 

The proof is derived by making minor modifications to that of Lemma 4.2 and 
so is omitted. We will always assume r/ > is small enough so that 

0<VQ + r]<u^:— r]<u*+r]<vi — r]<l. 



83 



The one-sided versions of the above Lemma also hold (recall Lemma 4.2 on which 
the proof is based is a one-sided result), that is, with only one-sided bounds on the 
densities in the hypothesis and conclusion. 

Theorem 1.15. Suppose Assumption 1 and (1.63), and let rj > 0. If e > is small 
enough, depending on rj, then coexistence holds for the voter model perturbation, the 
nontrivial stationary distribution v may be taken to be translation invariant, and 
any stationary distribution such that 

Yl ^(^) = E (1 - ^(^)) = O) = (6-6) 
satisfies ^{(^{x) = 1) G (u* — ri,u* + rf) for all x € eJJ^ . 

Proof. We use the block construction in the form of Theorem 4.3 of [14]. This result 
is formulated for D = 2 but it is easy to extend the proof to -D > 3, and we use this 
extension without further comment. Recall Q{r) = [— r, r]^ and Q^{r) = Q{r)r\eL'^ . 
Let \J = (Cg 2/2) log(l/e), L = wU/{aQD + 1), where ao > is a parameter to 
be chosen below, and /* = [u^, — 77/4, «* -|- r?/4]. Next we define the sets H and 
which appear in the above Theorem. Let 

H = {^e {0, ly^" : ^ has density in /* on Q{L)}, 

that is, if Qe = [0,0^)'^ n eZ'^ then the fraction of occupied sites in x + Qs is in 
/* = [u* - r//4, u* + 77/4] whenever x G a^Z^^ n [-L, if. If L' = L + 1, then G H} 
depends on Here we need to add 1 as the cubes of side with "lower 

left-hand corner" at x S [—L,L]'^ will be contained in [—L',L']'^. This verifies the 
measurability condition in Theorem 4.3 of [14] with L' = L + 1 in place of L which 
will affect nothing in the proof of Theorem 4.3. 

Let be the event on which (a) if = ^, then has density in /* on Q{wU) 
and (b) for all z G Q%wU + 1) and all t < U, X^'^ C Q((w + bo)U + 1), where 
60 = 16(3 + d)/CQ 2- Note that 

G^ GcT(^ASf|[o_^]x^Z''^^ox[o,i]'^jl>l[o,c/]x£Zd :y e {{w + bo)U) + (6.7) 

=g{Q{{bo + w)u + i)) X [o,u]) 

Informally, g{R) is the cj-field of generated by the points in the graphical represen- 
tation that lie in R. The above measurability is easy to verify using the duality 
relation (2.17). 

Consider now the Comparison Assumptions prior to Theorem 4.3 of [14]. In our 
context we need to show 

Lemma 6.3. For < e < eg s{t]): 

(i) if Co £ H, then on G^s, ^fj has density in I* on a^Lv^ + [— L,L]'^, 1 < i < D , 

(ii) if^GH, then P{G^) > 1 - e°-°l 
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Proof. By assuming e < £i{r]) we have U > and L > L2. Using the definition 
of L and the fact that \v',j\ < \\v[\\2 = \/D — 1 one easily checks that 

aoLv'i + [-L, Lf C [-wU, wUf for i = 1, . . . , D. (6.8) 

Part (a) of the definition of now gives (i). By Lemma 4.7 with parameters 
L = wU + 1, 2b = 60, c = Cq 2/2 and T' = U, and Lemma 6.2, for ^ G we have 
for e < EQ^iv), 

P{Gl) < e-°5 + c'^iwU + l^iU + i)^((^oCg 2)/i6)-2-rf 
< e-°5 + c(log(l/e))'^+ie < 

where the last two inequalities hold for small e. We may reduce Cq 2 to ensure that 
b = bo/2 satisfies the lower bound in Lemma 4.7. This proves (ii). □ 

Continue now with the proof of Theorem 1.15. Let e < eg 3 and define 

Vn = {{xyn) E Hn ■ CnU density in /* on aoLx + [—L,L]'^}. 

(To be completely precise in the above we should shift uqLx and aoLv'^ to the point 
in eZ"^ "below and to the left of it" but the adjustments become both cumbersome 
and trivial so we suppress such adjustments in what follows.) If we let 

Ry,n = {yaoL, nU) + Q((6o + w)U + 1) x [0, U], for {y, n) G Co 

and 

^ ^ r 2(6o + ^;)(aoZJ + l) - ^ 
aQW 

then Ry-i^^m H Ry2,n = if |(yi> "i) — {y2,n)\ > M. Since G{Ri), 1 < i < k are inde- 
pendent for disjoint i?i's, Lemma 6.3 allows us to apply the proof of Theorem 4.3 
of [14]. This shows there is an M-dependent (in the sense of (5.5)) oriented perco- 
lation process {W^} on Cd with density at least 1 — e'^^ such that Wq = Vq and 
Wn C Vn for all n > 0. We note that although a weaker definition of M-dependence 
is used in [11] (see (4.1) of that reference), the proof produces {VFn} as in (5.5). By 
Lemma 5.1 with r = rr^ i and 6 = e''^^, if e < si{rj), then 

lim inf PiiniJ has density in I* on a^Lx + [—L,Vf) (6.9) 

> (l-|)p(OGVo). 

We will choose different values of uq to first prove the existence of a stationary 
law, and then to establish the density bound for any stationary distribution. For 
the first part, set ao = 3 and take {£,q{x) : x G eZ'^} to be iid Bernoulli variables 
with mean u = (n* +u*)/2. The weak law of large numbers implies that if e is small 
enough 

P(Cg has density in /* on [-L,Lf) > (6.10) 
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Since ao = 3, L > 3 and \x — y\ > \\x — y\\2/VD > 1 for all x / y G H'^, 
{aoLx + [—L', L'Y" '■ X G "H^} is a collection of disjoint subsets of for each n. This 
and the measurability property of G H} noted above shows that if < e < eo(^) 
then {Vri} is bounded below by an M-dependent (as in (5.5)) oriented percolation 
process, {W,!^^}, with density > 1 — e'^'^ starting with an iid Bernoulli (1/2) field. 
Having established that our process dominates oriented percolation, it is now routine 
to show the existence of a nontrivial stationary distribution. We will spell out the 
details for completeness. 

Lemma 6.4. Assume = 3 and {5,q{x) : x G el/} are as above. There is an 
£()_^{r]) > so that for any e G (0, e6'^^(ry)) and any /c G N there are ti{k,e), 
Mi{k,e) >0 so that fort > ti, 

P{Y1 et{x)>kand i_^^-(x)>A;) >1-^. 



|x|<Mi \x\<M-_ 



Proof. As in Theorem A. 3 of [1 I] for G N there are no, ^O) Mq G N and zi, . . . , z^^^ G 
Q{Mq) satisfying \zi — Zj\ > 3M + 2£o + 1 for i / j, such that for n > uq with 
probability at least 1 — 

WU'' n Q(z„4) / for i = 1, . . . ,4A;. (6.11) 

The above implies there are (7(^^jj)-measurable y^ G Q^{zj,£Q) such that 

^f,lj has density in /* on SLyj + [-L, L]"^, j = 1,... ,4k. (6.12) 

This proves the result for t = nil. Intermediate times can be easily handled using 
Lemma 6.2 and the finite speed of the dual (Lemma 4.7). Those results show that for 
a fixed e < eg 2 * — ("'0 + l)f^) if we choose n > no so that t G [(n+l)f/, {n+2)U] 
(use < 2U = Cq 2 log(l/e) in applying Lemma 6.2), then on the event in (6.12) 
we have 



P(^f has density in /* on 3Lyj + [— L, L] , 



and * G Q{3Lyj,L' + boU) for all x G SLyj + [-L', L'Y and s G [0, t]\ilj^ 



>l_e-05_ci(log(l/e))'^+ie>i. 

where in the last we may have needed to make e smaller. 

Our separation condition on the {zj} and L > 3 implies that Q{3Lyj, L' + boU), 
j = 1, . . . ,4A: are disjoint and so the events on the left-hand side are conditionally 
independent as j varies. Therefore a simple binomial calculation shows that 

P{\{j < 4k : 5,1 has density in /* on 3Lyj + [-L, Lf}\ > k) 

^(^-^)('-l)^^-F 

Here the first 1 — comes from establishing (6.12) and the second 1 — ^ comes from 
the binomial error in getting fewer than k points with appropriate density at time 
t. Since the above event implies the required event with Mi = 3L(Afo + ^o) + -^^ we 
are done. □ 
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Fix e < eg 4. By Theorem 1.1.8 of [32] there is a sequence t„ — )• oo s.t. 

Jp" G ■) ds — )• z/ in law where z/ is a translation invariant stationary distri- 
bution for our voter perturbation. Lemma 6.4 easily shows that there are infinitely 
many O's and I's z/-a.s., proving the first part of Theorem 1.15. 

Turning to the second assertion, by Lemma 6.1 and symmetry it suffices to show 
that for e < e2(^) and any given stationary i/ with infinitely many O's and I's a.s. 
then 

sup fJ^{(,{x) = 1) < u* + rj. 
Start the system with law v. We claim that 

Lemma 6.5. There is a a {^q) -measurable r.v. xq € eZ*^ such that = on 
Q^{xq,L) a.s. More generally w.p. 1 there is an infinite sequence {xi : i G Z+} 
of such random variables satisfying \xi — Xj\ > 2L + 3 for all i ^ j. 

Proof. To see this condition on ^q, choose xq so that Co(^o) = and note that if Rf 
is the first reaction time of the dual X^'^, the event "^g ^ on Q^{xo, Ly^ occurs if 
for all X e Xq + [—LjL]''-, Rf > 1, X^'^ = xq, and sup^<]^ iXg'"^ — x\ < 1. Call the 
last event A{xq). The last condition has been imposed so that if |xo — xi| > 2L + 3 
then the events ^(xq) and ^(xi) are (conditionally) independent. Clearly they have 
positive probability. Given our initial configuration with \{y : ^^{y) = 1}\ = oo a.s., 
we can pick an infinite sequence Xj, i G N, with ^o(^«) — ~ > 2L + 3 

when j > i, so the strong law of large numbers implies that at time 1 there will be 
infinitely many Xj with ^f(x) = for all x G Q^{xi,L). By stationarity this also 
holds at time 0. □ 

Now condition on ^q, shift our percolation construction in space by xq, set oq = 
(2D)-^ and only require the density to be at most u* + rj/A in our definition of Vn 
which now becomes 

Vn = {{x, n) G T-Ln ■ £,nU density at most u* + r]/4 on xq + c^x + [— L, L]"^}, 

where we recall from (5.2) that cl = L/{2D). (Here we are using the one-sided 
version of Lemma 6.2 mentioned above, after its statement.) Then G Vq and the 
one-sided analogue of (6.9) shows that if e < e3(r/), then 

lim inf P(x G K) > 1 - ?• 

n^oo (a;,n)eW™ 4 

Recall from (5.2) that Ua;g-^/^xo + CLX-|- [— L, LY = M'', so this implies for any x G M'' 
and n large enough, 

P{inU density at most -|- ^ on x -|- [— L, Vf) > 1 — ^, 

and so by stationarity 

v{C has density at most u* + - on x + [-L, Vf) > 1 - - for all x G M.'^. 
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To complete the proof, run the dual for time ts (te as in (4.6)) and apply Lemma 4.5 
with n = + I to see that for x G eZ'^ and e < 63(7]) A £4 ^{rj/S), 



Ke(^) = i) = p(ct(x) = i) 

< PiRi < Q+E{P{Ri > t„etS^) = ii^o)) 

< (1 - e-^*'^) + EiPieoiBtf) = l\eo)) 
<c% + | + n* + | + |<n*+r?, 

where e is further reduced, if necessary, for the last inequality. □ 
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7 Extinction of the process 
7.1 Dying out 

Our goal in this section is to show that if /'(O) < and |^q| is o{e~'^), then with 
high probabihty wih be extinct by time 0(log(l/e)). Throughout this Section we 
assume that < e < Eq and that (1.65) holds, i.e., g^O, . . . , 0) = for < e < Eq 

Recall from (3.30) the drift at ex in the rescaled state G {0, ly^'' (recall the 
notation prior to (1.2)) is 

d,{ex,C,) = (1 - ^{x))hl{x,0 - ax)hl{x,0, 

and define the total drift for l^^j < oo by 

M^s) = Y.de{ex,Ce). (7.1) 

X 

Recall from (1.7) and (2.6) that 

hlix, = Erigliax + F^), . . . C(x + F^"))), (7.2) 

where Ey denotes the expected value over the distribution of {Y^, . . . Y^°), and also 
that 

c* = c*{g)= sup ||5il|oo + ||ffolloo + l, Cb = c*No. (7.3) 

0<e<eo/2 

It will be convenient to write 

Csiex + eY) = (C(x + eY^), ...^{x + eF^")). 

If Tit is the right-continuous filtration generated by the graphical representation, 
then ^ 

\et\ = \eo\ + M!+ [ Mes)ds, (7.4) 

Jo 

where is a zero mean L^-martingale. This is easily seen by writing (x) as a 
solution of a stochastic differential equation driven by the Poisson point processes 
in the graphical representation and summing over x. The integrability required to 
show A-I^ is a square integrable martingale is readily obtained by dominating l,^"^! by 
a pure birth process (the rates Cg are uniformly bounded for each e) and a square 
function calculation. 

Lemma 7.1. For any finite stopping time S 

e-^o'm < E{\es+t\\ns) < e^^'iesl 

Proof. By the strong Markov property it suffices to prove the result when S = 0. 
The fact that de{£x,^e) ^ "Ibolloo implies ^'^(Cf) ^ "llffollool^f I- It follows from 
(1.65) and (7.2) that 

de{ex,ie)<\\g{\\oo E^(^' = 2/-^)- 

y:li{y)=l i=l 



89 



Summing over x and then y, we get V'e('^f) ^ -^olbillool'^f | and (recalling (7.3)) the 
desired result follows by taking means in (7.4) and using Gronwall's Lemma. □ 

Let i^^'" be the voter model constructed from the same graphical representation 
as by only considering the voter flips. We always assume .^q''^ = .^q. 

Lemma 7.2. //cy.g = 4(2iVo + l)c* then 

Proof. Let (ex + eY) = (^f (ex + eY°), . . .,Cl{ex + eY^°)), where Y^ = 0, Y is 
independent of and note that in contrast to Y, Y contains 0. Let 

Deivo, m> • • • VNo) = -vogoivi, • • • , Wo) + (i - • • • , vno), 

and note that 

Eimes) - Mes'')\) <e(j2 imtiex + sy) - DM''iex + em) 

X 

< 2||L>,||oo^(V[ max il{ex + eYi)y^f{ex + eYi)] (7.5) 

X 

xl{4^(ex + ey)/e°(ex + ey)}), 



because for fixed x if the latter summand is zero, so is the former, and if the latter 
summand is 1, the former is at most 2||L)£||oo. 

Let Xt = Xl''^ ^ t S [0, s] be the dual of starting at (zq, . . . , zatq) = ex + eY at 
time s and let Rm, m > 1 be the associated branching times. We claim that 

£;( [ max esiex + eY')] (ex + eY) + if {ex + eY)]] (7.6) 



To see this, note that: 



(i) if i?i > s, then there are no branching events and so {Xt,t < s) is precisely the 
coalescing dual used to compute the the rescaled voter model values ^|''^(ex + eY). 

(ii) In the case Ri < s, if ^g(Xf) = for all £ G J(s) then (ex + eY') = for 
< i < Nq because working backwards from time to time s, we see that no site 
can flip due to a reaction, and again we have .^''(ex + eY) = ^^'^{ex + eY). 

Similar reasoning and the fact that the dual {X^'-' ,j G J^{t)) of the voter model 
i"'^ with the same initial condition z satisfies J'^{t) C J{t) for all t < s a.s., shows 
that 

max if{ex + eY')]m{ex + eY)^if{ex + eY)}) (7.7) 



<e[y1 eo{x!)i{Ri < s} 
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If Eq denotes expectation with respect to the law of ' when x = then, using 
(7.6) and (7.7), we may bound (7.5) by 

Bounding by the dominating branching random walk X, and using | J(^i)| = 2A'^o + l 
and P{Ri < s) = 1 — q-^* (^0+^)3 ^ ggg ^j-^g expected value in the last formula is 
at most 

\eo\Ei\Jis)\l{R, < s}) < |e§|e^*^o^S(|J(i?i)|l{i2i < s}) 
< {2No + l)|e§|e'=*^o"(l - e-^*(^o+i)^) < {2No + l)|e§|(e"*(^°+^)'^ - 1), 

which proves the desired result. □ 

For the next step in the proof we recall the notation from Section 1.8. We assume 
Y is independent from the coalescing random walk system {B^ : x S Z'^} used to 
define t{A) and t{A,B). RecaU from (1.88) and (1.85) that under (1.65) 

0^nO)= /3(5)P(r(y^)<oo,T(y^,{0}) = oo)-5(5)P(T(y^U{0})<oo). 

For M > define 

0%i= E /3e(5)P(r(y^)<M<r(y^{0}))-4(5)P(r(y^U{0})<M) . 

It follows from (1.89) that (with or without the e's) 

mS)\ + 14(5)1 < 22^"(ll5f Hoc + ll^^lloo) < 22^°c*(5) (7.8) 

(recall here that gf = gf by our ei = 00 convention). It is clear that limA//^oo,e-5>o ~ 
6, but we need information about the rate. 

Lemma 7.3. There is a £'j,3{M) ]^ (independent of the gf) so that 

\9i - 9i\\oo + Ibo - ^olloo + c*{g)e7.3{M) 



91j -9\< 22^" 



Proof. Define 

ej.3{M) = ^ sup {P{M < r < 00) : r = T(y^), t{Y^ U {0}),t{Y^, {0})}, 
seVNo 
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and note that ej,i,{M) J, as M t oo. Using (1.90) and (7.8) (the latter without the 
e's) we have 

\e%i-e\< \j3^{s)-Ks)\ + \Us)-ks)\ 

+ [\kS)mr{Y^)<M<r{Y',{0}))-P{r{Y')<^ = r{Y',m)\ 

+ \6{S)\P{M < t{Y^ U {0}) < oo) 

<22^°[||5f - <7l||oo + U - 5ol|oo] + 22^0c*(5)e7.3(M). 
The result follows. □ 

To exploit the inequality in Lemma 7.2 we need a good estimate of -E'(V'e(^f''^)) 
for small s. 

Lemma 7.4. There is a constant cy,^ (independent of gf) such that for e,5 > 0, 

E{Me/)) = +r?7.^(e,5), (7.9) 

where \v7,^{e, 5)\ < C74C* {9)5-"/^ ICol^e''. 

Proof. As usual we assume Y is independent of Summability issues in what 
follows are handled by Lemma 7.1 (and its proof) with gf = 0. The representation 
(1.83) and (7.1) imply that 

Ernes')) =Y.MS)E^oiS) - 1{S)EI{S) 
s 

E^,{S) = e{{1 - es'\ex))lles'\ex + eY^)) 

EfiS) = E[e/iex)Ylcf{ex + eY')). (7.10) 

We will use duality between ^^'"^ and {13^} (see (V.1.7) of [.'>2]) to argue that 
E^oiS) « |eg|P(r(y^) < < r{Y^^ {0})), / 5 C {1, . . . , iVo} 
EliS) « \eo\P{r{Y'', {0}) < fe-2) all 5 C {1, ... , iVo}- 

Beginning with the first of these, note that duality implies (recall Y^ = 0) 

E^s) = e[{i - eo^(eBf,_.)) n^o(^^r- + y', M) > 

= E E(jleo{eB's^-V)l{r{x + Y'Ax})>6e^'y 

-E^( n )i{T(x + y^w)>fe-2}) (7.11) 

x&Z'i iesu{o} 
= Si — II2. 
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If t{x + Y^) > ^i^gre are i / i G 5 so that t{{x + Y'}, {x + Y^}) > fe"^. 
If we condition on the values of the Y^,Y^ in the next to last line below, 

E(Jleo{eB^:-S)l{r{x + Y^) > de^'}] 
^ E E EieoieB^s^SmeB^s^Z'nM^ + Y'}, {x + Y^}) > 5e-'}) 

x&Z'i l<i<j<No 

< E E eo{ew)eo{ez) E ^(^r- = ^r- = 

«)6Z'' zeZ'' l<i<j<No x&Zd 

T{{x + Y'},{x + Y^})>Se'^}) 

< E E ^o(^w^)eg(e^) E = «^ - 2 - + IS) 

< NoiNo - l)|^gpc(l + 2fe-2)-d/2^ (7,12) 

where the local central limit theorem (e.g. (A. 7) in [ii]) is used in the last line. A 
similar calculation shows that 

^2 < |^§Pc(l + 2fe-2)-^/2. (7.13) 

To see this, note that r(x + Y^ , {0}) > 6e~^ implies that for io G S (this is where 
we require S non-empty) t{{x + Y^°},{x}) > fe~^) and we may repeat the above 
with i = io and j = 0. Returning to the study of Si, taking any io £ S we have 



E(j[eoieB^^-V)l{r{x + Y'')<6e-'<T{x + Y',{x})}) 
Y E[eo{ex + sBj;^,)l{r{Y') < 6e~' < r{Y^,m}] 



= \eo\P{r{Y') < 5e-2 < ^(y5 |o|))_ (7 

Together (7.12) and (7.14) bound Si. Using this with (7.13) in (7.11), we con- 
clude that 

E'oiS) = \eo\P{r{.Y^) < fe-' < r(y^{0})) +r?i(e,5,5), (7.15) 

where \rii{e,6, S)\ < cA^q |^Qp5~'^/^e'^. A similar, and simpler, argument shows that 
for Sc{l,...,Aro}, 

EliS) = |e§|^'(r(y^U {0}) < 5e~^) + r^2{e,6,S), (7.16) 

where \r]2{e,6,S)\ < cNo{No + l)|^g|2<5-'^/2e'^. 

Now use (7.15), (7.16) and the fact that /3e(0) = (by (1.65)), to obtain (7.9) 
with 

|r/7.4(^,5)| < E(l/^^(^)l + MS)\)cNo{No + 1)6-'^' Idol's'- 
s 

Finally use (7.8) to complete the proof. □ 
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For < ?7i < 1, let r(r?i) = T,{t]i) = mf{t > e'^i : l^f.^^J > (e-^+^)'^e''i} and 
note that T(r/i) — e''^ is an ('H()-stopping time. 

Lemma 7.5. There is a C7.5 so that if rn G (0, 1), then for all s > e''^ 

Ei^iCsms-e^.) < [ei,,,-2 + crse^'] a.s. on {r(7?i) > s}. 

Proof Let 6 = e'^K If < (e"^+^)'^e''i, then Lemmas 7.2 and 7.4 imply 

^(V'(a)) = iai^L-2+v(£) with 

For the second term we used the bound on |,^g|. The result now follows from the 
above by the Markov property and the definition of T{rji). □ 

Lemma 7.6. Let /3,ri2 G (0, 1]. There is an ej.eiP^V^) G (0, 1), so that if < e < 
fTy.g and 9 = /'(O) < —rj2, then |^q| < e~'^~^^ implies 

P{\C!\ > 0) < Ge^^^'e^/^ ^ > Miog(i/^), 

Proof. Let A < 7/2/2, r/i = /3(2 + T = Te{i]i) and 5 = An integration by 

parts using (7.4) shows that for t > 6, 

e'^^^^^le^ATl = e''m + l{r < T}[Ae^ne.1 + e'^A^dr + iVf , 
where A^"^ is a mean martingale. Since {r < T} G Tir-s, we have for (5 < s < t 

= / ii;(l{r <T}e^^i5;(A|e,^|+Ve(e.")l^r-5)|^s-5)dr. 

J s 

Using Lemmas 7.1 and 7.5 the above is at most 

E{f^l{r < T}e'-j{e)\C.s\dr\ns-s), (7.17) 

where 7(e) = Ae^"*^ + 6l|^_2 + 07.5^. Recall <5 = e"^^ and 6 = /'(O) < -m- By 
Lemma 7.3 and the uniform convergence of the to gi there is a ei(/3,?72) > so 
that if < e < £1, then 



7(^) < f e-^^^" - r?2 + 22^0 - g^\\oo] + 2'^^c*e,.3{e'>^e-') + cr.,e^'^ 



4 = 



< -r?2/4 < 0. 
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We assume < e < ei in what follows. Since the bound in (7.17) is therefore non- 
positive and our assumption on |^q| implies T > 5, we may use Lemma 7.1 and the 
fact that (5 < 1 to see that for t > 5, 

EiietATle^^'""^^) < e^'^^(l^ll) < e(^+'='')^|eol < e^^+'^'-e^"'^. (7.18) 
Now l^f I > 1 if it is positive so 

P{\C!\ > 0) < £;(|eMT|e^^'^^*^l{r > *})e"^* + P{T < t). (7.19) 

Let t > (2d/r]2)^og{e~^) and use (7.18) with A = r/2/2 to see that the first term is 
at most 

To bound P{T < t), we note that ICt^sI > {e~^^^fe'^^ if T < cx), so 

P{T <t)< Eil^T^slHT < t})(e~^+^)~V''i. 

By making ei smaller, depending on /3, we can assume that (2^/772) log(e^^) > e"^^ = 
6. Let S = {T — 6) A{t — 6), note {T < t} £ Tis, and use the lower bound in Lemma 
7.1 with A = to conclude the first inequality in 

mr-slHT < t}) < e'^>''E{\eTM\) < e^'^^+'^^e^"'^. 

The second inequality comes from (7.18) with A = (recall that t > 2d/r]2) log(e^^) > 
e^i = (5) and 5 < 1. Using the last two equations with (7.20) in (7.19), we conclude 
that 

where the definition of i]i is used in the last line. The result follows. □ 



7.2 The Dead Zone 

For the remainder of this Section we suppose (1.63), (1.65) and Assumption 2 are 
in force and — /'(O) > t/2 S (0,1]. We also assume that p{-) and q{-) have finite 
supports. More specifically, Rq £ N satisfies 

{x G Z'^ : p{x) > 0} C [-Ro, Ro]"^ and {x G Z'^^" : q{x) > 0} C [-Rq, Rof^". 

(7.21) 

In order to connect with the percolation results from Section 5 we need certain space- 
time regions suitable for applying Lemma 4.2 to decrease particle density, Lemma 4.7 
to control the spread of duals, and Lemma 7.6 to actually kill off particles. Recall 
that Q%r) = [-r,r]'^ n (eZ^). For Jq < Ji G N, < ij; < 1, A, K > 1, and T > 
define regions ^^{Jo, Ji,w, A, K) = Dq U Pi, where 

^0 = u/li(g^((K - 3A)T) X [(i - 1)T, jT]), 

^1 = ^i=jl{Q'{{K + jw -{w + A)Jq)T) X [jT, [j + l)T]). 

For help with the definition consult the following picture: 
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[K + {Ji-l- Jo)w - JoA]T 



Ji = 7 



Pi 



[K - JoA]T 



Vn 



[K - A]T 



The speed it) > is as in Assumption 2 (and may be assumed to be < 1), and 
T = A,i,2 log(l/e) is the same as in (4.6). For the regions VqjVijS we take 



Jn 



2d 



mAi.2 



+ 1, 



, 8(2d + 3)) , 2n 



A. 



4.2 



K = 2 + AJq, and Ji = Jq + 1 + 



K + A,h 



w 



(7.22) 



The choice K = 2 + AJq impUes that Q^{2T) x { JqT} is the "top" of Vq and 
the "bottom" of T>i. The choice of Ji imphes 

the top of T> contains Q^{2KT) x { JiT} and is contained in (7.23) 
Q%{2K + l)r) X { JiT}, and V contains the region Q^{2T) x [0, JiT]. 

Recah from Section 5 that a region C in X M+ is e-empty iff (x) = for 
ah {t,x) G C, where is our voter model perturbation as usual. If A C M"^ let 

Lemma 7.7. There exist £7.7, C77 > depending on ui,U2,w, C2, C2 (from Assump- 
tion 2) and ro,74.2,?72 such that such that if < e < £7,7 and 



miKT)) = 0, 



then 



PyV{Jo, Ji,w, A, K) is e-empty j > 1 
and with probability at least 1 — cq^^e'^ , 

for allj = l,...,Ji,{x,u) G {Q' {{K + Jiw + A{Ji - j))T) x [{j-l)T,jT]) 



(7.24) 
(7.25) 



and te[0,u- {j - 1)T],X^''' C Q'{{K + Jiw + ^( Ji - j + l))r). 



(7.26) 
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Proof. We begin with some notation for describing events in which the dual process 
is confined to certain space-time regions. For j > 1 and < r < s let TT{j,r,s) be 
the event 

X^'"" C Q%sT) V X e Q'{rT), u G [(j - l)r, jT], and te[0,u- {j - 1)T]. 
On rj'(j, r, s), duality and (1.65) imply 

?(,_i)t(Q'(sT)) = implies Q%rT) x [(j - l)r, jT] is e-empty. (7.27) 

Step 1. We first check that Vq is empty with high probability. For j S {1, . . . Jq} 
we bound the probability of Tt^J, K — jA, K — [j — V)A) by using Lemma 4.7. If we 
set c = ^4,2, = jT, L = (K — jA)T and 2b = A, then evaluating q in the lemma 
we obtain 

g = (^-2)Ae-2>2d+l 

o 

if > (2(i + 1). Hence the bound on in Lemma 4.7 gives us 

PiVrij, K - jA, K - {j - l)A)) > 1 - c'^(ir - jA)''e'' (7.28) 

for £ < £4.6(^/2) such that e^^ >2d+l and £(^4.3 log(l/e) + 1)'' < 1. 
By (7.27), on the intersection 

^jUTT{j,K-jA,K-{j-l)A), 

for each j G {1, . . . Jo}, if ef,_i)T(Q'((^ " (i " 1)^)^)) = then Q'{{K - jA)T) x 
[(j — 1)T, jT] is e-empty. Iterating this, (7.24) and (7.28) imply that for some positive 

P {Vo is e-empty) > 1 - c^JoK'^e'^ e < Eq. (7.29) 

Here, and throughout the proof, Eq will denote a positive constant depending only 
on our fixed parameters including vq. 

Step 2. By taking e small enough we may assume that (recall Lq is as in As- 
sumption 2) 

2 + Lq<KT<{K + wJi)T < £-■001/'^. (7.30) 

For j S {1, . . . , Ji — Jo}, on account of (7.24), we may apply Lemma 4.2 and the 
Markov property Ji — Jq times and conclude that for e < £4.2, 

P(^|r has density at most e^'^-S in Q'{{K + wj)T) for j = 1, . . . , Ji - Jq)) 
> 1 - (Ji - Jo)e-°^ 
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When the above event occurs, for any j G {1,2, . . . , Ji — Jq}, (recall that = 

^^T{Q-{{K + wj)T)) 

< E ^jTix + Qs)+ E ^JTix) (7.31) 

xea^ZdnQ''{{K+wj)T) x(iQ<^{{K+wj)T)-Q{{K+wjT-a^) 

< |Qe|e^4.2card(aeZ'^ n Q%{K + wj)T) + Cde'^^HK + wj)T]''~^a, 

< c,|Q,|eT4.2a,-'^((ir + wj)Tf + c^e-^K + wj)TY-\, 

< Cd{K + wJifAi2{log{e-^)y[e^i-2-'^ + e(Vi6d)-d] 

< e74.2/2-<i^ (7,32) 

for small enough e, where we have used 74.2 < (Wd)"^ in the last line. We have 
shown that for all e smaller than some positive £0, 

P{e,T{Q"(iK+wj)T)) < e'>4.2/2-<i for J = 1, ... , J,- Jo) > l-(Ji- Jo)e-°^ (7.33) 
Step 3. Fix j G {1, . . . , Ji — Jq}, and define {Ci'^,t > jT) by setting 

■''^ I otherwise, 



and then using our Poisson processes {A^,A^ : x G eZ } to continue constructing 
(^^'^ in the same way as |f is constructed. By Lemma 7.6, if e < £7.6(74.2/2, ??2), 

ejTiQ'HK + iw^)T) < e(T4.2/2)-rf implies 

Pdl^'^l > 0|^|r) < 6e2^''e^4.2/4 for t > {j + - 1)T. (7.34) 

Using we will show that with high probability, 

^jTiQ'iiK + jw)T) < e(T4.2/2)-rf implies 

Q%iK + jw - JoA)T) X [(j - 1 + Jo)T, (j + Jo)T] is e-empty. (7.35) 

To do this, define the event 

^tU) = n^li^Tij + i,K + wj- iA, K + wj-{i- l)A). 

Using Lemma 4.7 as in Step 1 we have for small enough e 

PiTrU + i,K + wj- iA, K + wj - {i - 1)A)) > 1 - ca{K + wj - lAfe'^ 

and thus 

P{Tt{3)) > 1 - c'^iSKfJoe'' (7.36) 

for small enough e. 
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Observe that on the event Txij) we have 

Xt'"" c Q'iiK + jw)T) V X G Q'iiK + jw - JoA)T), 

ue[{j-l + Jo)T, (j + Jo)T], and te[0,u- jT]. (7.37) 

Therefore, by duahty, on TT{j), 

etix) = i'/{x) for ah 

(x, t) G Q'{{K + jt^ - JoA)T) X [(j - 1 + Jo)T, (j + Jo)r]. 

Combining this observation with (7.34) and (7.36) we see that the event in (7.35) 
has probabihty at least 

1 - 6e2'=''eT4.2/4 _ c'a{3K)^Joe'^ 

for e smaller than some Eq. 

Step 4 We can now sum the last estimate over j = 1, . . . , Ji — Jq and use (7.33) 
to obtain 

P{Vi is e-empty) > 1 - Ji(e-°^ + 6e2'^''eT4.2/4 + c'a{3K)'^ Jqe"^) (7.38) 

for small enough e. (Actually we get a slightly larger set than Pi.) Combine (7.29) 
and (7.38) to obtain (7.25). 

Step 5 Finally, using the notation from Step 1, the event in (7.26) is just 

ri-UTrU, K + wJi + ( Ji - j)A, K + wJi + ( Ji - j + 1)A). 

As in Step 1, we can use Lemma 4.7 to bound the probability of this intersection 
by 1 — c!^J\{K + Ji{w + A)Ye'^ for small enough e, so we are done. □ 

Let Ki = K + Ji{w + A). For ^ G {0, 1}^^'', let G| be the event, depending 
on our graphical representation, on which V = 'D{Jq, Ji,w, A, K) is e-empty if 
= ^, and on which (7.26) holds. Note that (7.26) implies all the duals starting at 
(x,n) G P and run up until time u remain in Q(KiT). Hence duality implies that 
G| is Q{Q{KiT) X [0, JiT])-measurable, where we recall from (6.7) that Q{R) is the 
cT-field generated by the Poisson points in the graphical representation in the region 
R. By the inclusion (7.23) we have 

on G|, Q^{2T) x [0, JiT] is e - empty, and CjMQ'i^KT)) = 0, (7.39) 

providing that £ H = e {0,1}^^'' : i{Q^{KT)) = 0}. Adding the bounds in 
Lemma 7.7 we see that 

if ^ G if, then P{G\) > 1 - 2c7.7e-°^^^ if e < £7.7- (7.40) 
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7.3 Proof of Theorem 1.16 



Proof of Theorem 1.16. We continue to take T = A4.2 log(l/e), and with Ji from 
(7.22) we define 

L = T, T' = JiT, 

and set cl = L/{2D) as before. We set ^^{y) = l{\y\ > L),y € eZ'^, and cj^, z € eZ'^ 
denote the translation operators on {0, 1}^^ . For (x,n) Cd let 

1 otherwise, 

and define the percolation variables 

r]{x,n) = l(G|i,n occurs in the graphical representation in which the Poisson 

processes are translated by —clx in space and —riT' in time). (7.41) 

In the percolation argument which follows it is the first part of the definition of 
that will matter; the is really only a place-holder which allows us to define 
r] when the translated configuration is not in H. As in the proof of Theorem 1.15 
in Section 6, we are actually translating in space by the "lower left hand corner" in 
e'L'^ associated with —clx and as before suppress this in our notation. In Section 6 
we used Theorem 4.3 of [14]; here we copy the key definition in its proof. Using 
the measurability of G|, the independence of Q{R) for disjoint regions R, and (7.40) 
one can check that for any any M > ADKi, the family {ri{z), z £ Cd} satisfies the 

"''4 2 

modified M-dependent condition (5.5) with d = 207,78'^^^^^ . To see this argue 
exactly as in the proof of Theorem A. 4 of [14]. 

Using the percolation results from Section 5, we will show 

Lemma 7.8. There exists f > such that for e small enough ^q{Q^{KT)) = 
implies 

P(A) = P{il{Q^{ft)) = for all large t) > 1/2. 

Intuitively this is an immediate consequence of Lemmas 5.1 and 5.5. The first 
result implies that on then for large n, the wet sites satisfy W°n'H™ = WnCiTiJ^. 
The second result shows that if Bn is the collection of dry sites in T-CT^^ connected 
to the complement of U" _ in.'H"^^'^ by a path of dry sites on the graph with edge 
set then Bn = 9 eventually. Wet sites in V.^ will correspond to space-time 
blocks that are empty of I's while dry sites (i.e. not wet sites) in T-Ll[^ correspond 
to space-time blocks which may contain a 1. If a dry site in TilT^^ corresponds to 
a block containing a 1 there must be a dual path of I's leading from this 1 to a 
site outside of ^m=n/2^"^^'^ • This corresponds to a path of dry sites in and so 
cannot happen for large n since Bn = 9 for large n. Thinking of the corresponding 
space time regions are being filled with concrete, and the dry sites as air spaces, 
we see that there cannot be a 1 in T-CrT^'^ unless some air space reaches outside of 
U" _„ /o^JtT''^- We now give a formal proof. 

Tit — Th / Z 
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Proof. Fix as above and recall W^,Wn,Bn and Vr^ from Section 5. In particu- 
lar M^'^jTy are constructed from an iid Bernoulli field which is bounded above by 
r]{z,n), z £ T-ln,n > 1. By (7.40) and our condition on (which implies ^^'^ = 
in the definition of r]{0,0) = 1) we see that 

P(7?(0,0) = 1) > 1 - 2c7.7e-°^''(^4.2/4) > 3/4, 

for e small enough. By working with P{-\rj{0,0) = 1) in place of P we may assume 
77(0, 0) = 1 at a cost of proving (under our new P) that 

P(A) > 3/4. (7.42) 

Assume n > 1 and (y, n) € W^. Then for some i, letting y' = y — vl, {y' , n — 1) G 
W^n-i with 7]{y',n — 1) = 1 (if n = 1 we use ?7(0, 0) = 1 here). Continue to trace 
back the set of open sites y' = y'n-i-, ■ ■ ■ ,yo = 0. Proceeding through the y'^ values, 
using the second part of (7.39) and clv- + [-L, if C [-2L, 2Lf for i = 1, . . . , D, 
we see that ^^f^ = a_c^y'XCfrp,) in the definition of r]{y^,i) = 1. Therefore (7.39) and 
translation invariance, show that ri{y',n — 1) = 1 implies 

CticLy' + Q'{2L)) = for ah t G [(n - l)T',nT']. 

Since CLy + Q^{L) C ciy' + Q^{2L) we obtain 

{y,n) G implies (cl^ + Q%L)) = for ah t € [{n - l)T',nT']. (7.43) 

This confirms (5.8) in Section 5. 

Next by Lemma 5.1 we may assume e is small enough (independent of the choice 
of Co) so that P(O^) > 3/4 and 9 < 65,1 A 6*5.4. Let r = r5.i and assume w G 
By Lemma 5.1 there is an no G N so that 

Wfc° n ni'' = Wk r\'Hf'ik>no. (7.44) 

Let f = and assume uj ^ A. The latter implies that for infinitely many n > 2no 

there are t G [{n — l)T',nT'] and x G Q^{rt) with ^j(x) = 1. We claim that this 
implies 

Bn ^ 9 for n as above. (7.45) 

Lemma 5.5 implies the above is a null set, so it follows that Pi^'^ \ A) = and so 
(7.42) would be proved (recah P{nl^) > 3/4). 

To prove (7.45) fix such an n and trace backward in time a path of I's that 
leads to ^ti^) — 1- Here we are using (1.65) to show if all the inputs are the dual 
process will produce a at a given site. By (5.2) there must exist some (y, n) G T-Ln 
such that X G ciVy C c^y + Q'^{L) and a bit of arithmetic using the definition of f 
gives 

\x\ L fnT' ^ rn 

\y\ <— + — < + 2D<—, 

CL CL CL 4 

and we have taken no big enough for the last inequality. Hence (y, n) G T-Cn^^ and 
so (7.43) and (7.44) imply (y,n) ^ Wn, i.e., (y,n) is dry. By duality and the finite 
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range assumption (recall (7.21)), there must exist x' € eL^ and t' G \{n — l)T',t) 
such that \x — x'\ < RqE and £,fi{x') = 1. That is, t' > {n — 1)T' is the first time 
below t that the dual jumps or t' = {n — 1)T' if there is no such time in which 
case x' = X. We may assume e is small enough so that Rqe/cl < C5.3, in which 
case by (5.3) x' G c^Vy/ for some y' of the form y + v[ — Vj {y = y' is included). 

If {y',n) € Tin"'^'^ C H™, it follows from (7.43) that {y',n) must be dry, and thus 
(y',n) G £'(j;,„). 

Continue the above construction until either we reach a point {y" ,n) G {'HlT^'^Y 
with all earlier points in our path from {y,n) being dry, or we obtain x",y" such 

that Cln-i)T'^^"^ ~ (y"'"'^) ^ ^{y,n) 'HIT^'^ and x" G c^Vy". In the former 
case Bn ^ (recall the precise definition prior to Lemma 5.5). In the latter case if 
{y" — v[^n — 1) ^ T'C^-i ^''^^ for some i, then (7.45) holds. If not, then as one easily 
checks \cL{y" — v[) — x"\ < L, and so arguing as above, we see that {y" — v[,n — 1) 
is dry. Therefore the iteration can be continued until it stops as above or continues 
down to time (^ — 1)T', again forcing (7.45) in either case. □ 

Having established Lemma 7.8 the rest of the proof of Theorem 1.16 is routine. 
The proof of Lemma 6.5 shows that if we start from an initial configuration with 
infinitely many O's then at time 1 there will be infinitely many cubes of the form 
clx + Q'^{L) with X G Hq that are e-empty. By the Markov property this will hold at 
all times G N a.s. The above shows that if xq is chosen so that ^\{xq + Q{L)) = 1, 
then w. p. at least 1/2, = on clXq + Q{ft) for all large t. If this fails at some 
time we can try again at a later time N by the above and after a geometric (1/2) 
number of trials we will succeed and produce a linearly growing set of O's starting 
at some space-time location. Therefore the O's take over. □ 
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